nr  pi c  rqpy  ADA115168 


AFWL-TR-79-72 


AFWL-TR- 

79-72 


ENERGETICS  AND  RADIATIVE  PROPERTIES 
OF  EXCIMER  LASERS 

H.  H.  Michels 
R.  H.  Hobbs 

United  Technologies  Research  Center 
Silver  Lane 

East  Hartford,  CT  06108 


1/ 


1/ 


January  1981 


Final  Report 


Approved  for  public  release;  distribution  unlimited. 


DTIC 

ELECTEf! 
JUN  7  1982 


J 


AIR  FORCE  WEAPONS  LABORATORY 
Air  Force  Systems  Command 
Kirtland  Air  Force  Base,  NM  87117 


82  06  07  iOi 


AFWL-TR-79-72 


SlL“aS.PreP2red,by  ‘he  Unltetl  Technologies  Research  Center, 
“51  KrtT°rd*  Connectlcut,  under  Contract  F29601 -78-C-0031 ,  Job  Order  34320801 
with  the  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base,  New  Mexico.08 
Cap  am  Lawrence  A.  Wright  (NTY)  was  the  Laboratory  Project  Officer-in-Charge. 

When  US  Government  drawings,  specifications,  or  other  data  are  used  for  anv 
purpose  other  than  a  definitely  related  Government  procurement  operation  the 
Sje?SnLtJeSyrinCUrS  ^  ^Ponsibimy  nor  anyobl  igatiSn  whatsoe^;,  and 
suDofied  5 °vernment  may  have  formulated,  furnished,  or  in  any  way 

supp  led  the  said  drawings,  specifications,  or  other  data,  is  not  to  be  reoarded 

npr!!n  1Catl°n  or  otherwise,  as  in  any  manner  licensing  the  holder  or  any  other 

o'r  sen  arnCv0oTtprn  :H°nin0r  ?nve£nJ  any  ”9hts  or  Permission  to  man^a^re  use, 
y  patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  has  been  authored  by  a  contractor  of  the  United  States 
overnment.  Accordingly,  the  United  States  Government  retains  a  nonexclusive 

alfn^nthp6!  to  Publish  or  reproduce  the  material  contained  herein,  or 

allow  others  to  do  so,  for  the  United  States  Government  purposes. 

m  Iepor?  been  reviewed  by  the  Public  Affairs  Office  and  is  releasable 

to  the  National  Technical  Information  Service  (NTIS).  At  NTIS,  it  will  be 
available  to  the  general  public,  including  foreign  nations 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


|  \  .  ..  "\  ' ,  t  JL 

vv C( tYI 

LAWRENCE  A.  WRIGHT 
Captain,  USAF 
Project  Officer 


Chief,  Advanced  Concepts  Branch 


FOR  THE  DIRECTOR 

**  ]/)/ t 

THOMAS  W.  CIAMBRONE 
Colonel,  USAF 
Applied  Physics  Division 


DO  NOT  RETURN  THIS  COPY.  RETAIN  OR  DESTROY. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  of  THIS  page  (When  Pete  En(»f ed) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  U.  GOVT  ACCESSION  NO. 

AFWL-TR-79-72  ttfrA/LC  K08 

3.  RECIPIENT’S  CATALOG  NUMBER 

A.  TITLE  (mnd  Subtitle) 

ENERGETICS  AND  RADIATIVE  PROPERTIES  OF  EXCIMER 
LASERS 

5.  TYPE  OF  REPORT  A  PERIOO  COVEREO 

Final  Report 

6.  PERFORMING  ORG.  REPORT  NUMBER 

924270 

7.  AUTHORf*; 

H.  H.  Michels 

R.  H.  Hobbs 

8.  CONTRACT  OR  GRANT  NUMBER^*.) 

F29601  -78-00031 

9  performing  organization  name  ano  aooress 

United  Technologies  Research  Center 

Silver  Lane 

East  Hartford,  CT  06108 

10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  6  WORK  UNIT  NUMBERS 

62301  E/34320801 

tl.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Air  Force  Weapons  Laboratory  (NTY) 

Kirtland  Air  Force  Base,  NM  87117 

12.  report  oate 

January  1981 

13.  number  of  pages 

168 

14.  MONITORING  AGENCY  NAME  4  AOORESSf//  dillerent  Irom  Controlling  Ottice) 

15.  SECURITY  CLASS,  (of  thi»  rmport ) 

Unclassified 

IS*.  OECLASSI  FI  CATION/  OOWN  GRADING 
SCHEOULE 

16.  OISTRIBUTION  ST  ATEMENT  (ol  thie  /import) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (ol  the  ebetrmct  entered  in  Block  JO,  il  ditlotenl  tram  Report ) 

16.  supplementary  NOTES 

This  research  was  sponsored  by  the  Defense  Advance  Research  Projects  Agency 
(ARPA)  under  ARPA  Order  Number  3432. 

19.  KEY  WOROS  (Continue  on  reveree  tide  II  nocoooory  ond  identity  by  block  number) 

ArF  Ar2+  Ar^"1"  Excimer  Molecules 

Ar?F  j.  +  Noble  Gas  Dimer  Ions 

+  %  Electronic  Structure  Noble  Gas  Trimer  Ion5 

Ne?  Xe2 

Z(^  ABSTRACT  (Continue  on  reveree  etde  it  neceeeery  end  Identity  by  block  nuirlwr) 

This  research  program  was  concerned  with  a  theoretical  study  of  the  energetics 
of  the  noble  gas  halides  and  noble  gas  dimer  and  trimer  molecular  ions.  The 
electronic  structure,  potential  energy  curves,  and  radiative  transition 
probabilities  of  excimer  systems  such  as  ArF  and  Ar2F  have  been  examined 
using  quantum  mechanical  methods. 

A  systematic  .study  of  the.  electronic  structure  in  the  noble  gas  dimer  ion 
sequence,  Ne2+,  Ar2+,  Kr2+  and  Xe2+,  was  carried  out.  This  study  included 

DD  ,j?nMtj  1473  _ UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  'Wh en  Pete  Entered) 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGEQVhfi  Jaf  Enfrmd) _ 

20.  ABSTRACT  (Continued) 

J/ 

detailed  calculations  of  the  pertinent  potential  energy  curves  and  an 
analysis  of  the  calculated  spectroscopic  properties  of  the  bound  states 
of  these  ions. 

Finally,  calculations  were  carried  out  on  the  electronic  structure  and 
photoabsorption  characteristics  of  the  noble  gas  trimer  ’on  Ar3+.  This 
study  included  detailed  calculations  of  the  potential  energy  surfaces 
for  Arj-  ;  an  analysis  of  the  spectroscopic  properties  of  the  ground  state 
of  this  ion  and  prediction  of  the  strong  absorption  bands. 

T 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAG ZfWhvt  Dttm  Enfn4) 


SUMMARY 


The  electronic  structure,  potential  energy  curves,  and  radiative  tran¬ 
sition  probabilities  of  excimer  systems  have  been  examined  using  quantum 
mechanical  methods.  These  molecules  are  characterized  by  repulsive  or 
weakly  bound  ground  state  potential  curves  and  by  bound,  strongly  ionic,  or 
Rydberg,  excited  states.  They  constitute  a  very  interesting  class  of 
molecules  which  offer  the  possibility  for  high  power,  high  efficiency  UV 
laser  operation.  Calculations  have  been  carried  out  using  the  density 
functional  SCF-Xa  method,  modified  extensively  to  correct  for  well  known 
errors  arising  from  approximations  to  the  potential  and  exchange  terms. 

A  limited  number  of-  ab  in it io  calculations  were  also  carried  out  for  com¬ 
parison  purposes.  For  the  ArF  system  we  find  that  the  lowest  excited  ionic 
state  has  symmetry  ^1/2 >  and  that  the  dominant  laser  transition  observed 

at  1933  X  should  be  assigned  to  B  2^  2j:1/2-  The  C  ^1/2  tran“ 

sition  is  calculated  to  be  two  orders  of  magnitude  smaller  in  emission  inten¬ 
sity  than  the  dominant  transition,  thus  ruling  out  this  assignment  for  the 
observed  laser  line  in  ArF.  Preliminary  calculations  carried  out  for  Ar2F 
indicate  that  the  bound  upper  ionic  state  has  B£  symmetry  and  that  the 
most  probable  ground  state  also  has  “B2  symmetry.  This  polyatomic  system 
is  predicted  to  exhibit  a  broadband  emission  spectrum  with  a  relatively 
long  radiative  lifetime. 
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In  addition  to  calculations  performed  on  noble  gas  halide  excimer 


molecules,  a  systematic  study  of  the  electronic  structure  and  chemical 
binding  in  the  dimer  ion  seai't-ice,  Ne2»  A^,  K^,  and  Xe2»  has  been  carried 
out  using  density  functional  methods.  For  comparison,  ab  initio  configuration 
interaction  calculations  were  also  performed  for  the  A^  ion.  These  studies 
include  detailed  calculations  of  the  pertinent  potential  energy  curves  and 
an  analysis  of  the  calculated  spectroscopic  properties  of  the  bound  states 
of  these  ions.  A  regular  progression  is  found  in  the  spectroscopic  properties 
for  the  ground  A  state  which  leads  to  some  remarkably  simple  con¬ 

clusions  concerning  the  nature  of  the  binding  and  the  size  of  these  dimer 
ions.  For  the  heavier  systems,  Kr*  and  Xe*,  spin-orbit  coupling  becomes 
important,  resulting  in  a  strong  mixture  of  the  A-S  coupled  E  and  II  states. 
This  mixing  affects  the  strength  of  the  binding  in  the  ground  state.  A 
comparison  with  other  ab  initio  studies  and  analsyis  of  the  asymptotic  be¬ 
havior  at  large  internuclear  separations  is  given. 

As  a  final  study  under  this  research  program,  calculations  were  carried 
out  on  the  electronic  structure  and  photoabsorption  characteristics  of  the 
noble  gas  trimer  ion  Ar*-  This  study  included  detailed  density  functional 
calculations  of  the  potential  energy  surfaces  for  Ar*,  an  analysis  of  the 
spectroscopic  properties  of  the  ground  state  of  this  ion  and  predicition  of 

the  strong  absorption  bands.  Studies  in  D^,  ^2v’  an<^  ^°°h  sXmmetr^es  were 

+  2  ' 
carried  out.  The  lowest  state  of  Ar^  corresponds  to  the  degenerate  E 

symmetry  with  an  indicated  dissociation  energy  of ^0.17  eV  relative  to 
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+  2  + 

Ar2  (A  Su  )  +  Ar  Clg0) •  This  trimer  ion  exhibits  a  very  small  Jahn-Teller 
distortion  which  does  not  significantly  alter  our  conclusions  concerning 
the  ground  state  stability  or  absorption  characteristics.  A  strong  and  very 

2  t  2  ' 

broad  photoabsorption  band  (1=426-670  nm)  is  predicted  for  the  E  -*■  ~A^ 
transition  of  this  trimer  ion.  This  suggests  that  such  species  may  be  very 
important  in  the  analysis  of  loss  mechanisms  in  visible  excimer  lasers, 
operating  under  high  pressure  conditions  with  mixtures  containing  the  noble 
gases . 
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SECTION  I 


INTRODUCTION 

The  noble  gas  monohalides  constitute  an  interesting  and  important  class 
of  molecules  for  laser  applications.  These  molecules  are  characterized  by 
a  repulsive  or  very  weakly  bound  ground  state  potential  curve  and  bound, 
strongly  ionic,  excited  states.  The  character  of  the  electronic  structure 
of  these  excited  states  is  similar  to  that  of  the  ground  state  alkali  halide 

salts . 

Emission  spectra  from  several  noble  gas  halides  have  been  obtained  in 
low  pressure  discharge  flow  experiments  and  were  first  reported  by  Velazco 
and  Setser  (ref.  1)  and  by  Golde  and  Thrush  (ref.  2).  More  recently  emission 
spectra  for  several  systems  which  have  exhibited  lasing  have  been  reported 
by  Brau  and  Ewing  (ref.  3).  These  spectra  are  consistent  with  a  model  where¬ 
by  the  lasing  transition  is  between  a  strongly  bound  upper  electronic  state 
and  a  nearly  flat  repulsive  ground  state.  A  group  theoretical  analysis  of 

the  possible  electronic  configurations  for  these  states  does  not  lead  to  a 

2  + 

unique  assignment  since  both  the  upper  and  lower  levels  could  be  of  I  or 
^ H  symmetry.  Even  a  detailed  experimental  spectroscopic  study  of  XeCl 
(ref.  4)  failed  to  uniquely  assign  the  symmetries  of  the  lasing  transition. 

Theoretical  approaches  to  this  problem  have  consisted  of  analogies  with 
alkali  halide  systems  (ref.  5),  predictions  based  on  an  electron-gas  model 
(ref.  6)  and  ab  initio  calculations  for  the  KrF  system  (ref.  7).  Unfortu¬ 
nately,  the  theoretical  methods  utilized  to  date  have  been  unable  to 


positively  define  the  electronic  states  involved  in  the  noble  gas  halide  laser 


systems.  The  ab  initio  methods  would  appear  to  be  applicable  only  to  the 
lighter  systems  such  as  ArF  and  KrF  owing  to  greatly  increased  computational 
complexity  for  larger  numbers  of  electrons.  In  addition,  past  experience 
indicates  that  the  _ab^  initio  calculations  would  not  be  of  uniform  quality  in 
going  from  a  simple  system  such  as  ArF  to  a  heavier  system  such  as  XeCl.  The 
electron-gas  model  calculations  suffer  from  the  lack  of  self-consistency  in 
the  molecular  potential.  Recent  density  functional  calculations  performed 
in  this  Center  indicate  self-consistency  in  the  potential  is  crucial  to 
determining  optimum  electronic  wave func t ions  for  these  systems.  Thus  a 
theoretical  study  of  the  electronic  structure  of  the  noble  gas  halide  systems, 
capable  of  uniformly  describing  XeBr  as  well  as  ArF,  and  relevant  to  the 
understanding  and  interpretation  of  current  experimental  programs  in  noble 
gas  halide  laser  systems  seems  desirable  at  the  present  time. 

The  primary  information  to  be  obtained  from  a  theoretical  study  of 
noble  gas  halide  systems  is  the  electronic  structure,  energetics  and  life¬ 
times  of  the  low-lying  excited  states  which  are  involved  in  the  chemical 
reaction  mechanisms  leading  to  population  of  the  excited  laser  levels. 
Knowledge  of  the  partition  of  energy  between  the  products  of  a  chemical 
reaction  is  fundamental  to  the  understanding  of  the  dynamics  of  the  process, 
and  to  the  design  of  chemical  lasers.  Atom-atom  and  atom-molecule  reactions 

that  yield  electronically  excited  molecules  are  of  particular  scientific  and 
military  interest.  From  the  point  of  view  of  the  theoretical  chemist,  an 
understanding  of  the  dynamics  of  the  simplest  of  these  reactions  is  a  pre¬ 
requisite  for  understanding  more  complicated  gas  phase  reactions.  From  the 
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viewpoint  of  the  engineer  interested  in  the  design  of  chemical  lasers, 
studies  of  these  reactions  yield  kinetics  information  and  scaling  parameters 
that  enable  one  to  design  devices  with  the  potential  of  having  very  high 
specific  power.  A  great  deal  of  work  in  these  two  areas  of  science  and 
engineering  has  been  stimulated  by  the  recent  development  of  a  number  of 
highly  sophisticated  techniques  for  measuring  the  nonequilibrium  distribu¬ 
tion  of  the  vibrational  and  rotational  states  of  the  products  of  a  chemical 
reaction.  Over  the  past  decade,  many  atom-molecule  reactions  yielding 
vibrationally  excited  hydrogen  halides  have  been  studied  by  infrared  chemi¬ 
luminescent  techniques.  These  studies  have  furnished  a  great  deal  of  informa¬ 
tion  on  how  the  exothermicity  of  a  chemical  reaction  may  be  distributed  among 
the  energy  levels  (vibration,  rotation,  and  translation)  of  the  products 
(refs.  8-12).  More  recently,  it  has  been  possible  to  use  laser  technology 
and  molecular  beam  techniques  to  measure  the  internal  energy  distribution  of 
reactions  yielding  vibrationally  excited  hydrogen  halides  (refs.  13-15). 

A  knowledge  of  at  least  limited  regions  of  potential  energy  surfaces  is 
of  importance  in  analyzing  the  electronic,  vibrational  and  rotational  relaxa¬ 
tion  rates  in  chemical  laser  systems.  These  relative  relaxation  rates  govern 
both  the  feasibility  of  lasing  and  the  power  that  is  potentially  available. 
Such  relaxation  rates  are  governed  by  long-range  forces  between  atomic  and 
molecular  fragments  and  thus  knowledge  of  the  dissociation  behavior  repre¬ 
sented  by  the  potential  energy  surfaces  for  the  laser  system  can  be  used  in 
calculations  of  the  kinetic  behavior.  A  recent  example  is  the  use  of  such 
potential  energy  data  in  analyzing  deactivation  in  CC>2  (OO'l)  through 


The  discussion  above  has  indicated  the  way  in  which  knowledge  of  poten¬ 
tial  energy  surfaces  can  be  useful  in  analyzing  the  behavior  of  an  existing 
chemical  laser  system.  More  recently  it  has  become  apparent  that  detailed 
knowledge  of  such  surfaces  can  be  utilized  in  a  predictive  fashion.  Michels 
and  Broida  (ref.  17)  have  recently  outlined  a  mechanism  for  the  chemilumines¬ 
cence  in  the  visible  region  of  the  BaO  system.  This  suggestion  was  based  on 
the  combined  knowledge  of  the  results  of  calculations  on  the  potential  energy 
curves  for  BaO  and  measured  lifetime  data  for  excited  electronic  states  of 
this  system.  A  similar  parallel  study  of  A10  has  also  been  reported  by 
Michels  (ref.  18)  and  Johnson,  et  al.  (ref.  19). 

The  initial  impetus  for  these  studies  came  from  a  practical  problem  in 
aeronomy.  The  release  of  certain  chemical  species  into  the  upper  atmosphere 
results  in  luminous  clouds  that  display  the  resonance  electronic-vibrational- 
rotation  spectrum  of  the  released  species.  Such  spectra  are  seen  in  rocket 
releases  of  chemicals  for  upper  atmospheric  studies  and  upon  reentry  into  the 
atmosphere  of  artificial  satellites  and  missiles.  Of  particular  interest  in 
this  connection  is  the  observed  spectra  of  certain  metallic  oxides  and  halides. 
From  band  intensity  distribution  of  the  spectra,  and  knowledge  of  the  f-values 
for  electronic  and  vibrational  transitions,  the  local  optical  conditions  of 
the  atmosphere  can  be  determined. 

Because  of  inherent  difficulties  in  the  experimental  determination  of 
transition  probabilities  for  certain  metal  oxide  and  halide  systems  and 
because  of  recent  progress  in  the  _ab_  initio  calculation  of  electronic  wave- 
functions  for  diatomic  systems,  theoretical  programs  for  calculating  band 
absorption  and  emission  coefficients  were  undertaken  at  UTRC  (Contracts 
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F29b01— 71— C— 0 1 19  and  F44620-73-C-0077) .  These  programs  based  on  current 


quantum  mechanical  techniques  and  capabilities  developed  at  UTRC  were  direc¬ 
ted  toward  determining  the  vibrational  and  electronic  transition  probabilities 
and  integrated  absorption  coefficients  for  the  metal  oxide  band  systems  SiO, 
A10,  TiO,  FeO,  BaO  and  UO.  Additional  studies  of  halide  systems  have  included 
LiCl,  LiF,  NaF ,  HC1  and  HF .  Progress  on  the  studies  of  certain  of  these  sys¬ 
tems  has  been  reported  (refs.  16-21).  These  studies  formed  a  very  sound 
basis  for  more  extensive  analysis  of  metal  oxides  and  metal  halides  as  candi¬ 
date  laser  systems.  To  the  best  of  our  knowledge,  no  other  laboratory  in  the 
country  has  such  an  extensive  catalogue  of  theoretical  calculations  of  the 
electronic  structure  of  diatomic  oxides  and  halides  or  our  computational 
experience  with  these  types  of  molecules.  In  addition,  UTRC  has  fully  devel¬ 
oped  sophisticated  computer  programs  for  studying  both  diatomic  and  poly¬ 
atomic  molecules.  These  programs  have  been  extensively  used  at  UTRC  and  at 
the  Air  Force  computational  laboratories  at  AFWL,  Kirtland  Air  Force  Base 
and  AFGL,  Hanscom  Field. 

The  research  program  presented  herein  was  devoted  to  a  theoretical  study 
of  the  energetics  of  noble  gas  halides,  the  prediction  of  the  radiative  life¬ 
times  of  electronically  excited  species,  and  studies  of  photon  energy  loss 
paths  such  as  photoabsorption  or  ionization  of  the  upper  excited  electronic 
states.  The  goal  of  this  research  program  was  to  develop  sufficient  informa¬ 
tion  to  enable  one  to  determine  promising  candidates  for  UV  laser  systems  and 
to  develop  an  information  base  of  kinetic  data  useful  for  analyzing  loss 
process  in  excimer  lasers. 
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SECTION  II 


CURRENT  STATUS  OF  QUANTUM  MECHANICAL 
METHODS  FOR  DIATOMIC  SYSTEMS 

The  application  of  quantum  mechanical  methods  to  the  prediction  of 
electronic  structure  has  yielded  much  detailed  information  about  atomic  and 
molecular  properties  (ref.  22).  Particularly  in  the  past  few  years,  the  avail¬ 
ability  of  high-speed  computers  with  large  storage  capacities  has  made  it 
possible  to  examine  both  atomic  and  molecular  systems  using  an  ab  init io 
approach,  wherein  no  empirical  parameters  are  employed  (ref.  23).  Ab  initio 
calculations  for  diatomic  molecules  employ  a  Hamiltonian  based  on  the  nonrela- 
tivistic  electrostatic  interaction  of  the  nuclei  and  electrons,  and  a  wave- 
function  formed  by  antisymmetrizing  a  suitable  many-electron  function  of 
spatial  and  spin  coordinates.  For  most  applications  it  is  also  necessary  that 
the  wavefunction  represent  a  particular  spin  eigenstate  and  that  it  have 
appropriated  geometrical  symmetry.  Nearly  all  the  calculations  performed 
to  date  are  based  on  the  use  of  one-electron  orbitals  and  are  of  two  types: 
Hartree-Fock  or  configuration  interaction  (ref.  24). 

Hartree-Fock  calculations  are  based  on  a  single  assignment  of  electrons  to 
spatial  orbitals,  following  which  the  spatial  orbitals  are  optimized,  usually 
subject  to  certain  restrictions.  Almost  all  Hartree-Fock  calculations  have  been 
subject  to  the  assumption  that  the  diatomic  spatial  orbitals  are  all  doubly 
occupied,  as  nearly  as  possible,  and  are  all  of  definite  geometrical  symmetry. 
These  restrictions  define  the  conventional ,  or  restricted ,  Hartree-Fock  (RHF) 
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method  (refs.  25  and  26).  RHF  calculations  can  be  made  with  relatively  large 


Slater-type  orbital  ( STO)  basis  sets  for  diatomic  molecules  with  first  or 
second-row  atoms,  and  the  results  are  convergent  in  the  sense  that  they  are 
insensitive  to  basis  enlargement.  The  RHF  model  is  adequate  to  give  a  qualita¬ 
tively  correct  description  of  the  electron  interaction  in  many  systems,  and 
in  favorable  cases  can  yield  equilibrium  interatomic  separations  and  force 
constants.  However,  the  double-occupancy  restriction  makes  the  RHF  method 
inappropriate  in  a  number  of  circumstances  of  practical  interest.  In  particu¬ 
lar,  it  cannot  provide  potential  curves  for  molecules  dissociating  into 
odd-electron  atoms  (e.g.,  NO  at  large  intemuclear  separation),  or  into  atoms 

3  _  3 

having  less  electron  pairing  than  the  original  molecule  [e.g.,  O2  ♦  0(JP)]; 

3 

it  cannot  handle  excited  states  having  unpaired  electrons  (e.g.,  the  Z 
states  of  O2  responsible  for  the  Schuraann-Runge  bands);  and,  in  general,  it 
gives  misleading  results  for  molecules  in  which  the  extent  of  electron  correla¬ 
tion  changes  with  internuclear  separation. 

Configuration-interaction  (Cl)  methods  have  the  capability  of  avoiding 
the  limitations  of  the  RHF  calculations.  If  configurations  not  restricted  to 
doubly-occupied  orbitals  are  included,  a  Cl  can,  in  principle,  converge  to  an 
exact  wavefunction  for  the  customary  Hamiltonian.  However,  many  Cl  calculations 
have  in  fact  been  based  on  a  restriction  to  doubly-occupied  orbitals  and 
therefore  retain  many  of  the  disadvantages  of  the  RHF  method  (ref.  24).  The 
use  of  general  Cl  formulations  involves  three  considerations,  all  of  which 
have  been  satisfactorily  investigated:  the  choice  of  basis  orbitals,  the 
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choice  of  configurations  (sets  of  orbital  assignments),  and  the  specific  cal¬ 
culations  needed  to  make  wavefunctions  describing  pure  spin  states  (ref.  27). 
The  first  consideration  is  the  art  associated  with  quantum  mechanical  elec¬ 
tronic  electronic  structure  calculations.  Many  methods  (iterative  NSO, 
perturbation  selection,  first  order  Cl,  etc.)  have  been  advocated  for  the 
optimum  choice  of  configurations.  There  are  no  firm  rules  at  present  and  the 
optimum  choice  is  a  strong  function  of  the  insight  of  the  particular  research 
investigator.  The  last  consideration,  proper  spin  and  symmetry  projection, 
has  proved  difficult  to  implement,  but  computer  programs  have  been  developed 
for  linear  projection  algebra  at  this  Center,  and  the  Cl  method  has  been  found 
of  demonstrable  value  in  handling  excited  states  and  dissociation  processes 
which  cannot  be  treated  with  RHF  techniques. 

Either  of  the  above  described  methods  for  ab  initio  calculations  reduces 
in  practice  to  a  series  of  steps,  the  most  important  of  which  are  the  evalua¬ 
tion  of  molecular  integrals,  the  construction  of  matrix  elements  of  the 
Hamiltonian,  and  the  optimization  of  molecular  orbitals  (RHF)  or  configuration 
coefficients  (Cl).  For  diatomic  molecules,  these  steps  are  all  comparable  in 
their  computing  time,  so  that  a  point  has  been  reached  where  there  is  no 
longer  any  one  bottleneck  determining  computation  speed.  In  short,  the  inte¬ 
gral  evaluation  involves  the  use  of  ellipsoidal  coordinates  and  the  introduc¬ 
tion  of  the  Neumann  expansion  for  the  interelectronic  repulsion  potential 
(ref.  28);  the  matrix  element  construction  depends  upon  an  analysis  of  the 
algebra  of  spin  eigenfunctions  (ref.  29);  and  the  orbital  or  configuration 
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optimization  can  be  carried  out  by  eigenvalue  techniques  (refs.  30  and  31). 

All  the  steps  have  by  now  become  relatively  standard,  and  can  be  performed 
efficiently  on  a  computer  having  32,000  to  65,000  words  of  core  storage,  a 
cycle  time  in  the  microsecond  range,  and  several  hundred  thousand  words  of 
peripheral  storage. 

Both  the  RHF  and  Cl  methods  yield  electronic  wavefunctions  and  energies 
as  a  function  of  the  interauclear  separation,  the  RHF  method  for  one  state, 
and  the  Cl  method  for  all  states  considered.  -.The  electronic  energies  can  be 
regarded  as  potential  curves,  from  which  may  be  deduced  equilibrium  inter- 
nuclear  separations,  dissociation  energies,  and  constants  describing  vibra¬ 
tional  and  rotational  motion  (including  anharmonic  and  rotation-vibration 
effects).  It  is  also  possible  to  solve  the  Schrodinger  equation  for  the  motion 
of  the  nuclei  subject  to  the  potential  curves,  to  obtain  vibrational  wave- 
functions  for  use  in  transition  probability  calculations.  The  electronic 
wavefunctions  themselves  can  be  used  to  estimate  dipole  moments  of  individual 
electronic  states,  transition  moments  between  different  electronic  states, 
and  other  properties.  While  all  of  the  calculations  described  in  this  para¬ 
graph  have  been  carried  out  on  some  systems,  the  unavailability  of  good 
electronic  wavefunctions  and  potential  curves  has  limited  actual  studies  of 
most  of  these  properties  to  a  very  small  number  of  molecules. 

It  is  not  always  possible  or  practical  to  perform  an  ab  initio  or  mathe¬ 
matically  rigorous  calculation  of  the  electronic  structure  of  a  diatomic  mole¬ 
cule  containing  a  large  number  of  electrons  and  it  is  sometimes  useful  to 
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resort  to  approximations  based  on  certain  physical  insights.  One  such  approxi¬ 
mation  is  the  density  functional  method  (ref.  32)  as  applied,  to  molecules. 

In  one  implementation  of  this  approach,  namely  the  Xa  model,  the  one-electron 
Schrodinger  equation  is  set  up  within  an  approximation  to  the  true  potential 
which  is  spherically  symmetric  within  spheres  surrounding  the  several  nuclei, 
constant  in  the  region  between  adjacent  spheres  and  spherically  symmetric 
outside  an  outer  sphere  which  encloses  the  entire  collection  of  nuclei.  The 
resulting  equations  are  solved  by  a  multiple-scattering  method,  equivalent  to 
the  KKR  method  (ref.  33)  often  used  for  crystals.  From  the  resulting  charge 
densities,  one  can  compute  a  total  energy,  using  a  statistical  approximation 
for  the  exchange  correlations.  This  approximation  is  defined  apart  from  a 
multiplicative  factor,  a,  hence  the  name  Xa  method. 

The  Xa  method  has  been  very  successful  in  predicting  excitation  energies 
in  molecules  and  molecular  ions  but  has  met  with  only  limited  success  in  cal¬ 
culation  of  total  energies.  The  basic  deficiency  appeared  to  lie  in  the 
muffin-tin  approximation  to  the  true  potential  rather  than  in  the  statistical 
treatment  of  the  correlation  energy.  However,  recent  studies  at  this  Center 
have  shown  that  when  certain  variational  constraints  are  applied  to  the  X 

o 

wave funct ions ,  significant  improvements  in  the  calculated  total  energies  can 
be  obtained.  At  this  time  the  method  appears  to  have  utility  for  constructing 
regions  of  potential  energy  surfaces. 


16 


SECTION  III 


METHOD  OF  APPROACH 

I  .  QUANTUM  MECHANICAL  CALCULATIONS 

Central  to  the  proposed  studies  are  the  actual  quantum-mechanical 
calculations  which  must  be  carried  out  for  the  molecular  species.  For  added 
clarity,  various  aspects  of  these  calculations  are  discussed  in  individual 

subsections . 

a.  Levels  of  Approximation 

Much  evidence  on  diatomic  and  polyatomic  systems  indicates  the  near 
adequacy  of  a  minimum  Slater-type-orbital  (STO)  basis  for  constructing 
molecular  wavefunct ions  (refs.  34  and  35).  This  means  inner-shell  and 
valence-shell  STO's  of  quantum  numbers  appropriate  to  the  atoms  (Is,  2s,  2p 
for  C,  N,  0;  Is  for  H) .  The  main  deficiency  of  the  minimum  basis  set  is  in 
its  inability  to  describe  polarization  of  »  orbitals  in  atoms  adjacent  to  H 
atoms,  and  successful  calculations  usually  result  if  one  (or  a  set)  of  P 
orbitals  is  supplied  for  each  H  atom.  Values  of  the  screening  parameters  c 
for  each  orbital  can  either  be  set  from  atomic  studies  or  optimized  in  the 
molecule;  the  later  approach  is  indicated  for  studies  of  maximum  precision. 
When  high  chemical  accuracy  is  required,  as  for  the  detailed  studies  of  the 
ground  state  of  a  system,  a  more  extended  basis  should  be  used.  Double-Zeta 
plus  polarization  functions  or  optimized  MO’s  usually  are  required. 

The  chosen  basis  sets  give  good  results  only  when  used  in  a  maximally 
flexible  manner.  This  implies  the  construction  of  Cl  wavefunct ions  with  all 
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kinds  of  possible  orbital  occupancies,  so  that  the  correlation  of  electrons 
into  overall  3tates  can  adjust  to  an  optimum  form  at  each  geometrical  confor¬ 
mation  and  for  each  state.  Except  when  well-defined  pairings  exist  for  as 
many  electrons  as  possible,  a  single-configuration  study  (even  of  Hartree-Fock 
quality)  will  be  inadequate, 
b.  Spin  and  Symmetry 

Proper  electronic  states  for  systems  composed  of  light  atoms  should 
possess  definite  eigenvalues  of  the  spin  operator  S2  as  well  as  an  appro¬ 
priate  geometrical  symmetry.  The  geometrical  symmetry  can  be  controlled  by 
the  assignment  of  orbitals  to  each  configuration,  but  the  spin  state  must  be 
obtained  by  a  constructive  or  projective  technique.  Formulas  have  been 
developed  (ref.  36)  for  projective  construction  of  spin  states  from  orthogonal 
orbitals,  and  programs  implementing  these  formulas  have  for  several  years 
been  in  routine  use  at  UTRC. 

One  of  the  least  widely  appreciated  aspects  of  the  spin-projection 
problem  is  that  the  same  set  of  occupied  spatial  orbitals  can  sometimes  be 
coupled  to  give  more  than  one  overall  state  of  given  S  quantum  number.  It  is 
necessary  to  include  in  calculations  all  such  spin  couplings,  as  the  optimum 
coupling  will  continuously  change  with  changes  in  the  molecular  conformation. 
This  is  especially  important  in  describing  degenerate  or  near-degenerate 
excited  electronic  states. 

c.  Method  of  Ab  Initio  Calculation 

A  spin-free  nonrelativist ic  Hamiltonian  will  be  employed  in  the  Born- 

Oppenheimer  approximation.  For  a  molecule  containing  n  electrons  and  m 

18 


nuclei,  the  approximation  leads  to  an  electronic  Hamiltonian  depending 
parametrically  on  the  nuclear  positions: 


n  m 
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1  =  1  H 
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where  and  are  the  charge  and  position  of  nucleus  p,  r^  is  the  posi- 
tion  of  electron  i,  and  7^  is  the  Laplacian  operator  for  electron  i.  The 
Hamiltonian  H  is  in  atomic  units  (energy  in  Hartrees,  length  in  Bohrs). 

Electronic  wavefunct ions  i(i,  depending  parametrically  on  the  nuclear 
positions  are  made  to  be  optimum  approximations  to  solutions  of  the 
Schrodinger  equation 

J{ R)*(R)=  E  (R)*(R)  (2) 

by  involving  the  variational  principle 

,  0  /*'(  R)JftR)*(R)etT  (3) 

8w(R)=S  - - - 

The  integrations  in  equation  (3)  are  over  all  electronic  coordiates  and  the 
stationary  values  of  W  are  approximations  to  the  energies  of  states  described 
by  the  corresponding  iji.  States  of  different  symmetries  are  studied  by 
restricting  to  the  appropriate  form,  and  excited  states  of  any  symmetry  can  be 
handled  by  simultaneously  determining  the  i|i  for  the  excited  state  and  all 
lower-lying  states  of  the  same  symmetry. 
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The  specific  form  for  i|>  may  be  written 


*(R)  =  I  C„^^(R) 

H- 


(4) 


where  each  is  referred  to  as  a  configuration,  and  has  the  general 
structure 


^{r)  =  /?l9s  n  (Li,R)eu 


(5) 


where  each  i|)^£  is  a  spatial  orbital,  /J  is  the  antisymmetrizing  operator, 

&s  the  spi-n~Pr°ject ion  operator  for  spin  quantum  number  S,  and  0M  is  a 
product  of  a  and  8  one-electron  spin  functions  of  magnetic  quantum  number  M. 
No  requirement  is  imposed  as  to  the  double  occupancy  of  the  spatial  orbital, 
so  equations  (4)  and  (5)  can  describe  a  completely  general  wavefunction. 

In  Hartree-Fock  calculations  i i  is  restricted  to  a  single  ^  which  is 
assumed  to  consist  as  nearly  as  possible  of  doubly-occupied  orbitals.  The 
orbitals  are  then  selected  to  be  the  linear  combinations  of  basis 

orbitals  best  satisfying  equation  (3).  Writing 


V'  =  2  <Vi  X. 

r  V 


(6) 


t^e  avi  are  determined  by  solving  the  matrix  Hartree-Fock  equations 

l  F\v°v\  =  *\l  Sx^i  X)  (7) 

where  C£  is  the  orbital  energy  of 

The  Fock  operator  F^y  has  been  thoroughly  discussed  in  the  literature 
(ref.  37)  and  depends  upon  one-  and  two-electron  molecular  integrals  and  upon 
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the  avi‘  T*1^'s  malces  equation  (7)  nonlinear  and  it  is  therefore  solved 
iteratively.  UTRC  has  developed  programs  for  solving  equation  (7)  for  both 
closed-  and  open-shell  systems,  using  basis  sets  consisting  of  Slater-type 
atomic  orbitals.  Examples  of  their  use  are  in  the  literature  (ref.  38). 

In  configuration  interaction  calculations,  the  sumation  in  equation  (4) 
has  more  than  one  term,  and  the  c^  are  determined  by  imposing  equation  (3), 
to  obtain  the  secular  equation 


Z  (H^„- wsp)cv  =  0  (each/j.) 


(8) 


where 

H/W*;  [R)J  (R)v„(R)arT 

sMv=  / V  {R)*U  (R )dr 


(9) 


Equation  (7)  is  solved  by  matrix  diagonalization  using  either  a  modified 
Givens  method  (ref.  30)  or  a  method  due  to  Shavitt  (ref.  31). 

The  matrix  elements  and  may  be  reduced  by  appropriate  operator 
algebra  to  the  forms 
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where  P  is  a  permutation  and  its  parity.  The  sum  is  over  all  permutations. 

p 

<0MI  ®  S  I ^M>  a  "Sanibel  coefficient"  and  the  remaining  factors  are 
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spatial  integrals  which  can  be  factored  into  one-  and  two-electron  integrals. 
If  the  i|i^  are  orthonomal,  equations  (10)  and  (11)  become  more  tractable 
and  the  H  and  may  be  evaluated  by  explicit  methods  given  in  the 
literature  (ref.  36).  Computer  programs  have  been  developed  for  carrying 
out  this  procedure,  and  they  have  been  used  for  problems  containing  up  to  40 
total  electrons,  10  unpaired  electrons,  and  several  thousand  configurations. 

The  Cl  studies  described  above  can  be  carried  out  for  any  orthonormal 
set  of  ^  for  which  the  molecular  integrals  can  be  calculated.  Programs 
developed  by  URTC  make  specific  provision  for  the  choice  of  the  i|i  ^  as 
Slater-type  atomic  orbitals,  as  symmetry  molecular  orbitals,  as  Hartree-Fock 
orbitals,  or  as  more  arbitrary  combinations  of  atomic  orbitals. 

d.  Molecular  Integrals 

The  one-  and  two-electron  integrals  needed  for  the  above  described  method 
of  calculation  are  evaluated  for  STO's  by  methods  developed  by  the  present 
investigators  (ref.  39).  All  needed  computer  programs  have  been  developed 
and  fully  tested  at  UTRC. 

e.  Configuration  Selection 

Using  a  minimum  basis  plus  polarization  set  of  one-electron  functions,  a 
typical  system  can  have  of  the  order  of  10^  configurations  in  a  full  Cl 
(that  resulting  from  all  possible  orbital  occupancies).  It  is  therefore 
essential  to  identify  and  use  the  configurations  describing  the  significant 
part  of  the  wavefunct ion .  There  are  several  ways  to  accomplish  this  objec¬ 
tive.  First,  one  may  screen  atomic-orbital  occupancies  to  eliminate  those 
with  excessive  formal  charge.  Alternatively,  in  a  molecular-orbital  framework 
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one  may  eliminate  configurations  with  excessive  numbers  of  antibonding 
orbitals.  A  third  possibility  is  to  carry  out  an  initial  screening  of 
configurations,  rejecting  those  whose  diagonal  energies  and  interaction 
matrix  elements  do  not  satisfy  significance  criteria.  Programs  to  sort 
configurations  on  all  the  above  criteria  are  available  at  UTRC. 

Other,  potentially  more  elegant  methods  of  configuration  choice  involve 
formal  approaches  based  on  natural-orbital  (ref.  40)  or  multiconfiguration 
SCF  (ref.  41)  concepts.  To  implement  the  natural-orbital  approach,  an 
initial  limited-CI  wave-function  is  transformed  to  natural-orbital  form,  and 
the  resulting  natural  orbitals  are  used  to  fcrm  a  new  Cl.  The  hoped-for 
result  is  a  concentration  of  the  bulk  of  the  Cl  wavefunction  into  a  smaller 
number  of  significant  terms.  The  multiconfiguration  SCF  approach  is  more 
cumbersome,  but  in  principle  more  effective.  It  yields  the  optimum  orbital 
choice  for  a  pre-selected  set  of  configurations.  This  approach  works  well 
when  a  small  number  of  dominant  configurations  can  be  readily  identified. 

f.  Density  Functional  Approach  -  X^  Model 

The  X  model  (ref.  42)  for  the  electronic  structure  of  atoms,  molecules, 
a 

clusters  and  solids  is  a  local  potential  model  obtained  by  making  a  simple 
approximation  to  the  exchange-correlation  energy.  If  we  assume  a  non- 
relativistic  Hamiltonian  with  only  electrostatic  interactions,  it  can  be 
shown  that  the  total  energy  E  of  a  system  can  be  written  exactly  (ref.  43) 
(in  atomic  units)  as 
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This  expression  is  exact  provided  the  u^  are  natural  orbitals  and  n^  are 
their  occupation  numbers  (i.e.,  eigenfunctions  and  eigenvalues  of  the  first 
order  density  matrix).  The  first  term  in  equation  (12)  represents  the 
kinetic  and  electron-nuclear  energies.  The  second  term  is  the  nuclear 
repulsion  energy.  The  sums  (p,v)  are  over  all  the  nuclear  charges  in  the 
system.  The  third  term  is  the  electron-electron  repulsion  term,  which 
represents  the  classical  electrostatic  energy  of  the  charge  density  p  inter¬ 
acting  with  itself,  where 


pU)  =X  n.u*(l)  u.(l) 


(13) 


The  last  term  Exc  represents  the  exchange  correlation  energy  and  can  be 
expressed  formally  as 

E*c  =  T / P^'1  c*ri  I  P*c  ^,2^.  dr^  ,  (14) 

12 

where  pxc  (1,2)  represents  the  exchange-correlation  hole  around  an  electron 
at  position  1.  In  the  exact  expression,  pxc  is  dependent  on  the  second-order 
density  matrix.  In  the  Hartree-Fock  approximation  Exc  is  the  exchange 
energy,  pxc  represents  the  Fermi  hole  due  to  the  exclusion  principle  and 
depends  only  on  the  first-order  density  matrix.  In  the  method,  we  make  a 
simpler  assumption  about  Pxc*  If  we  assume  that  the  exchange-correlation 
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hole  is  centered  on  the  electron  and  is  spherically  symmetric,  it  can  be 
shown  that  the  exchange-correlation  potential 


U«C -f 


V1'21  df 

r  2 


12 


(15) 


is  inversely  proportional  to  the  range  of  the  hole,  rg,  where  rg  is 
defined  by 


47T  , 

3  r  s  /><')  =  1 

Therefore,  in  the  Xq  model,  the  potential  Uxc 
We  define  a  scaling  parameter  a  such  that 


is 


proportional  to 


(16) 


uX(y^  =  ~  ( 3/> (l)/87r)l/3 

The  expression  in  equation  (17)  is  defined  so  that  a  =  2/3  for  the  case  of  a 
free  electron  gas  in  the  Hartree-Fock  model  (ref.  44)  and  a  =  1  for  the 
potential  origially  suggested  by  Slater  (ref.  45).  A  convenient  way  to 
choose  this  parameter  fcr  molecular  and  solid  state  applications  is  to 
optimize  the  solutions  to  the  equations  in  the  atomic  limit.  Schwarz 
(ref.  46)  has  done  this  for  atoms  from  Z  =  1  to  Z  =  41  and  found  values 
between  2/3  and  1. 

In  the  "spin  polarized"  version  of  the  XQ  theory,  it  is  assumed  (as  in 
the  spin-unrestricted  Hartree-Fock  model)  that  electrons  interact  only  with  a 
potential  determined  by  the  charge  density  of  the  same  spin.  In  this  case  the 
contribution  to  the  total  energy  is  summed  over  the  two  spins,  s  =  +  1/2. 
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Exc:  2  ^  //°s(l)  Uxa,s(l,d^ 


(18) 


where  the  potential  is  sp in-dependent 


*df«s 


(l)=  (3^)  |  J^-jr)1^3 


(19) 


and  pg  is  the  charge  density  corresponding  to  electrons  of  spin  j.  The 
spin-polarized  Xa  model  is  useful  for  describing  atoms  and  molecules  with 

open-shell  configurations  and  crystals  which  are  ferromagnetic  or  anti- 
ferromagnet ic . 

Once  one  has  made  the  Xa  approximation  to  the  total  energy  functional  E 
in  equation  (12),  then  the  rest  of  the  theory  follows  from  the  application  of 
variational  principle.  The  orbitals  ^  are  determined  by  demanding  that  E 
be  stationary  with  respect  to  variations  in  This  leads  to  the  set  of 

one-electron  X^  equations 


+  1 

H- 


J  ri2  2  3  uxa 


ui  =€iui 


(20) 


where  ^  is  the  one-electron  eigenvalue  associated  with  Ui>  Since  p 
(?)  is  defined  in  terms  of  the  orbitals  Mi,  the  equation  (20)  must  be 
solved  iteratively,  until  self-consistency  is  achieved.  Empiricaly,  if  one 
takes  as  an  initial  guess  that  p  is  approximately  a  sum  of  superimposed 
atomic  charge  densities,  then  the  convergence  of  this  procedure  is  fairly 
rapid.  The  factor  of  2/3  multiplying  the  potential  is  a  result  of  the 
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nonlinear  dependence  of  Exc  on  p  This  also  has  as  a  consequence  that  the 


Xa  eigenvalues  e^  do  not  satisfy  Koopraans'  theorem,  i.e.,  they  cannot  be 
interpreted  as  ionization  energies.  However,  it  can  be  shown  that  the 
are  partial  derivatives  of  the  total  expression  of  equation  (12)  with  respect 
to  the  occupation  number, 


€. 


dt_ 

dn\ 


(21) 


If  E  were  a  linear  function  of  n^ ,  then  Koopmans 1  theorem  would  hold. 
However,  because  of  the  dominant  Coulomb  term,  E  is  better  approximated  by  a 
quadratic  function  in  n^ .  This  leads  to  the  "transition  state"  approxima¬ 
tion  which  allows  one  to  equate  the  difference  in  total  energy  between  the 
state  (n^,  nj)  and  (n^-1 ,  nj+1)  to  the  difference  in  the  one-electron 
energies  calculated  in  the  state  (n^-1/2,  nj+1/2).  The  error 

in  this  approximation  is  proportional  to  third-order  derivatives  of  E  with 
respect  to  n^  and  nj ,  which  are  usually  small  (ref.  47).  The  main 
advantage  of  using  the  transition  state  rather  than  directly  comparing  the 
total  energy  values  is  computational  convenience,  especially  if  the  total 
energies  are  large  numbers  and  the  difference  is  small. 

The  relationship  of  equation  (21)  also  implies  the  existence  of  a  "Fermi 
level"  for  the  ground  state.  This  can  be  seen  by  varying  E  with  respect  to 
n^  under  the  condition  that  the  sum  Z’n-  is  a  constant,  i.e., 

■i  ^ 


(22) 


irapli es  =  X , 


where  X  is  a  Lagrangian  multiplier. 


This  implies  that  the 
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total  energy  is  stationary  when  all  the  one-electron  energies  are  equal.  How 
e  occupation  numbers  are  also  subject  to  the  restriction  0  <  n^  £  1 . 
s  to  the  following  conditions  on  the  ground  state  occupation  numbers 

*  <X  )•  n.  =i 

•  i 

>  X  •)■  nj  =0  (23) 

°-  nj  -  1 

In  other  words,  the  ground  state  eigenvalues  obey  Fermi  statistics  with  X 
representing  the  Fermi  energy.  It  should  be  noted  that,  in  contrast  to  the 
Hartree-Fock  theory,  where  all  the  n^^  are  either  0  or  1 ,  the  Xa  model 
predicts,  in  some  cases,  fractional  occupation  numbers  at  the  Fermi  level. 

In  particular,  this  will  occur  in  a  system  (such  as  transition  metal  or 
actinide  atom)  which  has  more  than  one  open  shell. 

The  Xa  model  differs  in  other  significant  ways  from  the  Hartree-Fock 
method.  In  fact,  the  simplification  introduced  in  approximating  the  total 
energy  expression  introduces  several  distinct  advantages  over  Hartree-Fock: 

1.  The  primary  advantage  is  purely  computational.  The  one-electron 
potential  in  equation  (20)  is  orbital-independent  and  local,  i.e.,  it  is  the 
same  for  all  electrons  (except  in  the  spin-polarized  Xa  theory)  and  is  a 
multiplicative  operator.  On  the  other  hand,  the  Hartree-Fock  potential 
is  non-local,  or  equivalently,  there  is  a  different  local  potential  for  each 
orbital.  This  involves  a  great  deal  more  computational  effort,  especially 
for  systems  described  by  a  large  number  of  orbitals.  It  has  been  shown  (ref. 
48)  that  the  Xa  orbitals  for  the  first  and  second  row  atoms  are  about  as 
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accurate  as  a  double-zeta  basis  set,  and  are  probably  better  for  larger  atoms 
which  involve  electrons  with  4,  >  2. 

2.  The  orbital-independent  potential  leads  to  a  better  one-electron 
description  of  electronic  excitations  of  a  system.  Both  the  unoccupied  (n^  =  0) 
and  occupied  (n^  =  1)  eigenfunctions  are  under  the  influence  of  the  same 
potential  resulting  from  the  other  N-l  electrons.  The  Hartree-Fock  virtual 
orbitals  see  a  potential  characteristic  of  the  N  occupied  orbitals,  and 
therefore  are  not  as  suitable  for  describing  the  excited  state.  Actually, 
although  the  ground  state  virtual  eigenvalues  are  usually  a  good  description 

of  the  one-electron  excitations,  the  virtual  spectrum  of  the  transition  state 
potential  where  one-half  an  electron  has  been  removed  from  the  system  give  a 
much  better  first-order  picture  of  these  levels  (ref.  49). 

3.  As  has  been  shown  by  Slater  (ref.  50),  the  Xq  model  rigorously 
satisfied  both  the  virial  and  Hellman-Feynman  theorems,  independent  of  the 
value  of  the  parameter  a-  This  is  convenient  for  calculating  the  force  on  a 
nucleus  directly  in  terms  of  a  three  dimensional  integral,  rather  than  the 
six  dimensional  integrals  in  the  expression  for  the  total  energy  of  equation 
(12). 

g.  Computational  Aspects  of  the  X  Model 

a 

In  application  of  the  X^  model  to  finite  molecular  systems,  there  are 
two  practical  aspects  of  the  calculations  which  must  be  considered.  The 
first  concerns  the  choice  of  the  integration  framework  for  describing  the 
molecular  wavefunct ions  and  the  second  deals  with  the  choice  of  the  exchange 
parameter,  a,  in  different  regions  of  space. 
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In  computations  with  heteronuclear  molecules,  there  are  several  free 
parameters  that  must  be  chosen:  the  ratio  of  sphere  radii  for  the  atomic 
spheres  of  integration  at  a  given  internuclcar  separation,  the  degree  of 
sphere  overlap,  and  the  value  of  the  exchange  parameter  in  the  atomic 
spheres  and  the  intersphere  region. 

It  has  been  found  that  changing  the  ratio  of  the  sphere  radii  for  the  two 
atoms  in  a  heteronuclear  diatomic  molecule  introduces  changes  in  the  total 
energy  that  can  be  large  on  a  chemical  scale,  (v>  1  eV).  A  choice  for  sphere 
radii  based  on  covalent  bonding  radii  does  not  necessarily  provide  a  good 
estimate  for  these  calculations.  The  value  of  the  exchange  parameter,  a,  and 
the  sphere  radii  and/or  sphere  overlap  is  normally  fixed  in  Xa  calculations 
for  crystals  where  the  geometry  is  fixed.  However,  to  develop  a  potential 
curve,  the  molecular  description  needs  to  change  substantially  as  the  inter- 
nuclear  separation  varies  and  the  changing  sphere  radii  include  varying 
fractions  of  the  total  molecular  charge  (ref.  51).  Studies  made  at  UTRC  have 
shown  that  at  any  given  separation  the  total  energy  calculated  from  the  Xa 
model  is  a  minimum  at  the  radii  ratio  where  the  spherically  averaged  poten¬ 
tials  from  the  two  atomic  centers  is  equal  at  the  sphere  radius, 

VrSl)  =  V2^s2)  (24) 

This  relationship  between  the  potential  match  at  the  sphere  boundary  and  the 
minimum  in  the  total  energy  appears  to  hold  exactly  for  "neutral"  atoms  and 
holds  well  for  ionic  molecular  constituents.  In  the  case  of  two  ionic 
species,  the  long  range  tail  of  the  potential  must  go  like  +2/R  from  one  ion 
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and  -2/R  (in  Rydbergs)  for  the  other  ion  and  so  at  large  internuclear  separa¬ 
tions,  the  1/R  character  of  the  potential  does  not  invalidate  the  potential 
match  criterion  for  radii  selection.  This  match  for  the  atomic  potentials  is 
applied  to  the  sel. "-cons is tent  potentials. 

In  molecules  with  significant  charge  sharing  in  the  bonds,  the  radii  of 
the  atomic  spheres  is  frequently  increased  in  Xa  calculations  so  that  an 
overlap  region  appears  in  the  vicinity  of  the  bond  (ref.  52).  Studies  made 
at  UTRC  show  that  the  contribution  to  the  total  molecular  energy  from  the 
exchange  integral  shows  a  minimum  at  the  optimum  sphere  radius  or  sphere 
overlap.  This  provides  a  sensitive  criterion  for  selecting  these  parameters. 

The  values  of  the  exchange  parameters  in  the  spherical  integration  region 
around  each  atomic  center  are  frequently  set  at  the  atomic  values  both  for 
neutral  and  for  ionic  molecular  constituents.  However,  for  light  atoms,  the 
value  of  a  which  best  reproduces  Hartree-Fock  results  varies  substantially 
with  ionicity.  In  argon,  the  following  table  compares,  for  the  neutral  atom 
and  the  positive  ion,  the  HF  energy  and  the  Xa  energy  calculated  for 
several  values  of  a. 

TABLE  1.  COMPARISON  OF  CALCULATED  ATOMIC  ENERGIES  FOR  ARGON 


a  Energy  HF  Energy 


Ar° 

.72177 

-526.8176 

-526.8173 

Ar+1/2 

.72177 

-526.5857 

_ 

.72213 

-526.6007 

— 

Ar+1 

.72177 

-526.2447 

_ 

.72213 

-526.2596 

— 

.72249 

-526.2745 

-526.2743 
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The  optimum  value  of  a  changes  even  more  rapidly  in  the  fluorine  atom,  going 
from  0.73732  for  F°  to  0.72991  for  F  * .  Since  the  total  energy  depends 
linearly  on  a,  this  parameter  must  be  chosen  carefully. 

The  intersphere  exchange  coefficient,  is  chosen  to  be  a  weighted  average 
of  the  atomic  exchange  parameters  from  the  two  constituents.  At  small 
internuclear  separations,  the  optimum  radius  for  an  atomic  sphere  frequently 
places  significant  amounts  of  charge  outside  that  atomic  sphere  -  charge  that 
is  still  strongly  associated  with  its  original  center  rather  than  being 
transferred  to  the  other  center  or  associated  with  the  molecular  binding 
region.  To  best  account  for  these  cases  the  weighting  coefficients  are 

chosen  to  reflect  the  origin  of  the  charge  in  the  intersphere  (or  outersphere 
region) , 


or. 


intersphere 


\  \  -  ^  -  V> 

\  -  <0  *  \  -  «2°) 


(25) 


where  (Q  -  Qj0)  is  the  charge  lost  from  sphere  i  relative  to  its  atomic 
value  (or  ionic  value)  Qi°  and  ag^  is  the  atomic  exchange  parameter  for  sphere 
his  value  for  a£ntersphere  calcu^ateh  dynamically  -  it  is  updated  after 
each  iteration  in  the  self-consistent  calculation. 

While  for  heavy  atoms,  these  changes  in  the  exchange  parameter  would  be 
small,  the  a's  for  small  atoms  vary  rapidly  with  z  (and  with  ionicity).  The 
correct  choice  of  the  exchange  parameters  influences  not  only  the  total 
energy  calculated  for  the  molecule  but  also  in  some  cases  affects  the  distri¬ 
bution  of  charge  between  the  atomic  spheres  and  the  intersphere  region. 
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SECTION  IV 


DISCUSSION  OF  RESULTS 


The  theoretical  research  conducted  under  this  program  was  concerned 
with  a  study  of  the  energetics  of  the  noble  gas  halides  and  noble  gas 
dimer  and  trimer  molecular  ions.  The  emphasis  of  the  study  of  the  noble 
gas  halides  was  on  the  definition  of  potential  energy  curves,  prediction 
of  the  radiative  lifetimes  of  electronically  excited  states,  and  studies  of 
photon  energy  loss  paths  such  as  photoabsorption  and  photoionization  of  the 
upper  excited  electronic  states.  The  major  effort  of  the  study  of  the 
noble  gas  dimer  and  trimer  ions  was  the  definition  of  the  potential  energy 
curves  and  predicition  of  the  photoabsorption  cross-section  for  these  species. 

For  all  of  the  noble  gas  halide  systems  there  is  an  outer  crossing  of 
the  upper  ionic  states  arising  from  R+  +  X~  (R  =  noble  gas,  X  =  halide)  such 
that  these  upper  molecular  states,  at  their  equilibrium  intemuclear  separa¬ 
tions,  actually  lie  lower  than  any  of  the  Rydberg  molecular  states. 

The  possible  electronic  transitions  which  correspond  to  laser  activity 
in  these  systems  therefore  reduce  to  the  following  set: 

a2i+  -*  xh+ 

B2n  +  x2e+ 

2  +  2 
a  e  -*■  hi 

2  2 

B^n  -»■  m 


2  2  ? 

Further,  the  B  II  state  is  split  into  II  ^  and  ni/2  states  owing  to 
spin-orbit  coupling  effects  in  the  noble  gas  positive  ion.  This  splitting 
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is  small  (782  cm  for  Ne  but  is  appreciable  for  Xe  (10537  cm-^) . 

As  a  result  of  these  several  possibilities,  the  definite  assignment 

of  the  laser  transitions  becomes  difficult.  Brau  and  Ewing  (refs.  3  and 

53)  have  assigned  the  strong  laser  transitions  in  XeBr,  XeCl,  XeF  and  KrF 
2  +  2  + 

to  A  ^1/2  X  S i/ 2  in  accord  with  a  simple  ionic  model  for  the  potential 

curves.  Ab  initio  calculations  by  Dunning  and  Hay  (Ref.  7)  for  KrF  also 
2  +  2  "b 

suggest  the  A  ^1/2  ^  ^1/2  trans^t^on  corresponds  to  the  observed  emission 

o  2 

in  KrF  at  2480A.  Their  potential  curves,  however,  also  indicate  a  B 

^3/2  ^  transition  at  nearly  the  same  calculated  wavelength.  Since  this  is 

also  a  strong  parallel  transition,  a  unique  assignment  does  not  appear  to 

be  possible  on  the  basis  of  Dunning  and  Hay  calculations.  Further,  their 

2 

calculations  indicate  a  near  degeneracy  of  the  upper  ionic  A  ^1/2  and 
2 

B  ILj/2  states  which  is  difficult  to  understand  in  terms  of  simple  charge 

density  models  for  these  excited  states. 

Clugston  and  Gordon  (ref.  6)  suggest  that  the  dominant  laser  transition 
2  2  + 

is  B  II3/2  X  based  on  free-electron  gas  model  calculations.  Recent 

studies  carried  out  at  this  Center,  however,  question  the  reliability  of 
non-SCF  model  calculations  for  these  systems.  In  addition,  this  assignment 
of  a  perpendicular  transition  to  the  laser  emission  goes  against  the  experi¬ 
mental  observations  that  such  transitions  should  be  weak  compared  to  allowed 
transitions . 

In  order  to  shed  some  light  on  the  above  questions,  both  ab  in it io  and 
density  functional  X^  calculations  were  carried  out  for  the  lower  and 
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excited  levels  of  ArF  and  A^F.  The  ab  initio  calculations  were  carried  out 
in  a  manner  similar  to  that  described  previously  (ref.  57).  All  nonlinear 
parameters  of  the  molecular  orbital  basis  functions  were  individually  opti¬ 
mized  for  each  separate  electronic  state.  The  results  of  these  studies  have 
been  analyzed  and  published  as  technical  papers.  For  completeness  they  are 
presented  herein  as  Appendices  A  and  B. 

The  theoretical  principles  underlying  the  electronic  structure  of  the 
homonuclear  noble  gas  molecular  ions  have  been  discussed  by  Mulliken  (ref. 

54)  using  Xe^  as  a  prototype  and,  more  recently,  by  Gilmore,  Barr,  and  Dee 

(ref.  55)  for  the  case  of  Kr^.  For  these  dimer  ions,  four  electronic  states 
o+  2  2  2  +  1 

(A  £  ,  B  II  ,  C  H  ,  and  D  £  )  arise  from  a  ground  state  atom  (  Sn)  and 

2  2  + 

a  ground  state  ion  (  P^/2  3/2^'  The  lowest  state,  usually  labeled  A  Zu> 

+  +  +  “f’ 

is  bound  for  the  four  systems  Ne2,  Ar^,  Kr^,  and  Xe£  and  has  been  partially 

characterized  experimentally  through  scattering  experiments,  although  there 

has  been  much  controversy  concerning  the  dissociation  energies  of  these  ions. 

(ref.  54).  For  the  heavier  ions,  A-S  coupling  is  no  longer  strictly  valid 
2  2 

and  the  B  Tig  and  C  states  split  into  J-J  coupled  components  with  a 

resultant  mixing  of  £  and  n  symmetries.  This  splitting  is  small  for  Ne^ 

but  six  distinct  J-J  coupled  molecular  ion  states  arise  for  Xe^. 

With  the  exception  of  the  recent  calculations  by  Wadt,(ref.  56)  no 

attempts  at  systematic  studies  of  these  ions  have  been  reported.  Accurate 

ab  initio  methods  would  appear  to  be  practical  only  for  the  lighter  ions 

+ 

with  diminished  accuracy  for  heavier  systems  such  as  Xe2 .  Wadt  reports 
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such  a  trend  ln  hls  “lculatl°-  —b  -it-  that  hls  caussIan 

basis  is  less  flexible  or  less  complete  in  the  case  of  Xe^  than  .for  the 

corresponding  study  of  Ar+.  Recent  density  functional  calculations  of  noble 

*"  halide  molecules  (ref.  57,  ref.  58)  have  yielded  reliable  predictions 

of  geometry  and  chemical  binding,  provided  certain  optimisation  criteria 

ate  followed.  Such  calculations,  uhich  should  be  of  more  uniform  quality  than 

ab  tgitio  studies,  since  numerical  molecular  orbitals  not  limited  by  basis 

set  constraints  are  employed  (ref  fa\  .. 

pxoyed  (ref.  59),  constituted  the  principal  effort  in 

the  present  study.  A  limited  number  of  ab  initio  ,  7  7  • 

or  ab  initio  calculations  were  carried 

out  for  Ar*  for  comparison  purposes. 

Tbs  results  of  „ur  studies  of  the  noble  gas  molecular  ions  have  been  pre- 

publication  as  technical  papers.  For  completeness  our  results  are 
presented  herein  as  Appendices  C,  D  and  J. 

A  brief  summary  of  the  most  important  results  that  „ere  obtained  for  each 

of  the  molecular  species  that  were  examine  • 

wore  examined  is  given  below. 

ArF 

Our  predicted  emission  spectra  for  ArF  indicates  that  only  the 
2  +  2  + 

B  Zl/2  *  X  Zi/ 2  transition  exhibits  a  large  transition  moment.  Our  cal¬ 
culated  radiative  lifetime  is  5  ns  for  this  transition  and  our  calculated 
peak  emission  intensity  is  at  a  wavelength  of  196  nm,  as  compared  with 
the  experimental  laser  frequency  of  193  nm. 

Ar2F 

Calculations  for  ArpF  indicate  that  the  bound  upper  ionic  state  has 
2B2  symmetry  with  Ar-Ar  and  Ar-F  bond  lengths  similar  to  those  in  the  corres¬ 
ponding  diatomic  species.  The  terminating  state  of  the  emissl0n  from  this 
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trimer  species  has  2B2  symmetry  while  the  ground  state  has  \  symmetry. 
The  B2  -►  B2  transition  has  a  radiative  lifetime  of  200  ns  with  peak 
emission  at  280  nm.  The  observed  emission  for  Ar2F  is  very  broad  (50  nm) 
and  peaks  at  about  290  nm. 


Ne, 


We  rind  the  A  2^1/2u  ground  state  of  Ne2  to  be  bound  by  1.31  eV  with 
an  equilibrium  separation  of  1.69  %.  Our  calculated  fundamental  vibrational 
frequency  is  597.0  cm-1.  Our  predicted  maximum  photoabsorption  cross- 
section  at  300  K  is  0.20  A  at  a  wavelength  of  285  nm. 


Ar, 


The  A  2Z1/2u  ground  state  of  Ar*  is  bound  by  1.30  eV  at  an  equilibrium 

separation  of  2.43  A.  The  fundamental  vibrational  frequency  is  calculated 
-1 

be  297.9  cm  .  We  predict  a  maximum  photoabsorption  cross-section  at 
o  o 

300  K  of  0.48  A  at  a  wavelength  of  300  nm. 


Kr, 


The  dissociation  energy  of  the  aV^  ground  state  of  Kr*  is  1.18 
eV  when  spin-orbit  effects  are  properly  taken  into  account.  Our  calculated 
equilibrium  separation  is  2.75  X  and  the  fundamental  frequency  is  found 

to  be  176.7  cm  .  We  find  a  maximum  photabsorption  cross-section  at  300°K 
of  0.57  1  at  a  wavelength  of  320  nm. 

Xet 


The  A  ^1/2u  ground  state  of  this  ion  is  bound  by  1.06  eV  at  an 
equilibrium  separation  of  3.18  A.  The  calculated  fundamental  frequency 
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-1  0 
is  117.2  cm  .  Our  predicted  photoabsorption  spectrum  at  300  K  peaks  at 

o 

340  nm  at  a  value  of  0.68  A.  We  find  a  strong  temperature  dependence  of 
the  photoabsorption  spectrum  owing  to  the  smallness  of  the  vibrational 
spacing  of  the  ground  state. 


An  analysis  of  the  electronic  structure  of  Ar^  was  carried  out  in 
°3h»  ^2v  an<*  D°°h  symmetr:i-es  •  The  lowest  state  of  Ar  corresponds  to  the  degen¬ 
erate  2E'  symmetry  with  an  indicated  dissociation  energy  of  0.17  eV  relative  to 

Ar2+  (A  ZE^)  +  AR  (  SQ) .  This  trimer  ion  exhibits  a  small  Jahn-Tellcr  dis¬ 
tortion  from  symmetry.  We  find  a  wide,  shallow  well  for  the  ground 
state  with  an  indicated  Rg  of  2.9  X.  The  asymmetric  stretch  frequency  is 
calculated  to  be  390  cm 

A  strong  and  very  broad  photoabsorption  band  (>*500-670  nm)  is  pre- 
2  i  2 

dieted  for  tht  *"E  -*  A-^'  transition  of  this  trimer  ien.  This  absorption 

°2 

exhibits  a  peak  cross-section  of  0.3  A  at  a  wavelength  of  530  nm.  This 
suggests  that  such  species  may  be  very  important  in  the  analysis  of  loss 
mechanisms  in  visible  excimer  lasers. 
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APPENDIX  A 


THE  ELECTRONIC  STRUCTURE  OF  ArF  AND  Ar2F* 


L.  A.  Wright  (DYP) 

Kirtland  Air  Force  Base,  New  Mexico  87117 


Calculations  have  been  performed  on  the  electronic  structure,  potential 
energy  curves  and  radiative  transition  probabilities  of  ArF  and  Ar2F.  Our 
predicted  emission  spectra  for  ArF  indicates  that  only  the  B  ‘E  -*■  X  “E 
transition  exhibits  a  large  transition  moment  and  hence  a  short  5  nsec) 
radiative  lifetime.  Calculations  for  Ar2F  indicate  that  the  bound  upper  ionic 
state  has  ^B2  symmetry  with  Ar-Ar  and  Ar-F  bond  lengths  similar  to  those  in 
the  corresponding  diatomic  species.  The  terminating  state  also. has  2B2  symme¬ 
try  and  this  polyatomic  system  should  exhibit  a  relatively  long  radiative 
lifetime  (9/  200  nsec)  . 


Supported  in  part  by  AFOSR  under  Contract  No.  F44620-73-C-0077 . 
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1.  INTRODUCTION 


The  diatomic  molecules  formed  from  the  noble  gases ,  from  the  partially 
closed  shell  Group  II  B  elements  (Hg,  Cd,  Zn)  or  from  the  noble  gas -halogen 
combinations  constitute  an  interesting  and  important  class  of  molecules  for 
laser  applications.  These  molecules  are  characterized  by  a  repulsive  or  very 
weakly  bound  ground  state  potential  curve  and  by  bound  electronically  excited 
states.  These  excited  diatomic  molecules,  known  as  excimers,  radiate  in  a 
narrow  visible  or  UV  continuum  band  by  making  a  transition  to  a  repulsive  lower 
state  that  dissociates  to  ground  state  neutral  atoms.  Laser  behavior  is  con¬ 
trolled  by  competing  kinetic  processes  for  the  formation  and  destruction  of 
the  excited  molecular  levels .  The  character  of  the  electronic  structure  of 
the  excited  states  in  the  noble  gas  halide  systems  is  strongly  ionic  and  is 
similar  to  that  of  the  ground  state  alkali  halide  salts.  The  other  molecules 
exhibit  mainly  Rydberg  character  in  their  excited  electronic  states. 

Emission  spectra  for  several  noble  gas  halides  have  been  obtained  in  low- 

pressure  discharge  flow  experiments  and  were  first  reported  by  Velazco  and 
1  p 

Setser-1-  and  by  Golde  and  Thrush.  More  recently  emission  spectra  and  laser 
action  in  several  noble  gas  halide  systems  have  been  reported. The  spectra 
are  consistent  with  a  model  whereby  the  lasing  transition  is  between  a  strongly 
bound  upper  electronic  state  and  a  nearly  flat  repulsive  ground  state.  A  group 
theoretical  analysis  of  the  possible  electronic  configurations  for  these  states 
does  not  lead  to  a  unique  assignment  since  both  the  upper  and  lower  levels 
could  be  of  £  or  U  symmetry. 

Theoretical  studies  of  the  noble  gas  halide  molecules  have  consisted  of 
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analogies  with  alkali  halide  systems, 3  predictions  based  on  an  electron-gas 
model1  and  ab  initio  calculations  for  the  KrF  system. lfa  In  some  cases,  the 
theoretical  methods  utilized  to  date  have  been  unable  to  positively  define  the 
electronic  states  involved  in  the  noble  gas  halide  laser  systems.  Ab  initio 
methods  would  appear  to  be  practical  only  for  the  lighter  systems  such  as  ArF 
and  KrF  owing  to  greatly  increased  computational  complexity  for  larger  numbers 
of  electrons.  However,  detailed  ab  initio  Cl  calculations  on  systems  as  com¬ 
plex  as  XeCl  are  now  in  progress. ^7  The  electron-gas  model  calculations 
suffer  from  the  lack  of  self-consistency  in  the  molecular  potential.  Recent 
density  functional  calculations  indicate  that  self-consistency  in  the  potential 
is  crucial  for  determining  optimum  electronic  wavefunctions . ^  Such  calcula¬ 
tions,  coupled  with  a  limited  number  of  ab  initio  calculations  carried  out  for 
comparison  purposes,  constitute  the  principal  effort  in  the  present  study  of 
excimer  systems. 

2.  ELECTRONIC  STRUCTURE 


For  all  of  the  noble  gas  halide  systems  there  is  an  outer  crossing  of  the 
upper  ionic  states  arising  from  R+  +  X"  (or  R+)  (R  =  noble  gas,  X  =  halide) 
such  that  these  upper  molecular  states,  at  their  equilibrium  internuclear 
separations,  actually  lie  lower  than  any  of  the  Rydberg  molecular  states. 

The  possible  diatomic  electronic  transitions  which  correspond  to  laser 
activity  in  these  systems  therefore  reduce  to  the  following  set: 

B  2e+  X  2E+  B  V  -  A  \ 


C  2!T  -  X  2I+ 


c  2n  -  a  2n 


(1) 


Further,  the  state  is  split  into  antl  ~^\/2  sta^es  owing  to  spin-orbit 
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coupling  effects  in  the  noble  gas  positive  ion.  This  splitting  is  small  (702) 
cm-1)  for  Ne  but  is  appreciable  for  Xe  (10537  cm-1). 


As  a  result  of  these  several  possibilities,  a  definite  assignment  of  the 

laser  transitions  becomes  difficult.  Brau  and  Ewing^’^  have  assigned  the 

strong  laser  transitions  in  XeBr,  XeCl,  XeF  and  KrF  to  B  2E+  ,  -  X  2V+  ,  in 

1/2  h  1/2 

accord  with  a  simple  ionic  model  for  the  potential  curves.  Ab  initio  calcula¬ 
tions  by  Dunning  and  Hay16  for  KrF  also  suggest  the  B  2^+1j2  ~  X  2Z+±j0  transi¬ 
tion  corresponds  to  the  observed  emission  in  KrF  at  2480  A.  Their  potential 
curves  also  indicate  a  C  ^3/2  -*  X  £  l/2  transition  &t  nearly  the  same  calcu¬ 
lated  wavelength.  Further,  their  calculations  indicate  a  near  degeneracy  of 
the  upper  ionic  B  ant^  ^^3/2  s^a^es  which  is  difficult  to  understand  in 

terms  of  simple  charge  density  models  for  these  excited  states. 


Clugston  and  Gordon1^  suggest  that  the  dominant  laser  transition  is 
C  ”3/2  “*  ^  L  l/2 ’  "base(i  on  free  electron  gas  model  calculations.  Recent 
studies,  however,  question  the  reliability  of  non-SCF  model  calculations  for 
these  systems.  In  addition,  this  assignment  of  a  perpendicular  transition  to 
the  laser  emission  goes  against  the  experimental  observations  that  such  transi¬ 
tions  should  be  weak  compared  to  allowed  parallel  transitions. 


3-  CALCULATED  RESULTS  AND  DISCUSSION 


In  order  to  shed  some  light  on  the  above  questions,  both  ab  initio  and 
density  functional  calculations  were  carried  out  for  the  lower  and  excited 
levels  of  ArF  and  A^F.  The  ab  initio  calculations  were  carried  out  in  a 
manner  similar  to  that  described  previously. ^-9  All  nonlinear  parameters  of  the 
molecular  orbital  basis  functions  were  individually  optimized  for  each  separate 
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electronic  state#  The  resultant  ah  initio  results  are  shov*n  in  Fig#  1 
where  we  see  that  the  B  22+1/2  and  C  ^3/2  systems  are  nondegenerate  in  ArF. 
The  upper  C  state  is  less  strongly  polarized  and  exhibits  a  somewhat 

flatter  potential  curve  giving  rise  to  an  inner  crossing  of  A  £  ^  and 

c  a-bou^  A* 

These  calculations  indicate  the  following  long-range  behavior  for  the 


excited  state  potential  curves 

2^+  . 


B  2£  :  V(R)  (a.u. )  =  -  - 


1  1.16 


C  ^  :  V(R)  (a.u.)  = 


R  R- 


-  i  +  0^58 
R  R3 


(2) 


The  R“3  term  arises  from  the  quadrupole  contribution  to  the  interaction  poten¬ 
tial  at  long  range.  This  term  vanishes  for  the  alkali  halide  diatomic  mole¬ 
cules  M+X“  since  both  positive  and  negative  ions  are  in  the  ^  state.  However, 
for  the  noble  gas  halide  systems ,  the  noble  gas  ion  is  in  a  2P  state  with  a 
non-vanishing  quadrupole  moment.  Therefore,  simple  models  of  the  ionic  excited 
states  of  the  noble  gas  halides,  based  on  analogies  with  alkali  halide  systems, 
need  to  be  modified  to  properly  account  for  this  long-range  behavior.  The  cal¬ 
culated  value  for  the  coefficient,  3,  of  the  R  ^  term  has  the  correct  theoretical 
form  such  as  3  o  +  =  -23  P  .  Although  several  system  assignments  are  possible 

n 

within  the  error  of  the  potential  curves,  calculations  based  on  the  potentials 
shown  in  Fig.  1  clearly  indicate  that  only  the  B  22+1y2  -*  X  2Z  ^  transition 
exhibits  a  large  transition  moment  and  hence  a  short  (~  5  nsec)  radiative  life¬ 
time.  Our  calculated  peak  emission  intensity  is  at  a  wavelength  of  196  nm  as 
compared  with  the  experimental  laser  frequency  of  193  nm. 
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Potential  energy  curves  for  ArP  were  also  constructed  using  a  density 
functional  method.  As  described  in  Ref.  18,  our  program  has  been  exten¬ 
sively  developed  beyond  the  standard  techniques  described  in  the  literature. 
Modifications  to  minimize  the  muffin-tin  potential  error  and  to  insure  con¬ 
tinuity  in  the  molecular  potential  function  have  been  incorporated  into  our 
code  which  result  in  more  reliable  total,  energies.  The  results  for  ArP 
(without  including  spin-orbit  mixing)  are  shown  in  Fig.  2  which  indicate 
quantitative  agreement  with  our  ab  initio  results.  Again  we  see  that  the 
C  2n  excited  state  is  nondegenerate  with  the  upper  B  2S+  state  and  that  it 
is  somewhat  flatter  in  curvature.  The  origin  of  the  somewhat  steeper  repul¬ 
sion  exhibited  by  all  the  calculated  potentials  using  this  density  functional 
model  is  now  understood  as  a  manifestation  of  the  exchange  scaling  parameter 
variation  with  the  effective  nuclear  charge. 

Regions  of  the  potential  energy  surfaces  for  Ar2F  have  been  studied  to 
examine  the  possible  stability  of  this  triatomic  molecule  in  an  excited  elec¬ 
tronic  configuration.  It  has  been  suggested20  that  the  long  wavelength 

O  Q 

radiation  seen  at  2900  A  in  e-beam  excited  argon-fluorine  mixtures  (4100  A 
for  krypton-fluorine  mixtures)  may  be  due  to  a  triatomic  molecule  formed  in  an 
excited  ionic  state  from  Ar2  +  F",  similar  to  the  excited  ionic  state  formed 
from  Ar+  +  F“.  Such  a  molecule  would  radiate  at  longer  wavelengths  since 
Arg  is  stable  relative  to  Ar+  +  Ar. 

The  molecular  correlation  diagram  for  Ar2F,  given  in  Table  I  for  CQv 
symmetry,  should  represent  the  most  stable  geometry  for  both  the  excited  ionic 
and  ground  state  surfaces.  The  corresponding  orbital  designations  are  also 
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carried  out  preliminary 


given  for  the  linear  molecule  Ar-F-Ar.  We  have 
calculations  for  the  ground  and  first  excited  ionic  states  in  \  and  2B 
symmetries.  Our  calculated  results  are  shown  in  Fig.  3. 

We  fmd  a  flat  potential  surface  for  Ar2F  in  the  lowest  2 A±  state  for 
%  >  2.8  A  as  shown  in  Fig.  3.  This  results  from  forward  charge  transfer  from 
Ar  to  the  F  atom,  thereby  introducing  some  ionic  character  of  the  type  Ar+X- 
which  will  tend  to  cancel  out  neutral  atom  exchange  and  overlap  effects. 
For  heavier  systems  such  as  Xe2F,  some  small  long-range  binding  may  be  found 
for  this  lowest  2^  configuration.  The  lowest  2B2  state  exhibits  a  somewhat 

steeper  repulsion.  For  Ar2  -  F  separations  shorter  than  2.6  A,  both  the 
and  2B2  lower  3tates  exhih.t  a  rap.dly  increas.ng  repuls.on  („  5  eV/»}>  Thig 

character  leads  to  the  observed  broadband  emission  seen  in  these  systems. 

These  initial  studies  indicate  a  stable  Ar2F  molecule  in  the  excited  2B 


configuration  as  shown  in  Fig.  3.  Such  a  molecule  could  dissociate  symmetri¬ 
cally  to  Ar+  (2S+)  +  F“(1Sq)  or  unsummetrically  to  Ar+F"  (%+)  +  Ar  (lsD).  The 
proper  dissociation  (or  formation  path)  can  be  ascertained  by  more  detailed 
calculations  at  larger  internuclear  separations.  The  excited  2Aj_  configuration 
exhibits  stability  relative  to  Ar+F'  +  Ar,  but  not  to  Ar+  +  F”.  A  symmetric 

dissociation  leads  to  the  unbound  2Z+  state  of  Ar+  +  F"  ). 

S  2  o  * 

These  limited  calculations  suggest  that  the  upper  radiating  state  of  Ar2F 

'and  Kr^)  is  of  2B2  symmetry  with  equilibrium  bond  lengths  of  R,  <  2.6  A 
0 

and  Rat-F  <2.5  A.  Using  a  dissociation  energy  of  1.25  eV  for  X  2l£  of  Ar*  and 
our  calculated  ground  state  %,  repulsion  energy,  ve  predict  an  emission  wave¬ 
length  of  ~  280  rut.  The  observed  emission  for  Ar2F  is  very  broad  (~  50  nm)  and 
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peaks  at  about  290  ran. 

This  -»  %2  transition  corresponds  more  closely  to  the  C  -•  A  % 
transition  in  ArF  than  to  the  intense  -•  ^ £+  transition.  Scaling  our  cal¬ 
culated  transition  intensities  for  C  -*  A  we  predict  a  radiative  lifetime 
for  the  ionic  ^ B *  state  of  ~  200  nsec.  Further  detailed  studies  are  in  progress 
to  classify  the  state  assignments  and  predict  the  radiation  characteristics  of 
all  the  triatomic  noble  gas  halide  excimer  molecules. 
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TABLE  I.  MOLECULAR  ORBITAL  CORRELATIO 


2Ai  <R-,  =  3.18) 

-  2Ai  (Rt  =2.91) 


IONIC  __ 
Ar+  +  F“  +  Ar 


Ar+  (2P)  +  F~  ( ■  SQ)  +  Ar  0so) 


Ar2+  (2ZU+)  +  F~  (%) 


2Ai  (R-i  =  265) 


IONIC  Ar2+  +  F~ 


2B2  (R-,  =  291) 


2B2  (R-i  =  265) 


-  *2 
R1  - *  F 


2B2  (R-j  =  265) 


2Ai  (R-,  =2.65) 


Ar  (ISo)  +  Ar  (%)  +  F(2P) 


Ar2-F  SEPARATION,  R2(A) 

gure  3  Potential  Energy  Curves  for  ArjP 
(Density  Functional  Method) 
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APPENDIX  B 


THE  ELECTRONIC  STRUCTURE  OF  EXCIMER  MOLECULAR  LASERS* 


L.  A.  Wright  (DYP) 

Kirtland  Air  Force  Base,  New  Mexico  87117 


ABSTRACT 


The  electronic  structure,  potential  energy  curves  and  radiative  transition 
probabilities  of  excimer  systems  have  been  examined  using  quantum  mechanical 
methods.  These  molecules  are  characterized  by  repulsive  or  weakly  bound  ground 
state  potential  curves  and  by  bound,  strongly  ionic,  or  Rydberg,  excited  states 
They  constitute  a  very  interesting  class  of  molecules  which  offer  the  possibili¬ 
ty  for  high  power,  high  efficiency  UV  laser  operation.  Calculations  have  been 
carried  out  using  the  density  functional  SCF-Xa  method,  modified  extensively 
to  correct  for  well-known  errors  arising  from  approximations  to  the  potential 
and  exchange  terras.  A  limited  number  of  ab  initio  calculations  were  also  car¬ 
ried  out  for  comparison  purposes.  For  the  ArF  system  we  find  that  the  lowest 


•k 

^Supported  in  part  by  AFOSR  under  Contract  No.  F44620-73-C-0077 . 
Permanent  Address:  Quantum  Theory  Project 

University  of  Florida 
Gainesville,  Florida  32611 


59 


excited  ionic  state  has  symmetry  ^2+i/2>  anc*  that  the  dominant  laser  transition 
observed  at  1933  A  should  be  assigned  to  B  ^"Z  +  \/2  ♦  X  1/2*  The  C  II3/2  +  %■  2  1/2 

transition  is  calculated  to  be  two  orders  of  magnitude  smaller  in  emission  in¬ 
tensity  than  the  dominant  transition,  thus  ruling  out  this  assignment  for  the 
observed  laser  line  in  ArF .  Preliminary  calculations  carried  out  for  Ar2F 
indicate  that  the  bound  upper  ionic  state  has  2B2  symmetry  and  that  the  most 
probable  ground  state  also  has  2B2  symmetry.  This  polyatomic  system  is 
predicted  to  exhibit  a  broadband  emission  spectra  with  a  relatively  long  radi¬ 
ative  lifetime. 
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INTRODUCTION 


The  diatomic  molecules  formed  from  the  noble  gases,  from  the  partially 
closed  shell  Group  II  B  elements  (Hg,  Cd ,  Zn)  or  from  the  noble  gas-halogen 
combinations  constitute  an  interesting  and  important  class  of  molecules  for 
laser  applications.  These  molecules  are  characterized  by  a  repulsive  or  very 
weakly  bound  ground  state  potential  curve  and  by  bound  electronically  excited 
states.  These  excited  diatomic  molecules,  known  as  excimers,  radiate  in  a 
narrow  visible  or  UV  continuum  band  by  making  a  transition  to  a  repulsive  lower 
state  that  dissociates  to  ground  state  neutral  atoms.  Laser  behavior  is  con¬ 
trolled  by  competing  kinetic  processes  for  the  formation  and  destruction  of  the 
excited  molecular  levels.  The  character  of  the  electronic  structure  of  the 
excited  states  in  the  noble  gas  halide  systems  is  strongly  ionic  and  is  similar 
to  that  of  the  ground  state  alkali  halide  salts.  The  other  molecules  exhibit 
mainly  Rydberg  character  in  their  excited  electronic  states. 

The  system  Xe2,  shown  in  figure  1,  is  typical  of  the  noble  gas  excimer 
lasers.  Because  the  ground  state  potential  is  almost  purely  repuls ive ,' the 
emission  is  a  continuum  which  stretches  from  the  atomic  resonance  line  at  1470  A 

O 

to  2200  A.  This  emission  has  been  studied  extensively  (ref.  1)  and  has  been 
widely  used  as  a  background  source  for  VUV  absorption  spectroscopy.  Houtermans 
(ref.  2)  was  the  first  to  propose  the  use  of  bound-free  transitions  (in  Hg2)  for 
a  laser  system,  but  it  was  not  until  1970  that  Basov  (ref.  3)  demonstrated  gain 
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ac  1750  A  by  pumping  liquid  xenon  with  a  high  current  beam  of  electrons. 

Since  then  many  workers  have  reported  laser  action  in  the  gaseous  phase  in  xenon 
and  krypton  A  (ref.  4) . 

Emission  spectra  for  several  noble  gas  halides  have  been  obtained  in  low 
pressure  discharge  flow  experiments  and  were  first  reported  by  Velazco  and 
Setser  (ref.  5)  and  by  Golde  and  Thrush  (ref.  6).  More  recently,  emission 
spectra  and  laser  action  in  several  noble  gas  halide  systems  have  been  repor¬ 
ted  (refs.  7-18).  The  spectra  are  consistent  with  a  model  whereby  the  lasing 
transition  is  between  a  strongly  bound  upper  electronic  state  and  a  nearly  flat 
repulsive  ground  state.  A  group  theoretical  analysis  of  the  possible  electronic 
configurations  for  these  states  does  not  lead  to  a  unique  assignment  since  both 
the  upper  and  lower  levels  could  be  of  l  or  n  symmetry. 

In  addition  to  the  noble  gas  vi.olecules,  systems  such  as  Hg2,  HgCd,  and  HgTl 

have  been  proposed  (refs.  2,  19-21)  as  a  class  of  excimers  which  radiate  in  the 

visible  region  of  the  spectrum.  For  these  systems,  the  partially  filled  shell 
10  2 

of  Hg  (5d  6s*)  results  in  noble  gas-like  behavior  and  leads  to  potential 

curves  analogous  to  those  arising  in  Xe2 .  For  the  HgTl  system,  a  low-lying 
ionic  state  can  occur,  and  the  electronic  structure  is  probably  more  like  that 
of  the  noble  gas  halide  molecules. 

The  theoretical  principles  governing  the  electronic  structure  of  the 
homonuclear  noble  gas  molecules  have  been  discussed  by  Mulliken  (ref.  22)  using 
Xe2  as  a  prototype  and  more  recently  by  Barr,  Dee,  and  Gilmore  (ref.  23)  for  the 
case  of  Kr2 .  In  these  molecules,  the  ground  state  is  represented  by  a  repulsive 

X 


potential  curve  connecting  to  ‘so  states  of  the  separated  atoms.  The 
excited  states  constitute  a  Rydberg  series  terminating  on  the  molecular  ion  of 
symmetry  £  u  or  1/2U,  the  latter  being  more  appropriate  for  the  heavier 
systems  such  as  Xe2  where  L-S  coupling  is  no  longer  strictly  valid.  These 
excited  state  potential  curves  can  be  located  with  good  accuracy  using  quantum 
defect  theory  coupled  with  an  accurate  representation  of  the  ground  state  of 
the  molecular  ion.  the  radiative  decay  time  for  an  excited  Xe,  or  Kr, 
molecule  has  been  estimated  to  be  so  short  (  10  nsec),  that  fast  pumping  tech¬ 
niques  are  required  to  produce  significant  populations  in  their  excited  states. 

Theoretical  studies  of  the  noble  gas  halide  molecules  have  consisted  of 
analogies  with  alkali  halide  systems  (ref.  7),  predictions  based  on  an  electron- 
gas  model  (ref.  24)  and  ab  initio  calculations  for  the  KrF  system  (ref.  25). 

In  some  cases,  the  theoretical  methods  utilized  to  date  have  been  unable  to 
posit ively  define  the  electronic  states  involved  in  the  noble  gas  halide  laser 
systems .  Ab  initio  methods  would  appear  to  be  practical  only  for  the  lighter 
aystems  such  as  ArF  and  KrF  owing  to  greatly  increased  computational  complexity 
for  larger  numbers  of  electrons.  However,  detailed  ab  initio  Cl  calculations 
on  systems  as  complex  as  XeCl  are  now  in  progress  (ref.  26).  The  electron-gas 
model  calculations  suffer  fro.  the  lack  of  self-consistency  in  the  molecular 
potential.  Recent  density  functional  calculations  indicate  that  self-consis¬ 
tency  in  the  potential  is  crucial  to  determining  optimum  electronic  wavefunc- 
tions  (ref.  27).  Such  calculations,  coupled  with  a  limited  number  of  ab  initio 
calculations  carried  out  for  comparision  purposes,  constitute  the  principle 
effort  in  the  present  study  of  excimer  systems. 
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METHOD  OF  APPROACH 


Density  Functional  Approach  -  Xa  Model 

The  Xa  model  (ref.  28)  for  the  electronic  structure  of  atoms,  molecules, 
clusters  and  solids  is  a  local  potential  model  obtained  by  making  a  simple 
approximation  to  the  exchange-correlation  energy.  If  we  assume  a  nonrela- 
tivistic  Hamiltonian  with  only  electrostatic  interactions,  it  can  be  shown  that 
the  total  energy  E  of  a  system  can  be  written  exactly  (ref.  29)  (in  atomic 
units)  as 


E=  Z  n j  <  Uj 

i 


u.  >  + 


-r  I 


+  1. 

ij 


"i"j 


UiUi 


'12 


Uj  Uj  >  +EX( 


(1) 


This  expression  is  exact  provided  the  U£  are  natural  orbitals  and  nt 
are  their  occupation  numbers  (i.e.,  eigenfunctions  and  eigenvalues  of  the  first 
order  density  matrix).  The  first  term  in  equation  (1)  represents  the  kinetic 
and  e 1 ec t ron— nuc le ar  energies.  The  second  term  is  the  nuclear  repulsion  energy. 
The  sums  (u,  v)  are  over  all  the  nuclear  charges  in  the  system.  The  third  term 
is  the  electron-electron  repulsion  term,  which  represents  the  classical  electro¬ 
static  energy  of  the  charge  density  p  interacting  with  itself,  where 


^(l)  =1  n.u*(D  Uj(D 


(2) 
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The  last  term  pxc  (1,2)  represents  the  exchange-correlation  energy  and  can  be 
expressed  formally  as 


1  xC 


s  -Lfpl l)  dr7 /  fxc  (l^j  <j r2 

12 


(3) 


where  pxc  (1,2)  represents  the  exchange-correlat ion  hole  around  an 
electron  at  position  1.  In  the  exact  expression,  pxc  is  dependent  on  the 
second-order  density  matrix.  In  the  Hartree-Fock  approximation  Exc  is  the 
exchange  energy,  pxc  represents  the  Fermi  hole  due  to  the  exclusion  principle 
and  depends  only  on  the  first-order  density  matrix.  In  the  Xa  method,  we  make 
a  simpler  assumption  about  pxc .  If  we  assume  that  the  exchange-correlation 
hole  is  centered  on  the  electron  and  is  spherically  symmetric,  it  can  be  shown 
that  the  exchange-correlation  potential 


xc 


*/ 


<•  .<?;■»  drT 

'  12 


(4) 


is  inversely  proportional  to  the  range  of  the  hole,  rg ,  where  is  defined 


by 


47 T 
3 


p(  I)  =  I 


(5) 


1  /  Q 

Therefore,  in  the  Xa  model,  the  potential  Uxc  is  proportional  to  p  '  A 

scaling  parameter  a  is  defined  such  that 


UXa(D=  -  4?-  (3/>(l)/87r)l/3 


(6) 
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The  expression  m  equation  (6)  is  defined  so  that  a  =  2/3  for  the  case 


of  a  free-electron  gas  in  the  Hartree-Fock  model  (ref.  30)  and  a  =  1  for  the 
potential  originally  suggested  by  Slater  (ref.  31).  A  convenient  way  to  choos< 
this  parameter  for  molecular  and  solid  state  applications  is  to  optimize  the 
solutions  to  the  Xa  equations  in  the  atomic  limit.  Schwarz  (ref.  32)  has  done 
this  for  atoms  from  Z  =  1  to  Z  =  41  and  found  values  between  2/3  and  1. 
Computational  Aspects  of  the  Xa  Method 


To  implement  the  Xa  approximation  in  molecular  calculations,  it  is 
usually  convenient  to  break  up  the  volume  integrals  indicated  in  equations 
(3)  or  (4)  by  defining  spherical  regions  surrounding  each  atomic  center.  The 
molecular  potential  is  then  spherically  averaged  within  these  regions  and  as¬ 
sumed  constant  in  the  region  between  spheres.  The  potential  falls  to  zero  in 
the  usual  way  in  the  region  outside  of  an  overall  sphere  located  to  encompass 
those  spheres  surrounding  the  individual  atoms. 


There  are  two  other  practical  aspects  of  the  calculations  which  must  be 
considered  in  application  of  the  Xa  model  to  finite  molecular  systems.  The 
first  concerns  the  choice  of  the  integration  framework  for  describing  the 
molecular-orbital  wavefunct ions ,  and  the  second  deals  with  the  choice  of  the 
exchange  parameter,  a,  in  different  regions  of  space.  A  third  consideration, 
convergence  of  the  spherical  harmonic  expansion,  was  also  tested.  The  total 
energy  was  converged  when  all  terms  up  to  1  =  3  were  included. 
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In  computations  with  heteronuc lear  molecules,  there  are  several  free 
parameters  that  must  be  chosen:  the  ratio  of  sphere  radii  for  the  atomic 
spheres  of  integration  at  a  given  internuclear  separation,  the  degree  of  sphere 
overlap,  and  the  value  of  the  exchange  parameter  in  the  atomic  spheres  and  the 
intersphere  region. 

It  has  been  found  that  changing  the  ratio  of  the  sphere  radii  for  the  two 
atoms  in  a  heternuclear  diatomic  molecule  introduces  changes  in  the  total  energy 
chat  can  be  large  on  a  chemical  scale,  (~1  eV)  (ref.  27).  A  choice  for  sphere 
radii  based  on  covalent  bonding  radii  does  not  necessarily  provide  a  good  esti¬ 
mate  for  these  calculations,  the  value  of  the  exchange  parameter,  a,  and  the 
sphere  radii  and/or  sphere  overlap  is  normally  fixed  in  Xa  calculations  for 
crystals,  where  the  geometry  is  fixed.  However,  to  develop  a  potential  curve, 
the  molecular  description  needs  to  ch  inge  substantially  as  the  internuclear 
separation  varies  and  the  changing  sphere  radii  include  varying  fractions  of  the 
total  molecular  charge  (ref.  27).  Our  studies  have  shown  that  at  any  given 
separation,  the  total  energy  calculated  from  the  Xa  model  is  a  minimum  at  the 
radii  ratio  where  the  spherically  averaged  potentials  from  the  two  atomic 
centers  is  equal  at  the  sphere  radius, 

V,(rf|)- V2(rJ  (7/ 

This  relationship  between  the  potential  match  at  the  sphere  boundary  and  the 
minimum  in  the  total  energy  appears  to  hold  exactly  for  "neutral"  atoms  and 
holds  well  for  ionic  molecular  constituents.  In  the  case  of  two  ionic  species, 
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Che  long-range  tail  of  the  potential  must  go  like  +2/R  from  one  ion  and  -2/R  (in 
Rydbergs)  for  the  other  ion  and  so  at  large  internuclear  separations,  the  tails 
of  the  potential  cannot  match  well.  However,  at  intermediate  separations,  the 
1/R  character  of  the  potential  does  not  invalidate  the  potential  match  criterion 
for  radii  selection.  This  match  for  the  atomic  potentials  is  applied  to  the 
self-consistent  potentials. 

A  double  iteration  is  now  required  for  each  chosen  internuclear  geometry. 
Computationally,  an  initial  choice  for  the  sphere  radii  is  made,  and  the  calcu¬ 
lation  is  iterated  to  self-consistency  in  the  charge.  A  comparison  of  the 
values  of  the  potential  at  the  sphere  radii  is  then  made,  and  an  adjustment  in 
their  ratios  is  carried  out  until  equation  (7)  is  satisfied  to  within  a  few 


percent , 


In  molecules  with  significant  charge  sharing  in  the  bonds,  the  radii  of 
the  atomic  spheres  is  frequently  increased  in  Xa  calculations  so  that  an 
overlap  region  appears  in  the  vicinity  of  the  bond  (ref.  33).  Our  studies 
have  shown  that  the  contribution  to  the  total  molecular  energy  from  the  ex- 
change-correlation  term  shows  a  minimum  at  the  optimum  sphere  radius  or  sphere 
overlap.  This  provides  a  sensitive  criterion  for  selecting  these  parameters. 

The  values  of  the  exchange  parameters  in  the  spherical  integration  region 
around  each  atomic  center  are  frequently  set  at  the  atomic  values  both  for 
neutral  and  for  ionic  molecular  constituents.  However,  for  light  atoms,  the 
value  of  a  which  best  reproduces  Hartree  Fock  results  varies  substantially  with 
tonicity .  In  Ar  and  Ar+,  the  following  table  compares  the  HF  energy  and  the 
xa  energy  calculated  for  several  values  of  a. 
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a 

Xa  Energy  ( a. u. ) 

HF  Energy  (a. u. ) 

> 

f-1 

o 

.72177 

-526.8176 

-526.8173 

Ar*1'2 

.72177 

-526.5857 

- 

.72213 

-526.6007 

Ar*1 

.72177 

-526.2447 

- 

.72213 

-526.2596 

- 

.72249 

-526.2745 

-526.2743 

The  optimum  value  of  a  changes  even  more  rapidly  in  the  fluorine  atom,  going 
from  0.73732  for  F^  to  0.72991  for  F-^  .  Since  the  total  energy  depends 
linearly  on  a,  this  parameter  must  be  chosen  carefully. 

The  intersphere  exchange  coefficient  is  chosen  to  be  a  weighted  average  of 
the  atomic  exchange  parameters  from  the  two  constituents.  At  small  internuclear 
separations,  the  optimum  radius  for  an  atomic  sphere  frequently  places  signifi¬ 
cant  amounts  of  charge  outside  that  atomic  sphere  -  charge  that  is  still  strong¬ 
ly  associated  with  its  original  center  rather  than  being  transferred  to  the 
other  center  or  associated  with  the  molecular  binding  region.  To  best  account 
for  these  cases,  the  weighting  coefficients  are  chosen  to  reflect  the  origin  of 
the  charge  in  the  intersphere  (or  outersphere  region). 


^  inter, phere 


-Q?)+ «„(<?„ -Q?) 


(8) 


where  (Q  -  Q. °)  is  the  charge  lost  from  sphere  i  relative  to  its  atomic 
si  1 

value  (or  ionic  value  Q^°  and  ag^  is  the  atomic  exchange  parameter  for  sphere 
i.  This  value  for  aintersphere  *-s  calculated  dynamically  -  it  is  updated  after 
each  iteration  in  the  self-consistent  calculation. 
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While  for  heavy  atoms,  these  changes  in  the  exchange  parameter  would  be 
small,  the  a's  for  small  atoms  vary  rapidly  with  z  (and  with  ionicity).  The 
correct  choice  of  the  exchange  parameters  influences  not  only  the  total  energy 
calculated  for  the  molecule  but  also  affects  the  distribution  of  charge  between 
the  atomic  spheres  and  the  intersphere  region. 

Ab__Init_i£  Apj3roa£h_-_V£lenc_e_CI  Method 

The  ab  initio  calculation  procedure  chosen  for  these  studies  is  the 
valence-configuration-interaction  (VCI)  method  (refs.  34,  35).  A  spin- 
free  nonrelat ivistic  electrostatic  Hamiltonian  is  employed  within  the  Born- 
Oppenheimer  approximation.  For  a  diatomic  molecule,  this  approximation 
leads  to  an  electronic  Hamiltonian  depending  parametrically  on  the  inter- 
nuclear  separation  R.  Electronic  wave  functions  ip (R)  are  made  to  be  opti¬ 
mum  approximations  to  solutions  of  the  Schrodinger  equation 

H(R)<J,(R)  =  E(R)'HR)  (9) 

by  invoking  the  variational  principle.  The  specific  form  for  ip ( R)  may  be 
written  as 


=  I  cu^(R) 


(10) 


where  each  ip(R)  is  referred  to  as  a  configuration  and  has  the  general  structure 


n 


i l>*(R)  =  siCt  II  E)6m, 


(11) 
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where  <1>U£  is  a  spatial  orbital,  //is  the  ant isymmetrizing  operator,  is 

the  spin-projection  operator  for  spin  quantum  number  S,  and  0M  is  a  product 

ns 

of  a  and  g  one— electron  spin  functions  of  magnetic  quantum  number,  Mg .  The 
spatial  orbitals,  l/^,  have  the  general  form 

<t»  =  exp  {-<5(£ - (lV  +  (COS  0)  (12) 

where  for  —  £<5^,  the  standard  STO  form  is  required.  If  no  restriction  is 
imposed  as  to  the  double  occupancy  of  the  spatial  orbitals,  equations  (10)  and 
(11)  can  describe  a  completely  general  wave  function. 

The  basis  functions  chosen  for  describing  the  molecular  system  ArF  were 
opimized  elliptic  orbitals.  A  minimum  basis  was  chosen  for  the  K  and  L  inner 
shells,  and  both  minimum  and  elliptic  basis  were  used  to  describe  the  valence 
electrons.  A  typical  basis  for  R  =  3.5  a.u.  is  given  in  Table  I.  Configura¬ 
tion  interaction  was  taken  only  over  the  valence-shell  electrons.  Such  a  wave 
function  gives  an  accurate  representation  of  the  interatomic  potentials  at  large 
internuclear  separations  but  must  become  poorer  at  short  separations  where  L— 
shell  polarization  effects  are  important.  Both  o-  and  r-  coupled  valence  con¬ 
figurations  were  included  as  well  as  ionic  configurations.  The  configuration 
interaction  constituted  a  complete  Cl  over  the  valence  shell  electrons. 
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RESULTS  AND  DISCUSSION 


We  have  carried  out  both  detailed  ab  initio  and  density  functional  Xa 
calculations  for  the  ArF  and  A^F  systems  to  assess  problems,  both  conceptual 
and  computational,  in  defining  the  electronic  structure  of  the  noble  gas  halides. 

A  general  correlation  diagram  for  the  diatomic  systems  is  shown  in  Fig.  2. 
Here  we  have  indicated  the  detailed  coupling  of  atomic  states  to  form  the  ground 
(mainly  repulsive)  molecular  levels,  the  coupling  of  the  ionic  states  R+  and 
X-  to  form  the  upper  ionic  molecular  levels  and  the  very  complicated  excited 
Rydberg  state  structure  for  these  systems.  Only  the  lower-lying  Rydberg  states 
that  arise  from  excitation  in  the  noble  gas  or  excitation  in  the  corresponding 
halide  have  been  included  in  this  correlation  diagram.  The  left-hand  side 
corresponds  to  pure  (A,S)  coupling  which  is  mainly  applicable  to  light  systems 
such  as  ArF.  The  right-hand  side  corresponds  to  pure  (J,J)  coupling  which 
probably  is  never  quite  realized,  even  for  a  system  as  heavy  as  Xel.  Unfortun¬ 
ately,  intermediate  coupling  is  indicated  for  the  majority  of  these  noble  gas 
halide  systems  and  some  studies  of  molecular  spin-orbit  coupling  parameters  are 
clearly  indicated.  The  correlation  diagram  for  the  triatomic  system  Ar^F  is 
shown  in  Table  I  for  C2V  symmetry.  This  symmetry  dissociates  to  R2  +  F 
in  a  symmetrical  configuration.  Studies  of  the  potential  energy  surfaces  for 
a  lesser  symmetry  may  be  needed  to  examine  the  dissociation  channel  R+X  +  R. 

For  all  of  the  noble  gas  halide  systems  there  is  an  outer  crossing  of  the 
upper  ionic  states  arising  from  R+  +  X  (or  R£+  )  (R  =  noble  gas,  x  =  halide) 
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that  these  upper  molecular  states,  at  their  equilibrium  internuclear 

actually  Ue  lover  than  .„y  of  the  Rydbefg  _UcuUr  ^  ^ 

eX<"Ple-  the  l°“eSt  8>,db“S  *>'  ArF  arises  fr0„  Ar  3p(4s)  ,  p  2p  <t 

11.55  eV.  The  ionic  states  of  Ar+  +  F~  at  12  31 

at  12.31  eV  thus  cross  these  Rydberg 

states  at  R  =  36  bohr  radii. 

Possible  diatomic  electronic  transic.ons  uh.ch  correspoRd  ^  user 
activity  in  these  system  therefore  reduce  to  the  following  set: 


B  2Z+  -  x  2I+ 


B  2I  +  -  A  2n 


c  2n 


x  2z+ 


c  2n 


a  2n 


(13) 


Further  the  C  2H  state  is  spiit  into  2n3/2  a„d  ^  states  owing  to 

•Pin-orhit  ooupiing  effects  in  the  nohie  gas  positive  ion.  This  splitting  is 
small  (782  cm  1)  for  Ne  but  is  appreciable  for  Xe  (10537  cm~l). 

As  a  result  of  these  several  possibilities  a  definit 

J.L  ies ,  a  definite  assignment  of 

Che  laser  transitions  becomes  difficult.  Brau  and  E.ing  (refs.  7,  8)  have 
assigned  the  strong  laser  transitions  in  XeBr,  XeCl,  XeF  and  KrF  to 

1  1/2  *  X  ^  1/2  ^  aCC°rd  With  a  simPle  conic  model  for  the  potential 

curves.  Ab  initio  calculat ions  by  Dunning  and  Hay  (ref.  25)  for  KrF  also 

suggest  the  B  2I+.  +  v  2r  +  .  . 

1/2  1/2  trans it  ion  corresponds  to  the  observed 

emission  in  KrF  at  tu  • 

2480A.  Their  potential  curves,  however,  also  indicate 

a  C  2"3/2  *  X2  I  *1/2  transition  at  nearly  the  same  calculated  wavelength. 

Further,  their  calculations  indicate  a  near  degeneracy  of  the  upper  ionic 

1/2  "d  C2"3/2  StSteS  “hich  is  difficult  to  understand  in  terms  of 
simple  charge  density  models  for  these  excited  states. 
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Clugston  and  Gordon  (ref  24)  suggest  that  the  dominant  laser  transition 

*22+  • 

is  C  II3/2  ♦XI  i/2>  based  on  free-electron  gas  model  calculations.  Recent 

studies,  however,  question  the  reliability  of  non-SCF  model  calculations  for 
these  systems  (ref.  27).  In  addition,  this  assignment  of  a  perpendicular 
transition  to  the  laser  emission  goes  against  the  experimental  observation 
that  such  transitions  should  be  weak  compared  to  allowed  parallel  transitions. 

In  order  to  shed  some  light  on  the  above  questions,  both  ^b  initio  and 
density  functional  calculations  were  carried  out  for  the  lower  and  excited 
levels  of  ArF  and  A^F.  The  ab  initio  calculations  were  carried  out  in  a 
manner  similar  to  that  described  previously  (ref. 34).  All  nonlinear  param¬ 
eters  of  the  molecular  orbital  basis  functions  were  individually  optimized 
for  each  separate  electronic  state.  The  resultant  ab  initio  potential  curves 

are  shown  in  Fig.  3  and  Table  III  where  we  see  that,  unlike  the  previously  re- 

2  +  2 

ported  ab  initio  studies  of  KrF  (ref.  25),  the  B  and  C  U3/2  systems  are 

2 

non-degenerate.  The  upper  C  II3/2  state  1S  less  strongly  polarized  than  the 

2  +  •  •  •  ... 

B  I  y/2  state  and  exhibits  a  somewhat  flatter  potential  curve  giving  rise 

2  +  2  0 
to  an  inner  crossing  of  B  I  ^2  anc*  C  “3/2  at  a'3°ut  1*8  A. 

These  ab  initio  calculations  indicate  the  following  long-range  behavior 

for  the  excited  state  potential  curves: 


B  2I+:  V(R)( a.u.)=  --- ^ 
R  R 

C2n:  V(tf)(a.u.)=  -1+2^8 
R  R 


(14) 
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The  R"3  term  arises,  for  these  molecules,  mainly  from  the  quadrupole  contribu¬ 
tion  to  the  interaction  potential  at  long  range.  This  term  vanishes  for  the 
alkali  halide  diatomic  molecules  M+X"  since  both  positive  and  negative  ions  are 
in  the  h  state  with  a  non-vanishing  quadrupole  moment.  Therefore,  simple 
models  of  the  ionic  excited  states  of  the  noble  gas  halides,  based  on  analogies 
with  alkali  halide  systems,  need  to  be  modified  to  properly  account  for  this 
long  range  behavior.  The  calculated  value  for  the  coefficient,  3 ,  of  the  R 
term  has  the  correct  theoretical  form  such  that 

0(22*)=-2/3(2ri)  (15) 

Calculated  spontaneous  emission  spectra,  without  inclusion  of  spin  orbit 
Z-H  mixing,  are  shown  in  Fig.  4.  They  are  in  good  agreement  with  previously 
reported  emission  spectra  (ref.  10)  when  the  state  mixing  effects  are  properly 
accounted  for  and  suggest  a  mild  revision  in  the  ground  state  at  the  minimum  of 
the  excited  potential  curve  since  the  half-width  of  the  strong  transition  is 

o 

about  80  A. 

Although  several  system  assignments  are  possible  within  the  calculated 
error  of  the  potential  curves,  Fig.  4  clearly  indicates  that  only  the  B  2I 
-  X  V  transition  exhibits  a  large  transition  moment  and  hence  a  short  (~  5 
nsec)  radiative  lifetime.  Thus  the  original  assignment  of  Ewing  and  Brau  of  the 
B  V  -  xV  transition  to  the  strong  laser  line  in  these  noble  gas  halide 
systems  appears  to  be  confirmed  by  these  detailed  calculations. 

As  a  comparision,  potential  energy  curves  for  ArF  were  also  constructed 
using  a  density  functional  method.  As  described  above,  our  program  has  been 


76 


extensively  developed  beyond  the  standard  techniques  described  in  the  litera¬ 
ture.  Modifications  to  minimize  the  muffin-tin  potential  error  and  to  insure 
continuity  in  the  molecular  potential  function  have  been  incorporated  into  our 
code  which  results  in  more  reliable  total  energies.  The  results  for  ArF  are 
shown  in  Fig.  5  and  Table  IV  which  indicate  good  agreement  with  our  ab  initio 

results.  Again  we  see  that  the  C  excited  state  is  nondegenerate  with  the 
2  + 

upper  B  Z  state  and  that  it  is  somewhat  flatter  in  curvature.  The  ori¬ 
gin  of  the  somewhat  steeper  repulsion  exhibited  by  all  the  calculated  poten¬ 
tials  using  this  density  functional  model  is  now  understood  as  a  manifestation 
of  the  exchange  scaling  parameter  variation  with  the  effective  nuclear  charge. 
Thus  a  in  Eq.  (6)  should  be  modified  at  each  internuclear  separation  to  account 
for  internal  changes  from  Ar-F  to  Ar  -  F  By  introducing  a  concept  which 
we  call  "dynamic  scaling"  into  the  density  functional  calculations,  our  most  re¬ 
cent  data  for  ArF  are  now  almost  in  complete  agreement  with  the  ab  initio  re¬ 
sults.  The  procedure  is  somewhat  akin  to  uniform  scaling  of  a  wavefunction  in 
an  ab  initio  calculation.  The  significant  point  here  is  that  this  density 
functional  method  can  be  applied  to  a  system  such  as  Kr2F  with  only  about 
2-3  times  the  effort  that  is  required  for  KrF.  In  contrast,  ab  initio  calcu¬ 
lations  for  Kr2F  would  constitute  an  enormous  (>  100)  increase  in  effort  over 
that  required  for  a  study  of  KrF  of  comparable  accuracy.  The  density  functional 
method  thus  offers  promise  of  giving  a  uniform  description  of  the  electronic 
structure  of  all  the  viable  candidates  for  noble  gas  halide  laser  systems. 


Recently,  regions  of  the  potential  energy  surfaces  for  Ar2F  have  been 
studied  to  examine  the  possible  stability  of  this  triatomic  molecule  in  an 
excited  electronic  configuration.  It  has  been  suggested  (ref.  36)  that  the 
long  wavelength  radiation  seen  at  2900  A  in  the  e-beam  excited  argon-fluorine 
mixtures  (4100  A  for  krypton-fluorine  mixtures)  may  be  due  to  a  triatomic 
molecule  formed  in  an  excited  ionic  state  from  Ar2+  +  F",  similar  to  the 
excited  ionic  state  formed  from  Ar+  +  F~.  Such  a  molecule  would  radiate  at 
longer  wavelength  since  Ar2+  is  stable  relative  to  Ar+  +  Ar. 

The  molecular  correlation  diagram  for  Ar2F ,  given  in  Table  I  for  C2v 
symmetry,  should  represent  the  most  stable  geometry  for  both  the  excited  ionic 
and  ground  state  surfaces.  The  corresponding  orbital  designations  are  also 
given  for  the  linear  molecule  Ar-F-Ar.  We  have  carried  out  preliminary  calcula 
Cions  for  the  ground  state  in  \  and  %  symmetries  and  for  excited  electronic 
states  in  2B2  and  \  symmetry.  Our  calculated  results  are  shown  in  Fig.  6. 

We  find  a  flat  potential  surface  for  Ar2F  in  the  lowest  2Al  state  for 


R2  >  2.8  A  as  shown  in  Fig.  6.  This  results  from  forward  charge  transfer 
from  Ar  to  the  F  atom,  thereby  introducing  some  ionic  character  of  the  type 
Ar2  F  x  which  will  tend  to  cancel  out  neutral  atom  exchange  and  overlap 

effects.  For  heavier  systems  such  as  Xe2F,  some  small  binding  may  be  found 

for  this  lowest  \  configuration.  The  lowest  2B2  state  exhibits  a  some¬ 
what  steeper  repulsion. 

These  initial  studies  indicate  a  stable  Ar£F  molecule  in  the  excited 
2B2  configuration  as  shown  in  Fig.  6.  Such  a  molecule  could  dissociate 


symmetrically  to  Ar+2  (22u+)  +  F"(1Sq)  or  unsymmetrically  to  Ar+F~  (22  +  ) 

+  Ar  (■LS0)The  proper  dissociation  (or  formation  path)  can  be  ascertained  by 
more  detailed  calculations  at  larger  internuclear  separations.  The  excited  2A^ 
configuration  exhibits  stability  relative  to  Ar+F~  +  Ar,  as  shown  in  Fig.  6, 
but  not  to  Ar2+  +  F  .  A  symmetric  dissociation  leads  to  the  unbound  22 

g 

state  of  Ar2+  *  F~  (^SQ). 

These  limited  calculations  suggest  that  the  upper  radiating  state  of 
2 

Ar2F  (and  Kr2F)  is  of  B2  symmetry  with  equilibrium  bond  lengths  of 

•  o 

RAr-Ar  ^  2‘^  A  and  RAr-F  ^  2-^  A-  Using  a  dissociation  energy  of  1.25  eV 

for  X2^u+  of  Ar2+  and  our  calculated  ground  state  2B2  repulsion  energy, 

we  predict  an  emission  wavelength  of  ~  280  nm.  The  observed  emission  for  Ar2F 

is  very  broad  ( A  ~  50  nm)  and  peaks  at  about  290  nm. 

This  B2  B2  transition  corresponds  more  closely  to  the  C  * 

A  2H  transition  in  ArF  than  to  the  intense  2z+  2S+  transition.  Scaling  our 

calculated  transition  intensities  for  C  2n  +  A  2jI,  we  predict  a  radiative  life- 
•  .  .  9  * 

time  for  the  ionic  B2  state  of  ~  200  nsec.  Further  detailed  studies  are  in 
progress  to  classify  the  state  assignments  and  predict  the  radiation  characteris¬ 
tics  of  the  triatomic  noble  gas  halide  excimer  molecules. 
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TABLE  I 


OPTIMIZED 

ELLIPTIC 

0REITALS 

FOR  ArF 

AT  R  =  3. 

.5  a .u . 

Ar 

Vd) 

A 

A 

Van 

A 

A 

2nU) 

A 

A 

2rLWl 

A 

A 

3s* 

2.5635 

2.6390 

2.5795 

2.7532 

2.5519 

2.6489 

2.5853 

2.7706 

3Pl 

2.2528 

2.2619 

2.3049 

2.3731 

2.2500 

2.2372 

2.3114 

2.4049 

3Po 

2.2728 

2.3074 

2.2763 

2.5828 

2.2520 

2.3632 

2.3014 

2.5164 

2s 

2.5507 

2.5807 

2.5034 

2.4979 

2.5544 

2.5783 

2.4933 

2.4963 

2Pl 

2.5053 

2.5498 

2.3697 

2.3624 

2.5277 

2.5713 

2.3379 

2.3275 

2P0 

2.5996 

2.6163 

2.4338 

2.4254 

2.5562 

2.5957 

2.4198 

2.3920 

*  The  approximate  STO  form  is  used. here  to  label  these  elliptic  functions. 
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I 


Symmetry 


TABLE  II 

MOLECULAR  ORBITAL  CORRELATION  DIAGRAM  FOR  Ar2F 
Orbitals  (Valence  Shell  Elec t r ons ) 


Ar 

Ar 

F 

Ar 

Ar 

Ar 

Ar 

Ar 

Ar 

F 

F 

F 

Atomic 

3s 

3s 

2s 

3PZ 

3py 

3Px 

3py 

3Px 

3PZ 

2Px 

2PZ 

2p. 

D«°h 

lag 

lau 

3°g 

llru 

liru 

lirg 

lTg 

2au 

2*u 

3au 

2lr, 

C2v 

lal 

lb2 

2ai 

3al 

4ai 

Ibl 

2b  2 

1  c&2 

3b  2 

2bl 

4b  2 

5a 

TABLE  III 


Ab  INITIO  POTENTIAL  CURVES  FOR  ArF  (INCLUDING  SPIN-ORBIT  CORRECTION) 

Energy  (eV) 


R(a.u. ) 

Y2  + 

.2 

'2 

D2  + 

9 

2 

x  2 

A  II3/2 

A  H1/2 

B  2  i/2 

c  nj/2 

D  m/2 

3.50 

1.8635 

3.1808 

3.2309 

7.4317 

7.5699 

7.7474 

3.75 

1.0104 

1.8191 

1.8692 

6.4926 

6.7506 

6.9281 

4.00 

0.5236 

1.0125 

1.0626 

6.0937 

6.428A 

6.6059 

4.25 

0.2474 

0.5393 

0.5894 

6.0194 

6.3805 

6.5580 

4.50 

0.0980 

0.2718 

0.3219 

6.1318 

6.4899 

6.6674 

4.75 

0.0201 

0.1238 

0.1739 

6.3381 

6.6771 

6.8546 

5.00 

0.0185 

0.0449 

0.0950 

6.5854 

6.8970  , 

7.0745 

6.00 

0.0022 

0.0049 

0.0550 

7.6301 

7.8469 

8.0244 

7.00 

0.0 

0.0003 

0.0504 

8.3294 

8.4655 

8.6430 

8.00 

0.0 

0.0 

0.0501 

8.8472 

8.9373 

9.1148 

9.00 

-0.0003 

-0.0003 

0.0498 

9.2440 

9.3068 

9.4843 

10.00 

-0.0008 

-0.0008 

0.0493 

9.5580 

9.6035 

9.7810 

00 

0.0 

0.0 

0.0501 

12.311 

12.311 

12.488 

TABLE  IV 


DENSITY  FUNCTIONAL  POTENTIAL  ENERGY  CURVES 


Energy  (eV) 


R(  a. u. ) 

X2I+(I) 

3.75 

4.136 

4.00 

2.095 

4.25 

1.075 

4.50 

0.599 

4.75 

0.3  26 

5.00 

0.259 

5.25 

0.191 

5.50 

0. 163 

5.75 

— 

6.00 

0.136 

6.50 

— 

7.00 

0.122 

8.00 

0.122 

9.00 

0.  1  22 

CD 

0.0 

A2n(i)  bV(h) 


— 

11.477 

3.932 

8.076 

2.299 

7.123 

1.415 

6.715 

0.844 

6.647 

0.571 

6.715 

0.422 

6.851 

0.299 

7.055 

0.  272 

7.259 

0.163 

7.504 

— 

7.912 

0.122 

8.348 

0.122 

8.960 

0.122 

9.463 

0.0 

12.375 

FOR  ArF 


c2n ( ii) 


7.395 

7.082 

6.974 

6.987 

7.083 

7.232 

7.422 

7.572 

8.076 

8.484 

9.042 

9.504 

12.375 


R+  +  X- 


R  +  X 


(A,S) 


(J"|, 


UNIQUE  ASSIGNMENTS  FROM  (Jj,  J2)  TO  (A,  S)  COUPLING 
ARE  INDICATED  BY  SOLID  LINES 


Figure  2  -  General  Correlation  Diagram  for  Intermediate  Coupling  of  Noble-gas 
Halide  Diatomic  Molecules  Including  Low-Lying  Rydberg  States 
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Figure  3  -  Pote 


Ar  +  F 

's0+  2po1/2 
'S0+  2po3/ 


'2 


2Ai  (Ri  =3.18} 

-  2Al  (Ri  =2.91) 


IONIC  _ 
Ar+  +  F“  +  Ar 


2Ai  (R,  =  2.65) 


Ar+  (2P)  +  F-  (1Sq)  +  Ar  (“*S0) 


Ar2+  (22u+)  +  F~  0so) 


IONIC  Ar2+  +  F~ 


2B2  (Rt  =  2.91) 


2B2  (Rt  =  2.65) 


2B2  (Rt  =  Z65) 


r2 

R1  - - - •  F 


2Ai  (Ri  =2.65) 


Ar  (ISq)  +  Ar  (1s0)  +  F(2p) 


Ar2-F  SEPARATION,  R2(A) 

Figure  6  -  Potential  Energy  Curves  for  Ar?F  (Densitv  F  , 

2f  tensity  Functional  Method) 
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APPENDIX  C 


Electronic  Structure  of  the  Noble  Gas  Dimer  Ions* 


L.  A.  Wright 

Kirtland  Air  Force  Base,  New  Mexico  87117 

ABSTRACT 

A  systematic  study  of  the  electronic  structure  and  chemical  binding  in  the 
dimer  ion  sequence,  Ne2  ,  A^  ,  and  Xe2+,  has  been  carried  out  using  both 

density  functional  and  ab  initio  configuration  interaction  computational  approaches. 
This  study  includes  detailed  calculations  of  the  pertinent  potential  energy 
curves  and  an  analysis  of  the  calculated  spectroscopic  properties  of  the  bound 
^  Su  (1/2^1)  state  of  these  ions. 

Supported  in  part  by  ?ARPA  under  ARPA  Order  No.  3432  under  AFWL  Contract 
F29601-78-C-0031  and  AFOSR  under  Contract  F49620-77-C-0064. 
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-  i.i  ..  . 


INTRODUCTION 


The  diatomic  molecules  formed  from  the  noble  gases,  from  the  partially  closed 
shell  ^roup  IIA  and  I IB  elements  or  from  the  noble  gas-halogen  combinations  con¬ 
stitute  an  interesting  and  important  class  of  molecules  for  laser  applications. 
These  molecules  are  characterized  by  a  repulsive  or  weakly  bound  ground  state 
potential  curve  and  by  bound,  electronically  excited  states.  These  excited  diatomic 
molecules,  known  as  excimers  ,  radiate  in  a  narrow  visible  or  uv  continuum  band  by 
making  a  transition  to  a  repulsive  lower  state  that  dissociates  into  ground  state 
neutral  atoms.  The  resulting  laser's  operating  characteristics  are  controlled  by 
kinetic  processes  competing  for  the  formation  and  destruction  of  the  excited 
molecular  levels  and  by  absorption  losses  arising  from  noble  gas  dimer  ions  and 
high-lying  Rydberg  states  of  the  excimer  molecules.  This  last  observation  suggested 
the  need  for  the  present  study  of  the  electronic  structure  and  absorption  character¬ 
istics  of  the  noble  gas  dimers  which  rate  studies  show  to  be  present  in  large  con¬ 
centrations  in  the  candidate  noble  gas  excimer  laser  systems. 

1-1+ 

Recent  experimental  studies  have  drawn  attention  to  uhotodissociation  of  the 
noble  gas  dimer  ions  as  an  important  loss  mechanism  in  excimer  lasers.  Hunter1 
has  measured  the  absorption  cross-section  profiles  for  the  strong  1/2  I-*l/2g  II 
(a  Eu  -*■  D  Ig  )  transition  in  both  Ar2+  and  Kr2  .  He  finds  very  strong, 
broad  band  absorption  peaking  at  292  nm  for  Ar2+  and  at  325  nm  for  Kr2+  .  These 
wavelengths  are  near  the  operating  region  of  several  excimer  lasers:  XeBr  (282  nm), 
XeCl  (308  nm),  XeF  (353  nm) ;  thus  providing  seme  concern  for  possible  inherent 
low  efficiency  operating  characteristics  for  these  systems. 

Theoretical,  treatments  of  the  noble  gas  dimer  ions  have  consisted  mainly  of 
ab  initio  SCF  studies  of  the  ground  A  X^  state  and,  by  far,  most  of  these 
studies  have  been  devoted  to  He2  .  We  shall  consider  here  only  the  systems^ 

Ne2  ,  Ar2+,  Xr2  and  Xe^ which  are  chemically  similar,  with  binding  arising  from 
the  outer  shell  n-electrons . 


■4-  +  cr 

The  Ne2  and  Arg  ions  were  first  studied  by  Gilbert  and  Wahl"  using  SCF  wavefunc- 

tions  constructed  from  STO  basis  functions.  These  authors  have  also  retorted 

2  2  2  +  s. 

.Limited  SCF  studies  of  the  B  /7  ,  C  /7  and  D  UL  excited  states  of  Xe0+  and 

S  U  g  2 

Ar2 •  studies  of  the  ground  state  of  the  Ne2  ion  were  retorted  bv 

T  8 

Cohen  and  Schneider '  ’  as  part  of  a  study  of  the  Rydberg  states  of  this  system. 
bCF  calculations  for  Ar2  at  short  internuclear  separations  have  been  reported  by 
Sidis,  et  al  and  more  recently  by  Stevens,  et  al  ,  using  an  extensive  STO  basis. 

4*  4“  j. 

Finally,  a  uniform  study  of  the  systems  Arg  ,  Kr2  and  Xe2  has  been  retorted  by 
11  12 

Wadt  who  employed  a  POL  Cl  using  gaussian  basis  functions. 

The  theoretical  principles  underlying  the  electronic  structure  of  the  homo- 

nuclear  noble  gas  molecular  ions  have  been  discussed  by  Mulliken1^  using  Xe2+ 

as  a  prototype  and,  more  recently,  by  Gilmore,  Barr  and  Deell+  for  the  case  of 

Kr„  .  For  these  dimer  ions,  four  electronic  states  (A 2£  +  ,  B2/7  ,  C2/7  and  D2JC  +  ' 

i  u  g’  u  g 

arise  from  a  ground  state  atom  (  SQ)  and  a  ground  state  ion  (2F^,3/2).  The  lowest 

2  + 

state,  usually  labeled  A  ,  is  bound  for  the  four  systems  Ne2+  ,  Ar2+,  Kr2+ 

and  Xe2  and  has  been  partially  characterized  experimentally  through  scattering 

experiments, although  there  has  been  much  controversy  concerning  the  dissociation 

energies  of  these  ions.  For  the  heavier  ions, A  -  s  coupling  is  no  longer  strictly 

valid  and  the  B  II  and  C  states  split  into  J— J  coupled  components  with  a 

resultant  mixing  of  E  and  n  symmetries.  This  splitting  is  small  for  Ne^  but  six 

distinct  J-J  coupled  molecular  ion  states  arise  for  Xe2+  . 

With  the  exception  of  the  recent  calculations  by  Wadt,11  no  attempts  at  systematic 

studies  of  these  ions  have  been  reported.  Accurate  ab  initio  methods  would  appear 

to  be  practical  only  for  the  lighter  ions  with  diminished  accuracy  for  heavier 

systems  such  as  Xe2+ .  Wadt  reports  such  a  trend  in  his  calculations11  which 

indicate  that  his  gaussian  basis  is  less  flexible  in  the  case  of  Xe^+  than  for  the 

corresponding  study  of  Ar2  .  Recent  density  function  calculations  of  noble  gas 

15 

halide  molecules^’  have  yielded  reliable  predictions  of  geometry  and  chemical 
binding,  provided  certain  optimization  criteria  are  followed.  Such  calculations,  which 
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should  be  of  more  uniform  quality  than  ab  initio  studies,  constitute  the  principle 
effort  in  the  present  study.  A  limited  number  of  ab  initio  calculations  were 
carried  out  for  Ar2+  for  comparison  purposes.  The  overall  agreement  between  the 
two  methods  is  excellent. 


Method  of  Approach 

Density  Functional_A£proach  -  XjX_Mo_del_ 

17 

The  Xa  model  for  the  electronic  structure  of  atoms,  molecules,  clusters 
and  solids  is  a  local  potential  model  obtained  by  making  a  simple  approximation 
to  the  exchange-correlation  energy .  If  we  assume  a  non-relativistic  Hamiltonian 
with  only  electrostatic  interactions  ,  it  can  be  shown  that  the  total  energy 

1  O 

E  of  a  system  can  be  written  exactly  (in  atomic  units)  as 


E=  I  n,  <  Uj  -mg-  V,2  +  £  7^- 
1  I  V 


u,  >  +  +  I 

TfJ.V 


-rfj  ni"i<uiuil^h''i>+E.c 


(1) 


This  expression  is  exact  provided  the  u.  are  natural  orbitals  and  n  are 

1  i 

their  occupation  numbers  (i.e.,  eigenfunctions  and  eigenvalues  of  the  first  order 
density  matrix).  The  first  term  in  Eq.  (l)  represents  the  kinetic  and  electron- 
nuclear  energies.  The  second  term  is  the  nuclear  repulsion  energv.  The  sums 
(#,!')  are  over  all  the  nuclear  charges  in  the  system.  The  third  term  is  the  electron- 
electron  repulsion  term,  which  represents  the  classical  electrostatic  energy 
of  the  charge  density  p  interacting  with  itself,  where 


(2) 
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The  last  term  E  represents  the  exchange-correlation  energy  and  can  be  ex- 
xc  * 

pressed  formally  as 


e,c  di T/  5* 

ri2 


(3) 


vhere  p  ^  ( 1,2)  represents  the  exchange-correlation  hole  around  an  electron  at 

position  1.  In  the  exact  expression,  Pxc  is  dependent  on  the  second-order  density 

matrix.  In  the  Hartree-Fock  approximation  E  is  the  exchange  energy,  p  rep- 

xc  xc 

resents  the  Fermi  hole  due  to  the  exclusion  principle  and  depends  only  on  the  first- 
order  density  matrix.  In  the  Xa  method,  a  simpler  assumption  about  p ^  is 
made.  If  we  assume  that  the  exchange-correlation  hole  is  centered  on  the  electron 
and  is  spherically  symmetric,  it  can  be  shown  that  the  exchange-correlation  potential 


_  r  Px cC',2) 

*c  -J  r!2 


(b) 


is  inversely  proportional  to  the  range  of  the  hole,  r  ,  where  r  is  defined  by 

s  s 


47 r  ^ 

3  r  s  pd  =  1 

(5) 

1  /o 

Therefore,  in  the  Xa  model,  the  potential  Uxc  is  proportional  to  p  (r).  A 
scaling  parameter  Ct  is  defined  such  that 

Uxa^=  T~  ( 3/>(l)/87r)l/3  (6) 
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The  expression  in  Eq.  6  is  defined  so  that  a  =  2/3  for  the  case  of  a  free 

19 

electron  gas  in  the  Hartree-Fock  model  and  Ot  =  1  for  the  potential  originally 

20 

sugEested  by  Slater  .  A  convenient  way  to  choose  this  parameter  for  molecular 

and  solid  state  applications  is  to  optimize  the  solutions  to  the  Xa  eauations  in 

21 

the  atomic  limit.  Schwarz  has  done  this  for  atoms  from  Z  =  1  to  Z  =  kl  and 
found  values  between  2/3  and  1. 

Computational  Aspects.  £f_the_Xa_  Method 

To  implement  the  Xa  approximation  in  molecular  calculations  it  is  usually 
convenient  to  break  up  the  volume  integrals  indicated  in  Eqs.  (3)  or  (U) 
by  defining  spherical  regions  surrounding  each  atomic  center.  The  molecular  po¬ 
tential  is  than  spherically  averaged  within  these  regions  and  assumed  constant 
in  the  region  between  spheres.  The  potential  falls  to  zero  in  the  usual  way  in 
the  region  outside  of  an  overall  sphere  located  to  encompass  those  spheres  surrounding 
the  individual  atoms . 

There  are  two  other  practical  aspects  of  the  calculations  which  must  be 
considered  in  application  of  the  Xa  model  to  finite  molecular  systems.  The  first 
concerns  the  choice  of  the  integration  framework  for  describing  the  molecular- 
orbital  wavefunctions  and  the  second  deals  with  the  choice  of  the  exchange  par¬ 
ameter,  Ct  ,  in  different  regions  of  space.  A  third  consideration,  convergence  of  the 
spherical  harmonic  expansion,  was  also  tested.  The  total  energy  was  converged 
when  all  terms  up  to  l  =3  were  included.  In  molecules  with  significant  charge 
sharing  in  the  bonds,  the  radii  of  the  atomic  spheres  is  frequently  increased  in 
Xa  calculations  so  that  an  overlap  region  appears  in  the  vicinity  of  the  bond. 

Our  studies  have  shown  that  the  contribution  to  the  total  molecular  energy 
from  the  exchange-correlation  term  shows  a  minimum  at  an  optimum  sphere  radius 
or  sphere  overlap.  This  result  can  be  deduced  from  the  Xa  model  by  comnaring  the 
exchange  energy  of  the  electron  density  at  the  edge  of  the  spherical  integration 
region  with  that  across  this  boundary  in  the  intersphere  region.  This  provides  a 


100 


sensitive  criterion  for  selecting  these  parameters. 

The  values  of  the  exchange  parameters  in  the  spherical  integration  region 

around  each  atomic  center  are  frequently  set  at  the  atomic  values  both  for 

neutral  and  for  ionic  molecular  constituents.  It  is  known,  however,  that  the 

21 

value  of  a  which  best  reproduces  Hartree  Fock  results  varies  with  ionicity. 

While  for  heavy  atoms,  these  changes  in  the  exchange  parameter  would  be  small, 

the  fit's  for  small  atoms  varj  rapidly  with  z  (and  with  ionicity).  The  correct 

choice  of  the  exchange  parameters  influences  not  only  the  total  energy  calculated 

for  the  molecule  but  also  affects  the  distribution  of  charge  between  the  atomic 

spheres  and  the  intersphere  region. 

Ab_Initio_  Aj>£r£a£h_-_Valenc_e_CI_  Method 

The  ab_  initio  calculation  procedure  chosen  for  these  studies  is  the  valence- 

23  24 

configuration-interaction  (VCl)  method.  ’  .The  specific  form  for  i//(r)  may  be 
written  as 


■Z'8'"®  >  =  E=X’U(E>  CD 

g 

where  eachi/^(R)  is  referred  to  as  a  configuration  and  has  the  general  structure 

^’“(R)  =  AOsn  I  [(Z'1<b)'Aj'2(b)  •  •  •  e;  j  (8) 

*  (B)-  •  •  •  •  C-1(A)K 

and  where  is  a  spatial  orbital,  /7 is  the  antisymmetrizing  onerator, O 

s 

is  the  spin-projection  operator  for  spin  quantum  number  S,  and  a 

product  of  ct  and  yS  one-electron  spin  functions  of  magnetic  quantum  number  M  . 

s 

The  spatial  orbitals,  0^  have  the  general  form 


0i  = 


-«i*  -*i»>  •  --'i' 


tP77q-rni“1pJ"11  (cos 6  ) 


(9) 


*» 
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,  the  standard  STO  form  is  recovered. 


where  for  £  =  ±6± 

The  wavefunctions  represented  by  Eq.  (8)  correctly  dissociate  to  an  atom  and  ion 
valence  bond  product  wave  function  at  large  internuclear  separations  with  an  op¬ 
timized  description  of  the  charge  density  about  the  seoarate  atom  and  ion  nuclei. 

This  situation  is  to  be  contrasted  with  that  arising  from  the  construction  of  a 
Hartree-Fock  wavefunction  for  the  molecular  ion.  Within  the  HF  framework  we  write 


4> 


g»u 


An[u 


S.u 


(a)u^  (0  )  .  .  .  u 


g.U 

2n-l 


(10) 


where 


S,u 


ui  =  0i(A)  t  4>i  =  av.  (  +  6^) 


(11) 


This  symmetry  adaptation  of  the  orbital  basis,  represented  bv  Eo .  (ll)  ,  intro¬ 
duces  a  severe,  non-physical  correlation  error  into  the  calculations  which  manifests 
itself  at  large  internuclear  spearations.  The  wavefunction  represented  by  Eq .  (10) 

dissociates  into  an  atom  and  ion  pair  with  charge  density  of  the  hypothetical 

+l/2 

species  R  .  This  situation  has  been  noted  by  previous  authors^,  who  suggested, 

however,  that  this  correlation  error  arising  from  symmetry  adaptation  did  not 

vary  significantly  with  changing  internuclear  separations.  We  find  different 

results,  especially  for  Neg  ,  which  separates  into  an  atom  and  ion  pair  with 

significantly  different  spatial  charge  distributions.  This  symmetry  dilemma 

arising  from  the  restricted  Hartree— Fock  model  has  been  discussed  previously  by 

T~  25 

Lowdin  who  suggested  alternate  computational  routes. 

There  are  several  complications  introduced  by  the  choice  of  Eq.  (8)  for  the 
molecular  ion  configuration.  For  heavy  nuclei  systems  such  as  Xe„+ ,  the  atomic 
and  ionic  basis  functions  will  be  similar  and  Eq.  (ll)  may  be  a  satisfactory 
approximation.  For  light  nuclei  systems  such  as  Neg"**,  the  distinction  between 
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the  atomic  and  ionic  basis  functions  must  be  retained  which  introduces  a  non¬ 
orthogonality  problem  into  Eq.  (8).  One  solution  is  to  remain  within  a  non- 
orthogonal  framework.  This  leads  to  a  secular  equation  of  the  form 


(4  -  =  0 


(12) 


where  J"  is  the  non-orthogonality  matrix.  Calculations  of  this  tyre  have  been 

2^  27 

reported  by  Harris  and  co-workers.  ’  An  alternate  approach  is  to  expand  the 
atomic  and  ionic  functions  as  linear  combinations  of  symmetry  adapted  MO's  such  as 


(Vuu); 


(13) 


This  leads  to  a  restricted  Cl  among  symmetry  adapted  functions  which  can  be  cast 
in  an  orthogonal  framework.  Both  approaches  were  examined  for  the  calculations 
reported  herein. 

Sp in -0 rbit  Effects 


The  ground  state  of  a  noble  gas  ion  is  split  by  spin-orbit  coupling  into  two 
2  2 

components,  a  state,  which  lies  lowest,  and  a  state.  This  splitting 

is  small  in  the  case  of  Ne  (?8l  cm  ^")  but  is  of  the  same  magnitude  as  the  chemical 
binding  energy  in  the  case  of  Xe  (10537  cm  ^) .  Rigorous  inclusion  of  spin- 

28 

orbit  effects  in  molecular  calculations  is  still  computationally  prohibitive^ 

but  for  these  noble  gas  ions,  which  exhibit  small  changes  from  separated  atom- 

7  Q 

ion  pair  electron  densities,  an  asymptotic  atomic  approximation  is  valid. 

In  this  approximation,  the  hamiltonian  is  written  as 


M R)  =  J.  A,s  (R)  +  vS"° 


(1M 
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where  ,  is  taken  in  the  atom— ion  limit  of  infinite  internuclear  separation  and 
is  considered  to  be  R- independent .  For  the  0  =  1/2  components,  Eq.  (lU)  leads  to 
a  secular  equation  of  the  form 


EZ(R)  -  -V2  a 

-Vja  lEn(E)  +  a]  - 


o  (15) 


z  n 

where  E  (r)  and  E  (r)  are  the  pure  A-S  coupled  eigenvalues  and  Ct  is  tne  spin- 
orbit  splitting  parameter  (3a  =  4E2P  2P  )  which  can  be  obtained  from  atomic 

3/2"  PV2 

data.  For  the  fl  =  3/2  components,  there  is  no  mixing  of  the  Z  and  n  states  and 
Eq.  (1U)  leads  simply  to 


n 


[E  (r)- a]  -  X 


3/2 


=  0 


(16) 


For  illustration,  the  spin-orbit  components  for  Q  =1/2  take  the  form 


W  1  -  C8.u<u  * 


(IT) 


'2  II 

g>u 


Cg,u^g,u  +  Cg,u^g,u 


where  -- gjU  -*  V2/3  and  C^u  -*  Vl/^  for  large  internuclear  separations.  For 
the  lighter  systems  such  asNe2+,  c|>u  -  1  and  c|jU  -  0  for  internuclear 
separations  near  Re .  Inclusion  of  these  spin-orbit  effects  was  carried  out  for  all 
of  the  ions  studied  herein.  We  find  that  this  effect  fully  accounts  for  the 
systematic  decrease  of  the  dissociation  energy  of  the  ground  state  of  these 
molecular  ions  in  the  progression  from  Ie2+  to  Xe2+. 
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RESULTS  AND  DISCUSSION 


We  have  carried  out  uniform  quality  density  functional  Xa  calculations  for 

the  molecular  ions  Ne„+ ,  Ar2+,  Kr2+  and  Xe2+  to  define  the  gound  A  11  £  state, 

2  2 

and  its  spectroscopic  properties,  and  the  nature  of  the  repulsive  B  n  ,  C  H 

g  u 

2  + 

and  D  £  states  which  arise  from  a  ground  state  noble  gas  atom-ion  pair .  As 
g 

a  check  on  the  accuracy  of  these  calculations,  ab_  inrtio_  wave  functions  were  con¬ 
structed  for  Ar„+  in  both  orthogonal  and  non-orthogonal  frameworks.  The  ab-initio 

calculations  were  carried  out  using  computational  programs  which  have  been  previously 
23  27 

described.  ’  The  density  functional  calculations  were  carried  out  with  both 

l6 

overlap  and  a  parameters  optimized  for  each  internuclear  separation.  There  is 
a  small  variation  in  the  optimum  parameters  arising  from  molecular  ion  symmetry 
considerations.  These  changes  were  found  to  be  negligible  for  Kr2+  and  Xe2+  but 
were  important  in  the  case  of  Ne2+.  Although  we  constructed  separate  optimum 
functions  for  each  molecular  ion  symmetry,  we  observe  that  transition  state  cal¬ 
cinations  of  energy  differences  between  these  ion  states  would  suffice  for  Kr2+ 

+  4* 

and  Xe2  but  would  definitely  be  in  error  for  Ne2  . 

The  general  nature  of  the  potential  energy  curves  arising  from  a  ground  state 
atom-ion  pair  is  illustrated  in  Fig.  1  for  Ar2+ •  These  curves  were  constructed  from 
a  Cl  expansion,  as  described  above,  using  a  minimum  STO  basis  augmented  by  the  in¬ 
clusion  of  optimized  4s  and  4p  functions  to  allow  more  flexibility  in  the  valence 


shell  charge  distribution.  Spin-orbit  effects  are  not  included  and  the  asymptotic 

2 

limit  is  the  center  of  gravity  of  the  P,  ,  argon  ion  components.  The  ground 

72  72 

2  + 

A  state  is  seen  to  be  bound  by  1.38  eV  with  an  equilibrium  separation  o^ 

O 

2.4a.  This  is  in  excellent  agreement  with  the  latest  experimental  study  of  this 

30  + 

ion  by  Mosely,  et  alJ  who  conclude  that  the  dissociation  energv  of  Ar2  is  1.33 

2  0 

±.02  eV.  There  is  a  small  binding  ( .1—. 2  eV)  of  the  B  II  state  at  R~3A  but  this 

g 
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state  is  clearly  repulsive  in  the  Franck -Condon  region  of  the  ground  state.  The 
2  2  -f 

higher  C  “nu  D  £  states  indicate  no  binding  at  the  internuclear  separations 
o 

^b.Qh)  that  we  examined,  although  a  weak  polarization  (-a/R  )  term  must  be  found 
at  larger  separations . 

Systematic  studies  of  these  ions  were  carried  out  using  the  density  functional 
Xa  method,  as  described  above.  Separate  calculations  were  performed  for  each 
molecular  ion  symmetry  and  for  a  wide  range  of  internuclear  separations.  A  double 
iteration  computation  scheme  was  employed  since  each  separate  calculation  was 
optimized  for  both  the  charge  distribution  and  the  radii  of  the  spheres  of  inte-* 
gration.  The  results  are  illustrated  in  Figs.  2-5  for  the  ion  sequence  Ne2+ , 

*4*  *+•  «4* 

Ar2  ,  Kr^  and  Xe2  .  In  these  figures  we  have  included  the  spin-orbit  splitting  of 
2 

the  n  A-s  coupled  states  and  the  resultant  mixing  of  the  1/2  _  components  with 
g,u  g,u 

2  2 

the  corresponding  2  and  £  states.  As  previously  mentioned,  this  splitting  is 

8  ^ 

+  4* 

small  for  Ne2  but  is  very  significant  for  Xe2  . 

Comparing  the  results  of  our  calculations  in  the  sequence  Ne2+-»  Xe2+  ,  we 

find  a  regular  progression  in  the  calculated  dissociation  energy  and  the  eauil- 

2  <4* 

ibrium  internuclear  separation  for  the  ground  1/21  [A  ]  state.  The  calcula- 

*4*  4* 

ted  binding  energy  decreases  from  1.3  eV  for  Ne2  to  1.0  eV  for  Xe^  while  the 

equilibrium  separation  increases  smoothly  from  3.2  to  5.9  bohrs .  A  summary  of 

our  calculated  spectroscopic  constants  for  these  dimer  ions  is  given  in  Table  I. 

These  data  were  obtained  by  fitting  our  calculated  J-T  coupled  potential  curves 

31 

to  the  Hulb ert -Hi rschf elder  modification  of  the  Morse  function  using  a  least- 
squares  fit  to  the  calculated  values  near  Re. 

An  examination  of  the  calculated  constants  reported  in  Table  I  indicates  some 
remarkable  regularities.  We  observe  that,  in  the  absence  of  spin-orbit  coupling 
effects,  the  dissociation  energy  of  all  four  dimer  ions  is  identical ,  within  the 
error  of  our  calculations.  This  suggests  that  the  binding  arising  from  A-S 
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coupling  in  these  systems  is  independent  of  the  principal  quantum  number  of  the 
valence  shell  and  that  the  experimentally  observed  decrease  in  the  binding 
energies  in  the  sequence  Ne2  -*  We2  has  its  origin  in  the  increased  degree  Oi 
mixing  of  the  repulsive  ^  (H)  component  into  the  ground  electronic  state.  A 
second  remarkable  simplicity  m  the  structure  of  these  ions  can  be  found  by  ex¬ 
amining  the  progression  in  the  calculated  equilibrium  separations,  ^n  Table  I 
we  list  a  column  of  the  average  expectation  of  <r>  for  these  ions,  taken  from 
accurate  Hartree-Fock  calculations,  and  compare  this  with  our  computed  values 
for  Pe.  We  find  the  condition  of  additivity  of  radii  is  an  excellent  approxima¬ 
tion  for  these  dimer  ions  with  R  (a.u.)  =  2.Q<?>+  constant. 

e 

Thus,  in  the  absence  of  spin-orbit  effects,  these  molecular  ions  exhibit  text¬ 


book  simplicity  for  our  understanding  of  their, binding  and  geometry. 

Further  confirmation  of  the  results  reported  here  can  be  obtained  by  com¬ 
parison  of  the  calculated  data  given  in  Table  I  with  experimental  determinations 
of  dissociation  energies  of  these  ions  from  photoionization  or  rainbow  scattering 


collision  studies.30 ’33~36  For  all  of  the  ions  studied  here,  differences  between 
theory  and  exneriment  lie  within  the  reported  experimental  uncertainties. 

in  addition  to  studies  of  the  stable  ground  state  of  these  ions,  we  show  in 
Figs.  2-5  the  calculated  J-J  coupled  excited  molecular  ion  states  which  correlate 
with  a  ground  state  atom-ion  pair.  There  are  seme  experimental  data  available 


to  verify  the  location  of  these  excited  states, from  rainbow  scattering  experiments35 
and,  more  recently,  from  photoabsorption  studies  in  Ar2  and  Kr2  .  We 

find  peak  absorption  for  the  strong  1/2  1  -*  1/2  II  transition  at  298  nm  for 
Ar2+  and  at  323  nm  for  Kr2+  .  These  calculated  data  can  be  compared  with  the  recent 
experimental  profiles  obtained  by  Hunter1  who  obtained  peaks  at  292  nm  and  325  nm 

for  Ar  +  and  Krp+  ,  respectively.  The  long  wavelength  transitions  for  photo- 

2  2  3  6 

absorption  (1/2  I  -  1/2  I  and  1/21  -  3/21)  have  been  studied  by  Miller,  et  al 
*  U  g  u  g 


who  report  an  overlapped  absorption  profile  with  a  peak  at  720  nm.  This  can  be 
compared  with  our  estimate  of  720-770  nm. 

Further  studies  of  the  photoabsorption  characteristics  of  these  noble  gas 
ions  are  in  progress . 
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Figure  1.  Ab  initio  potential  energy  curves  for 

without  inclusion  of  spin-orbit  coupling 


78-03 


Figure  2.  Density  functional  potential  energy  curves  for  Neg* 
(spin-orbit  effects  included) 
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Figure  3.  Density  functional  potential  energy  curves  for  Ar2 
(spin-orbit  effects  included) 
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Electronic  Structure  of  the  Noble  Gas  Dimer  Ions.  I. 
Potential  Energy  Curves  and  Spectroscopic  Constants* 


L.  A.  Wright 

Kirtland  Air  Force  Base,  New  Mexico  87117 

ABSTRACT 

A  systematic  study  of  the  electronic  structure  and  chemical  binding  in  the  dimer 

•f*  *4*  *4*  *4“ 

ion  sequence,  Ne2  ,  A ,  Kr2  and  Xe2  ,  has  been  carried  out  using  both  density 
functional  and  ab  initio  configuration-interaction  computational  techniques.  This  study 
includes  detailed  calculations  of  the  pertinent  potential  energy  curves  and  an  analysis 
of  the  calculated  spectroscopic  properties  of  the  bound  states  of  these  ions.  A  regular 
progression  is  found  in  the  spectroscopic  properties  for  the  ground  A  2^‘1/2u  state 
which  leads  to  some  remarkably  simple  conclusions  concerning  the  nature  of  the  binding 
and  the  size  of  these  dimer  ions.  For  the  heavier  systems,  Kr2+  and  Xe2+,  spin-orbit 
coupling  becomes  important ,  resulting  in  a  strong  mixture  of  the  A  -S  coupled  E  and  IT 
states.  This  mixing  affects  the  strength  of  the  binding  in  the  ground  state.  A  com¬ 
parison  with  other  ab  initio  studies  and  an  analysis  of  the  asymototic  behavior  at 
large  intemuclear  separations  is  given.  These  dimer  ion  species  illustrate  the 
classic  Hartree-Fock  symmetry  dilemma  arising  from  improper  dissociation  character. 

The  nature  of  this  problem  for  ionized  homopolar  species  is  discussed. 


Supported  in  part  by  AFWL  under  Contract  F29601-78-C-0031  and  by  AFOSR  under 
Contract  ?b962Q-jj-C-006h. 
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I.  INTRODUCTION 


The  diatomic  molecules  formed  from  the  noble  gases,  from  the  partially  closed 
shell  Group  IIA  and  IIB  elements  or  from  the  noble  gas-halogen  combinations  con¬ 
stitute  an  interesting  and  important  class  of  molecules  for  laser  applications. 

These  molecules  are  characterized  by  a  repulsive  or  weakly  bound  ground  state 
potential  curve  and  by  bound,  electronically  excited  states.  These  excited  diatomic 
molecules,  known  as  excimers,  radiate  in  a  narrow  visible  or  uv  continuum  band  by 
making  a  transition  to  a  repulsive  lower  state  that  dissociates  into  ground  state 
neutral  atoms.  The  resulting  laser's  operating  characteristics  are  controlled  by 
kinetic  processes  competing  for  the  formation  and  destruction  of  the  excited 
molecular  levels  and  by  absorption  losses  arising  from  noble  gas  dimer  ions  and 
high-lying  Rydberg  states  of  the  excimer  molecules.  This  last  observation  suggested 
the  need  for  the  present  studies  of  the  electronic  structure  and  absorption  character¬ 
istics  of  the  noble  gas  dimer  ions  which-  rate  studies  show  to  be  present  in  large 

1— U 

concentrations  in  the  candidate  excimer  systems . 

In  this  paper,  we  present  the  results  of  electronic  structure  calculations  for 

+  +  +  + 

the  dimer  ion  sequence  Ne2  ,  Arg  ,  Kr^  and  Xe^  .  These  results  include  the  per¬ 
tinent  potential  energy  curves  and  calculated  spectroscopic  constants.  Further 
papers  in  this  series  deal  with  photodissociation  and  absorption  processes  of 

the  noble  gas  dimer  ions.  These  include  the  predicted  absorption  spectrums  for 

o  +  2  +  2  +  2 

the  strong  A  -*■  D  Z  g  system  and  the  weaker  A  E^  B  Hg  system. 

Theoretical  treatments  of  the  noble  gas  dimer  ions  have  consisted  mainly  of 

2  + 

ab  initio  SCF  studies  of  the  ground  A  state  and,  by  far,  most  of  these 

studies  have  been  devoted  to  He2+-  We  consider  here  only  the  systems,  Ne2  ,  Ar^"*" , 

Kr2+  and  Xe2+,  which  are  chemically  similar,  with  binding  arising  from  the  outer 
shell  p-electrons. 
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The  Ne2+  and  Ar2+  ions  were  first  studied  by  Gilbert  and  Wahl5  using  SCP 
wave functions  constructed  from  STO  basis  functions.  These  authors  have  also  re- 

g 

ported  limited  SCF  studies  of  the  B  2H  c  2II  and  D  2£  +  excited  states  of  Ne0 

+  8  u  g  2 

and  Ar2  .  Cl  studies  of  the  ground  state  of  the  Ne2+  ion  were  reported  by  Cohen 

T  8 

and  Schneider  5  as  part  of  a  study  of  the  Rydberg  states  of  this  system.  SCP 
calculations  for  Ar2  at  short  internuclear  separations  have  been  reported  by 
Sidis,  et  al9  and  more  recently  by  Stevens,  et  al10,  using  an  extensive  STO 
oasis.  A  uniform  study  of  the  systems  Ar2+ ,  Kr2+  and  Xe2+  has  also  recently 
been  reported  by  Wadt11  who  employed  a  POL  Cl12  calculation  using  gaussian  basis 
functions.  This  study,  like  the  present,  was  directed  toward  the  prediction  of 
the  absorption  properties  of  these  dimer  ions.  Calculations  for  Xe2+  using  an 
effective  core  potential  have  been  reported  by  Ermler,  et  al.1^  These  cal¬ 
culations  were  part  of  a  more  extensive  study  of  Xe2  which  included  the  Rydberg 
states  of  this  system. 

Finally,  a  semi -empirical  analysis  of  the  potential  energy  curves  for  Xe2+ 
has  been  reported  by  Dehmer  and  Dehmer.14  These  authors  constructed  potential 
energy  curves  based,  in  part,  on  their  photoelectron  spectrum  of  the  Xe2  molecule 
on  the  elastic  scattering  measurements  of  difference  potentials  reported  by 

Jones,  et  al,  5  and  on  the  ab  initio  calculations  of  Wadt,11  which  were  used  in 
reduced  potential  form. 


+ 


5 


The  theoretical  principles  underlying  the  electronic  structure  of  the  homo- 

nuclear  noble  gas  molecular  ions  have  been  discussed  by  Mulliken1^  using  X&2+ 

as  a  prototype  and,  more  recently,  by  Gilmore,  Barr  and  Dee17  for  the  case  of 

Kr2  .  For  these  dimer  ions,  four  electronic  states  (a2£  +,  B2]T  ,  C2H  and  D2£  +) 

1  u  S’  u  g 

arise  from  a  ground  state  atom  (  S  )  and  a  ground  state  ion  (2P,  3J.  The  lowest 

'2' 

state,  usually  labeled  A  2Z^  ,  is  bound  for  the  four  systems  Ne2+,  Ar0+,  Kr2+ 
and  Xe2  and  has  been  partially  characterized  experimentally  through  scattering 
experiments,  although  there  has  been  much  controversy  concerning  the  dissociation 
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l6 

energies  of  these  ions.  For  the  heavier  ions,  A-S  coupling  is  no  longer  strictly 

valid  and  the  B  and  C  2n  states  split  into  J-J  coupled  conroonents  with  a 
8  u 

resultant  mixing  of  £  and  II  sionmetries  .  This  splitting  is  small  for  Ne^+  but  six 


distinct  J-J  coupled  molecular  ion  states  arise  for  Xe,. 


11 


With  the  exception  of  the  recent  calculations  by  Wadt,  no  attempts  at 
systematic  studies  of  these  ions  have  been  reported.  Accurate  ab  initio  methods 
would  appear  to  be  practical  only  for  the  lighter  ions  with  diminished  accuracy 

4.  11 

for  heavier  systems  such  as  Xe^  •  Wadt  reports  such  a  trend  in  his  calculations 
which  indicate  that  his  gaussian  basis  is  less  flexible  or  less  complete  in  the 

+  4* 

case  of  Xe2  than  for  the  corresponding  study  of  Arg  .  Recent  density  functional 

calculations  of  noble  gas  halide  molecules1^ have  yielded  reliable  predictions 

of  geometry  and  chemical  binding,  provided  certain  optimization  criteria  are 

followed.  Such  calculations,  which  should  be  of  more  uniform  quality  than  ab  initio 

studies,  since  numerical  molecular  orbitals  are  employed,  constitute  the  principle 

effort  in  the  present  study.  A  limited  number  of  ab  initio  calculations  were 

carried  out  for  Ar2+  for  comparison  purposes. 

II.  Method  of  Approach 

A.  Density  Functional  Approach  -  X^  Model 
20 

The  X^  model  for  the  electronic  structure  of  atoms,  molecules,  clusters 

and  solids  is  a  local  potential  model  obtained  by  making  a  simple  approximation 

to  the  exchange-correlation  energy.  If  we  assume  a  non-relativist ic  Hamiltonian 

21 

with  only  electrostatic  interactions,  it  can  be  shown  that  the  total  energy  E 
of  a  system  can  be  written  (in  atomic  units)  as 


E=  Z  Hj  <  U| 


--fvf  +  I  7 


fJL  'I/A 


U;  >  +  2  £  r 

1  \i--£V  rfJ.V 
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This  expresrion  is  exact  provided  the  u^  are  natural  orbitals  and  n^  are 
their  occupation  numbers  (i.e.,  eigenfunctions  and  eigenvalues  of  the  first  order 
density  matrix).  The  first  term  in  Eq.  (l)  represents  the  kinetic  and  electron- 
nuclear  energies.  The  second  term  is  the  nuclear  repulsion  energy.  The  sums 
(y,v)  are  over  all  the  nuclear  charges  in  the  system.  The  third  term  is  the  electron- 
electron  repulsion  term,  which  represents  the  classical  electrostatic  energy 
of  the  charge  density  p  interacting  with  itself,  where 


j(l)  =1  n j u*( •>  UjO) 

i 


(2) 


The  last  term  Exc  represents  the  exchange-correlation  energy  and  can  be  expressed 
formally  as 


Exc  =  ~f/^(l)  dri  I  d r2  , 

12 


(3) 


where  PXCU>2)  represents  the  exchange-correlation  hole  around  an  electron  at 
position  1.  In  the  exact  expression,  pxc  is  dependent  on  the  second-order  density 
matrix.  In  the  Hartree-Fock  approximation  Exc  is  simply  the  exchange  energy  and 
p  represents  the  Fermi  hole  due  to  the  exclusion  principle  and  depends  only  on 
the  first-order  density  matrix.  In  the  method,  a  simpler  assumption  about  pxc 

is  made.  If  we  assume  that  the  exchange-correlation  hole  is  centered  on  the  electron 
and  is  spherically  symmetric,  it  can  be  shown  that  the  exchange-correlation  potential 


Uxc  '  J 


dr„ 


12 


(4) 
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is  inversely  proportional  to  the  range  of  the  hole,  r  ,  where  r  is  defined,  by 

s  s 


3  r*  pU)  =  \  (5) 

Therefore,  in  the  X  model,  the  potential  U  is  proportional  to  p^^(r).  A 

xc 

scaling  parameter  a  is  defined  such  that 

uxa(l,=  "  ( 3/> (I )/87r )'/3  (g) 

The  expression  in  Eq.  6  is  defined  so  that  a= 2/3  for  the  case  of  a  free 

22 

electron  gas  in  the  Hartree-Fock  model  and  ct=l  for  the  potential  originally 

23 

suggested  by  Slater  .  A  convenient  way  to  choose  this  parameter  for  molecular 

and  solid  state  applications  is  to  optimize  the  solutions  to  the  X  equations  in 

2k 

the  atomic  limit.  Schwarz  has  done  this  for  atoms  from  Z  =  1  to  Z  =  Ul  and 
found  values  between  2/3  and  1.  However,  the  optimum  value  for  a  is  dependent 
on  the  ionicity  of  the  system.  For  accurate  total  energies,  this  variation  must 
be  properly  described  as  is  discussed  below. 

B .  Computational  Aspects  of  the  X^  Method 

To  implement  the  X^  approximation  in  molecular  calculations  it  is  usually 
convenient  to  break  up  the  volume  integrals  indicated  in  Eqs .  (3)  or  (1+)  by  de¬ 
fining  spherical  regions  surrounding  each  atomic  center.  The  molecular  potential 
is  then  spherically  averaged  within  these  regions  and  assumed  constant  in  the 
region  between  spheres.  The  potential  falls  to  zero  in  the  usual  way  in  the 
region  outside  of  an  overall  sphere  located  to  encompass  those  spheres  surrounding 
the  individual  atoms. 

There  are  two  other  practical  aspects  of  the  calculations  wh_ch  must  be 
considered  in  application  of  the  XQ  modi....  to  finite  molecular  systems.  The  first 
concerns  the  choice  of  the  integration  framework  for  describing  the  molecular- 
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orbital  wavef uncti ons  and  the  second  deals  with  the  choice  of  the  exchange  par¬ 
ameter,  a,  in  different  regions  of  space.  A  third  consideration,  convergence  of 
the  spherical  harmonic  expansion,  was  also  tested.  For  the  systems  studied  here, 
the  total  energy  was  converged  when  all  terms  up  to  £=3  were  included. 

In  molecules  with  significant  charge  sharing  in  the  bonds,  the  radii  of 
the  atomic  spheres  is  frequently  increased  in  calculations  so  that  an  overlap 
region  appears  in  the  vicinity  of  the  bond.2^  Our  studies19  have  shown  that 
the  contribution  to  the  total  molecular  energy  from  the  exchange-correlation 
term  shows  a  minimum  at  an  optimum  sphere  radius  or  sphere  overlap.  This  result 
can  be  deduced  from  the  ^  model  by  comparing  the  exchange  energy  of  the  electron 
density  at  the  edge  of  the  spherical  integration  region  with  that  across  this 
boundary  in  the  intersphere  region.  This  provides  a  sensitive  criterion  for 
selecting  these  parameters. 

The  values  of  the  exchange  parameters  in  the  spherical  integration  region 

around  each  atomic  center  are  frequently  set  at  the  atomic  values  both  for 

neutral  and  for  ionic  molecular  constituents.  It  is  known,  however,  that  the 

P  il 

value  of  Of  which  best  reproduces  Hartree  Fock  results  varies  with  ionicity. 

While  for  heavy  atoms,  these  changes  in  the  exchange  parameter  would  be  small, 
the  a's  for  small  atoms  vary  rapidly  with  z  (and  with  ionicity).  In  Table  I 
we  show  the  optimum  exchange  parameters  for  the  ion  pairs  included  in  this  study. 

It  is  obvious  that  the  exchange  parameters,  and  the  corresponding  electron  densities, 
must  be  optimized  for  correct  ion-atom  pair  separation.  The  correct  choice  of 
the  exchange  parameters  influences  not  only  the  total  energy  calculated  for 
the  molecule  but  also  affects  the  distribution  of  charge  between  the  atomic 
spheres  and  the  intersphere  region. 

C.  Ab  Initio  Approach  -  Valence  Cl  Method 

The  ab  initio  calculation  procedure  chosen  for  these  studies  is  the  valence- 

2  7  2  8 

configuration-interaction  (VCl)  method.  ’  The  specific  form  for  ^(R)  may  be 
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written  as 


0S’U  (R) 


-  L 


^  (R) 


(T) 


where  each  0g,U(R)  is  referred  to  as  a  configuration  and  has  the  general  structure 


<’U(R)  -  /?0s/7  [*„!<*  )0„2  (a)-  •  •  0Ja)|[0;i(b)0;2(b)  . 


•  0'  ,(B) 

^n-!  J  (Q) 


K: 


(B)  0^2  (B  ) 


0  (B) 

i/n 


<f>  '  (A)  <p  '  (A) 

r  '  n2 


0  (A) 

>n-i 


M? 


and  where  <(>  is  a  spatial  orbital,  /7  is  the  antisymmetrizing  operator, 

is  the  spin-projection  operator  for  spin  quantum  number  S,  and  is  a 

s 

product  of  a  and  B  one-electron  spin  functions  of  magnetic  quantum  number,  Ms> 

The  spatial  orbitals ,  $  ^  have  the  general  form 

(oosy  )  (9) 

1 

where  for  ?.  =  ±6 . ,  the  standard  STO  form  is  recovered, 
l  l  ’ 

The  wavefunctions  represented  by  Eq.  (8)  correctly  dissociate  to  an  atom 
and  ion  valence  bond  product  wavefunction  at  large  internuclear  separations  with 
an  optimized  description  of  the  charge  density  about  the  separate  atom  and  ion  nuclei. 

This  situation  is  to  be  contrasted  with  that  arising  from  the  construction  of  a 

•> 

Hartree-Fock  wavefunction  for  the  molecular  ion.  Within  the  HF  framework  we  write 


!//s’u  =  An  [  u^,u 


(cOu ®’u(0  ) 


Ug’u  (a' 
2n  -1 


(10) 
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where 


g,u 

Ui  =  (A) 


0i(B  )  ;  4>i  =  av.  (  4>i  +  0-  ) 


(11) 


This  symmetry  adaptation  of  the  orbital  basis,  represented  by  Eq.  (ll),  intro¬ 
duces  a  severe,  non-physical  correlation  error  into  the  calculations  which  manifests 
itself  at  large  internuclear  spearations.  The  wavefunction  represented  by  Eq.  (10) 
dissociates  into  an  averaged  atom-ion  pair  with.  charge  density  of  the  hypothetical 
species  R  1^2 .  This  situation  has  been  noted  by  previous  authors^,  who  suggested, 
however,  that  this  correlation  error  arising  from  symmetry  adaptation  did  not 
vary  significantly  with  changing  internuclear  separations.  We  find  different 
results,  especially  for  Ne2+,  which  separates  into  an  atom  and  ion  pair  with 
significantly  different  spatial  charge  distributions.  This  symmetry  dilemma 

arising  from  the  restricted  Hartree— Fock  model  has  been  discussed  previously  by 
"  .  29 

Lowdin  who  suggested  alternate  computational  routes. 

There  are  several  complications  introduced  by  the  choice  of  Eq.  (8)  for  the 
molecular  ion  configuration.  For  heavy  nuclei  systems  such  as  Xe2+,  the  atomic 
and  ionic  basis  functions  will  be  similar  and  Eq.  (ll)  may  be  a  satisfactory 
approximation.  For  light  nuclei  systems  such  as  Ne2  ,  the  distinction  between 
the  atomic  and  ionic  basis  functions  must  be  retained  which  introduces  a  non¬ 
orthogonality  problem  into  Eq.  (8).  One  solution  is  to  remain  within  a  non- 
orthogonal  framework.  This  leads  to  a  secular  equation  of  the  form 

(  i  ~  Al)^  =  0  (12) 
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where  £  is  the  non-orthogonality  matrix.  Calculations  of  this  type  have  been 
reported  by  Harris  and  co-workers."^  5 An  alternate  approach  is  to  expand  the 
atomic  and  ionic  functions  as  linear  combinations  of  symmetry  adapted  MO's  such 
as 


KiU»); 


K  1  <) 


(13) 


This  leads  to  a  restricted  Cl  among  symmetry  adapted  functions  which  can  be 
cast  in  an  orthogonal  framework. 

D.  Spin-Orbit  Effects 

The  ground  state  of  a  noble  gas  ion  is  split  by  spin-orbit  coupling  into  two 
2  2 

components,  a  state,  which  lies  lowest,  and  a  state-  This  splitting 

is  small  in  the  case  of  Ne  (T8l  cm  but  is  of  the  same  magnitude  as  the  chemical 

binding  energy  in  the  case  of  Xe  (10537  cm  Rigorous  incluaior  of  spin- 

orbit  effects  in  molecular  calculations  is  still  computationally  prohibitive,"^2 

but  for  these  noble  gas  ions,  which  exhibit  small  changes  from  separated  atom— 

7  Q 

ion  pair  electron  densities,  an  asymptotic  atomic  approximation  is  valid. 

In  this  approximation,  the  hamiltonian  is  written  as 


J{  R) 


(R)  +  V 


S-0 


lib) 


where  v  is  taken  in  the  atom-ion  limit  of  infinite  internuclear  separation  and 
is  considered  to  be  R- in dependent .  For  the  n  =  1/2  components,  Eq.  (lU)  leads  to 
a  secular  equation  of  the  form 


Er(R)  -  A 


3, 

2 


-V2  ct 


-  V2  a 


[  e/7(R)  +  a]  -  \ 


(15) 
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-here  EE(r)  and  e"(h)  are  the  pur.  A-S  coupled  eigenvalues  and  a  la  the  snin- 
orbit  a  putting  parameter  (3d  =  d  E,  _  2p  ,  ^  oan  be  obtained  ^ 

data.  For  the  0  -  3/2  components  ,  there  is  no  mixing  of  the  £  and  H  states 
and  Eq.  (lL)  leads  simply  to 


[E/7(R)-a]  -  A3/2 


0 


(16) 


For  illustration,  the  spin-orbit  components  for  a  =  1/2  take  the  form 


1  g,u  1 

<u 

<u 

,  cn 

Lg,U 

^g,u 

1  II 

g’U 

-‘?.u 

rE 

+  Cg,u 

c 

(17) 


where  C  gjU  -  WT  and  -  /l/3  for  large  intemuclear  separations.  For 

the  lighter  systems  such  as  He/,  *  1  and  cj>a*0  for  intemuclear 

separations  near  R_.  Inclusion  of  these  spin-orbit  effects  was  carried  out  for 
all  of  the  ions  studied  herein.  tie  shall  see  that  this  effect  fully  accounts  for 

the  systematic  decrease  of  the  dissociation  energy  of  the  ground  state  of  these 
molecular  ions  in  the  progression  from  Ne2+  to  Xe2+. 

•  Long-Range  Interactions 

Fne  long-range  energy  of  interaction  between  a  noble  gas  atom  and  a  singly 
charged  ion  (with  the  exception  of  He-He+)  can  be  written  as33’3** 

3q,oQ+2 

-1]/  "17“  (18) 


V(  R-h») 
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2R 


TT 


3aoa* 
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Jhe  first  and  third  terms  represent  the  main  components  of  energy  of  induction 
between  the  charge  and  the  induced  dipole.  The  second  term  is  the  London  dis¬ 
persion  energy  based  on  the  approximate  Drude  formula.  A  small  contribution 
arising  from  the  quadrupole  polarizability  of  the  noble  gas  atom  has  been 
neglected  in  this  term.  The  term  in  R~8  is  unimportant  for  these  dimer  ions 
owing  to  the  smallness  of  the  permanent  moments  for  these  systems . 

In  Table  II  we  give  values  for  the  parameters  appearing  in  Eq.  (l8).  Since 
dipole  polarizabilities  for  the  noble  gas  ions  are  unavailable,  they  are  estimated 
from  those  for  the  neutral  atoms  using  expectation  values,  <r2>,  calculated  from 
corresponding  neutral  and  ionic  Hartree-Fock  atomic  wavefuncticns .  This  approx¬ 
imation  appears  to  be  reliable  for  the  corresponding  halogen-halogen  ion  pairs, 
where  there  are  known  experimental  data  for  comparison.  The  long-range  form 
of  the  interaction  potentials  can  be  represented(in  a.u.)  as 


Ne-Ne+ :  U(R->-“>) 


1.831  _  2.7I+9 

~lr  R  6 


(19) 


Ar-Ar 


U(R+« 


5.5^0 


39.945 


(20) 


Kr-Kr+:  U(FH-°°) 


8.368 


88.697 

IT 


(21) 


Xe-Xe+ ;  U(R*°°) 


13.645 

IT 


219.327 

"7~ 


(22) 


Values  of  the  long-range  interaction  calculated  from  Eq.  (i9)-(22)  are  presented 
in  Table  III.  These  data  are  compared  below  with  the  results  obtained  from 
the  density  functional  calculations. 


III. 


RESULTS  AND  DISCUSSION 


We  have  carried  out  density  functional  calculations  for  the  molecular  ions 

Ne,-  +  ,  Ar„+,  Kr„+  and  Xe2+  to  define  the  ground  A  ^u+  state,  and  its  spectro- 

2  2  2  + 

scopic  nroperties,  and  the  nature  of  the  mainly  repulsive  B  I  ,  C  I  and  D  E 

g  u  g 

states  which  arise  from  a  ground  state  noble  gas  atom-ion  pair.  As  a  check  on  the 

accuracy  of  these  calculations  ,  ab  initio  wave  functions  were  constructed  for  Ar0+  in 

both  orthogonal  and  non-orthogonal  frameworks.  The  ab  initio  calculations  were 

27  71 

carried  out  using  computational  programs  which  have  been  previously  described.  ’ 

The  density  functional  calculations  were  carried  out  with  both  overlap  and  a 

19 

parameters  optimized  for  each  ion  symmetry.  There  is  a  small  variation  in  the 
optimum  parameters  arising  from  molecular  ion  symmetry  considerations.  These  changes 
were  found  to  be  negligible  for  Kr2+  and  Xe2+  but  were  important  in  the  case  of  Ne2+. 
Although  we  constructed  separate  optimum  functions  for  each  molecular  ion  symmetry, 
we  observed  that  transition  state  calculations  of  energy  differences  between  these 
ion  states  would  suffice  for  Kr2+  and  Xe2+  but  would  definitely  be  in  error  for  Ne2+. 
A.  Ab  Initio  Results 

The  general  nature  of  the  potential  energy  curves  arising  from  a  ground  state 
noble  gas  atom-ion  pair  can  be  illustrated  using  Ar2+  as  a  prototype.  Since  our  main 
emphasis  was  on  density  functional  methods,  which  appeared  to  offer  more  promise  of 
uniform  quality,  Ar2+  represented  a  typical  system  with  relative  ease  of  &b  initio 
computations . 

A  Cl  expansion  was  constructed  as  described  above  using  a  minimum  elliptic  basis 
optimized  at  each  internuclear  separation  studied,  and  augmented  by  the  inclusion 
of  optimized  Us  and  Up  functions  to  allow  more  flexibility  in  the  valence  shell 
charge  distribution.  In  the  computational  framework  described  by  Eq.  (8),  a  partial 
orthogonalization  of  the  basis  set  is  possible.  This  involves  orthogonalization 
of  the  core  functions  among  themselves,  followed  by  separate  orthogonalization 
of  the  valence  shell  functions  to  the  core  functions.  The  antisymmetrizing 
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operator  then  spans  only  the  valence  shell  electrons  and  the  inner  (Ne-core)  electrons 
are  taken  to  he  pairwise  singlet-coupled.  An  alternate  approach,  within  a  totally 
orthogonal  framework,  is  to  construct  a  double  set  of  symmetry  orbitals  based 
separately  on  the  optimized  neutral  atom  and  ion  basis  functions.  This  new  orbital 
basis  (essentially  equivalent  to  a  double-zeta  basis  for  the  valence  electrons) 
can  then  be  projected  into  a  restricted  Cl  equivalent  to  Eq.  (8),  by  expansion 
of  the  separate  atomic  end  ionic  orbitals  in  terms  of  their  symmetrized  representa¬ 
tion  as  given  in  Eq.  (ll).  Except  for  some  small  higher  excitation  contributions 
in  the  orthogonalized  Cl  framework,  the  calculated  results  (and  computational  effort), 
using  either  approach,  were  virtually  identical.  The  trade-off  is  between  a 
large  Cl  in  an  orthogonal  framework  and  more  complicated  matrix  elements  with 
the  non-orthogonal  basis. 

The  resultant  ab_  initio  potential  curves  for  Ar2+  are  given  in  Table  IV  and  illus¬ 
trated  in  Fig.  1.  Spin-orbit  effects  are  not  included  and  the  asymptotic  limit 

2  SS 

is  the  center  of  gravity  of  the  P3/2  1/2  argon  ion  components.  The  ground 
2  + 

A  state  is  seen  to  be  bound  by  1.26  eV  with  an  equilibrium  separation  of 

o 

2.42A.  This  is  in  reasonable  agreement  with  the  latest  experimental  study  of  this 

.  36  + 

ion  by  Moseley, et  al  who  conclude  that  the  dissociation  energy  of  Ar2  is  1.33 

±.02  eV.  Inclusion  of  spin-orbit  effects  would  tend  to  reduce  the  calculated 

dissociation  energy  of  Ar2+  and  somewhat  worsen  the  agreement.  There  is  a  small 

binding  (^  .1  eV)  of  the  B  'TT  state  at  Rv-3  A  but  this  state  is  clearly  repulsive 

S 

2 

in  the  Franck-Condon  region  of  the  ground  state.  The  C  11^  state  potential  curve 

o  2 

suggests  little  or  no  binding,  at  least  for  separations. <  3  A.  The  D  Z  state 

g 

is  repulsive  over  the  entire  range  of  internuclear  separations  that  were  studied. 

We  shall  see  that  very  similar  results  are  found  for  Ar2+  using  the  density 
functional  method. 
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B.  Dsnsitv  Functional  X  Results 

- d - 3 - 

Systematic  studies  of  these  dimer  noble  gas  ions  were  carried  out  using 

•  19  V,  a 

the  density  functional  X^  method,  as  described  above.  Previous  studies  naa 
indicated  that  these,  nearly  closed-shell,  interacting  atom-ion  systems  should  be 
favorable  candidates  for  application  of  density  functional  methods.  This  is 
based  on  the  observation  that  the  binding  in  these  systems  is  dominated  by 
simple  exchange  forces  and  that  no  drastic  charge  redistribution  occurs  during 
molecular  ion  formation.  Separate  calculations  were  performed  for  each  molecular 
ion  symmetry  and  for  a  wide  range  of  internuclear  separations.  A  double  iteration 
computation  scheme  was  employed  since  each  separate  calculation  was  optimized 
for  both  the  charge  distribution  and  the  radii  of  the  spheres  of  integration. 

The  optimum  overlap  radii  vary  slightly  with  the  symmetry  of  the  species  but 
are  close  to  the  values,  2k. 5,  21.0,  20.5  and  20.0$,  respectively,  for  the 
sequence  Ne2+  -*  Xe2+.  This  optimization  is  very  important  if  accurate  total, 
energies  are  required.  The  one-electron  excitation  spectra  are  much  less  sen¬ 
sitive  to  overlap  optimization. 

The  resultant  density  functional  potential  energy  curves,  assuming  A-S  coupling, 
are  given  in  Tables  V-VIII  for  the  ion  sequence  Ne2+,  Ar2+,  Kr2+  and  Xe2+.  A  com¬ 
parison  can  be  made  between  these  calculated  potential  curves  and  the  long-range  ion- 
atom  potentials,  as  given  in  Table  III,  at  large  internuclear  separations.  Con¬ 
sidering  only  g-u  symmetry  effects,  the  B  II  and  C  states  have  asymptotic 

o  ^ 

behavior  of  the  form 

Ug(u)(R"”)  =  u  (R  ^  ®)  (7)  <  XA  XB  l^XB  XA  >  (23) 

where  U(R*»)  is  given  by  Eqs.  19-22  and  and  xj  refer  to  separated  atom  and  ion 
wavefunctions  located  on  centers  A  and  B,  respectively.  A  similar  argument  can  be 
made  for  states  of  I  symmetry  which  is  valid  at  somewhat  larger  internuclear 


131 


separations.  We  see  that  Eq.  (23)  closely  corresponds  to  our  density  functional 

results;  the  g-u  states  split  the  long-range  polarization  potential.  The  ab_  initio 

potential  energy  curves  reported  by  Wadt11  do  not  exhibit  this  behavior.  In  his 

studies,  the  B  "n  state  for  Ar 2+  -*•  Xe2+  closely  follows  Eqs.  20-22  indicating  too 

much  long-range  repulsion  in  the  calculated  potential  energy  curves . 

In  Tables  IX-XII  we  have  included  the  spin-orbit  splitting  of  the  A-S  coupled 

states  and  the  resultant  mixing  of  the  0  —  1/2  component  of  the  corresponding 
2  2 

I  and  H  states.  As  previously  mentioned,  this  splitting  is  small  for  Ne2+ 

but  is  very  significant  for  Xe2  .  Figures  2-5  illustrate  our  calculated  potential 

energy  curves  including  these  spin-orbit  coupling  effects. 

Summaries  of  our  calculated  spectroscopic  constants  are  given  in  Tables 

XIII-XVI  for  the  dimer  ion  sequence  Ne2  ->-Xe2  .  These  spectroscopic  data  were 

obtained  by  fitting  our  calculated  J-J  coupled  potential  curves  to  the  Hulbert- 

37 

Hirschfelder  modification  of  the  Morse  function  using  a  least-squares  fit  to 
the  calculated  values  near  Re>  Comparing  the  results  of  our  calculations  in  the 

4*  4* 

sequence  Ne2  ■*  Xe2  we  find  a  regular  progression  in  the  cal  ulated  dissociation 
energy  and  the  equilibrium  internuclear  separation  for  the  ground  A  l  ^2u  state. 

The  calculated  binding  energy  decreases  from  1.31  eV  for  Ne2+  to  1.06  eV  for 

4*  O 

Xe2  while  the  equilibrium  separation  increases  uniformly  from  1.7  to  3.2  A. 

A  comparison  of  our  calculated  data  for  the  A  l  i/2m  state  of  these  ions  is 
given  in  Table  XVII.  An  examination  of  these  calculated  data  reveals  some 
remarkable  regularities.  We  observe  that,  in  the  absence  of  spin-orbit  coupling 
effects,  the  dissociation  energy  of  all  four  dimer  ions  is  identical,  within 
the  error  of  our  calculations  (1.38  ±  .02  eV).  This  suggests  that  the  binding  arising  from 
A-S  coupling  in  these  systems  is  independent  of  the  principal  quantum  number  of  the 
valence  shell  and  that  the  experimentally  observed  decrease  in  the  binding  energies 
in  the  sequence  Ne2  -«  Xe2  has  its  origin  in  the  increased  degree  of  mixing  of 
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the  repulsive  It-^u  comP°nen1:  into  the  ground  electronic  state.  A  second  re¬ 
markable  simplicity  in  the  structure  of  these  ions  can  be  found  by  examining  the 
progression  in  the  calculated  equilibrium  separations.  In  Table  XVII  we  list  a 

column  of  the  average  expectation  of  <r>  for  these  atom-ion  pairs,  taken  from 

26  38 

accurate  Hartree-Fock  calculations,  ’  and  compare  this  with  our  computed  values 

for  .  We  find  the  condition  of  additivity  of  radii  is  an  excellent  approximation 

for  these  dimer  ions  with  R  (a.u.)  =  2.0  <r>  +  constant.  Thus,  in  the  absence 

of  spin-orbit  effects,  these  molecular  ions  exhibit  text-book  simplicity  for  our 

understanding  of  their  binding  and  geometry. 

A  comparison  of  the  density  functional  potential  energy  curves  given  in 

Tables  V-VIII  can  be  made  with  the  recent  ab  initio  calculations  reported  by  Wadt11 

for  the  ions  Ar2  ,  Kr2  and  Xe^  and  the  ab  init io  with  effective  core  notential 

13 

calculations  recently  reported  by  Ermler,  et .  al  for  Xe2  .  The  ab  initio  potential 

curves  of  Wadt11  all  tend  to  exhibit  more  repulsive  energy  than  the  corresponding 

density  functional  curves.  These  differences  increase  as  the  intemuclear  separation 

decreases,  with  the  result  that  the  ab  initio  potential  curves  for  the  ground 
2 

A  £]_/2u  stages  underestimate  the  true  dissociation  energies.  The  differences  are 

correspondingly  greater  for  Xe2+  than  for  Ar 2+  which  suggests  a  deficiency  in  the 

analytic  basis  functions  chosen  for  the  ab  initio  studies.  The  studies  by 
13 

Ermler,  et .  al  point  out  the  significant  savings  in  the  use  of  core  potentials 
but  suffer  from  a  common  deficiency  in  the  basis  functions  chosen  for  the  valence 
shell  electrons. 

At  short  internuclear  separations  (<Re'1,  the  differences  between  the  ab  initio 
curves  and  our  density  functional  curves  appear  to  diminish.  This  suggests  an 
accumulating  error  in  the  density  functional  results  at  shorter  separations,  nrobably 
arising  from  core  overlap  effects  which  are  not  fully  represented  in  our  calculations. 
Qualitatively,  the  ab  initio  and  density  functional  results  are  remarkably  similar. 
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A  further  comparison  of  our  calculated  dissociation  energies  can  be  made  with 

lit  It 5 

the  recent  data  reported  by  Dehmer  and  Dehmer  *  for  the  excited  states  of 

Ar2+,  Kr0+  and  Xe0+  and  previous  experimental  studies  of  the  ground  state.  This 

comparison  is  shown  in  Table  XVIII  which  indicates  good  agreement  for  all  of 

lit 

these  dimer  ions.  Dehmer  and  Dehmer  also  report  semi-empirical  potential  energy 
curves  for  Xe0+  based,  in  part,  on  their  photoelectron  spectrum  of  the  Xe„  molecule, 
on  the  elastic  scattering  measurements  of  difference  potentials  reported  by 
Jones,  et.  ■al’'-'*,  and  on  the  ab  initio  calculations  of  Wadt,‘'“L  which  were  used  in 
reduced  potential  form.  These  semi -empirical  curves  all  qualitatively  disagree  both 
with  our  results  for  Xe,_,+  and  with  the  ab  initio  results  reported  by  Wadt.  The 
curves  reported  by  Dehmer  and  Dehmer  exhibit  much  flatter  shapes  with  softer 
repulsive  branches  for  all  symmetries  of  Xe2+ .  Part  of  the  difference  arises 
from  their  assumption  of  a  constant  scaling  error  in  using  Wadt's  potential 
curves  in  reduced  (U(R)/De)  form.  Ab  initio  potential  curves,  in  general,  cannot 
be  characterized  by  an  error  term  of  their  functional  form. 

Partial  confirmation  of  the  results  reported  here  can  also  be  obtained  by  com¬ 
parison  of  the  calculated  data  given  in  Tables  XIII-XVI  with  experimental  determina¬ 
tions  of  dissociation  energies  of  these  ions  from  photoionization  or  rainbow 

36  39  — i+3 

scattering  collision  studies.  ’  For  all  of  the  ions  studied  here,  differences 

between  theory  and  the  most  recent  experimental  values  for  the  dissociation  energies 

lie  within  the  reported  experimental  uncertainties.  There  are  also  some  experimental 

data  available  to  verify  the  location  of  the  excited  states  of  these  dimer  ions. 

l+l 

These  data  are  from  rainbow  scattering  experiments  and,  more  recently,  from 

+  +  1  36  42-44 

photoabsorption  studies  in  Ar 2  and  Kr2  .  ’  ’  We  find  a  peak  calculated 

2  2 

absorption  for  the  strong  A  D  Sl/2g  transition  at  298  nm  for 

Ar~,  and  at  323  nm  for  Kr2+.  These  calculated  data  can  be  compared  with  the  recent 
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experimental  profiles  obtained  by  Hunter1  who  obtained  peaks  at  292  run  and  325  run 

■4*  + 

for  kx 2  and  Kr2  ,  respectively.  The  long  wavelength  transitions  for  photo- 
2  2  2  2 

absorption  (A  ^p/2u->  ®  ni/2g  811(1  ^  ^  ^/2g^  have  been  studied  in  detail  L'y 

36  42-43  +  +  + 

Moseley  et  al,  ’  for  Ar2  ,  Kr2  and  Xe2  .  Further  papers  in  this  series  will 

report  more  detailed  comparisons  of  these  photoabsorption  spectra  with  our  calculated 
profiles . 

The  apparent  success  of  the  density  functional  method  in  describing  the  electronic 
structure  of  these  noble  gas  dimer  ions  should  not  be  interpreted  as  suggesting  that 
this  method,  in  general,  is  to  be  preferred  over  ab  initio  methods.  As  has  been 
pointed  out  above,  closed  shell  system  interactions,  where  there  are  no  drastic  charge 
rearrangemerts  and  no  degeneracy  complications,  are  ideal  candidates  for  the  density 
functional  method.  This  method,  when  properly  applied,  is  capable  of  yielding  results 
of  uniform  accuracy  for  both  light  and  heavy  atom  systems  owing  to  numerical  representa¬ 
tion  of  the  valence  shell  orbitals.  The  well-known  analytic  expansion  problems  of 
ab  initio  basis  functions  are,  in  a  sense,  circumvented.  In  addition,  there  is  a  large 
body  of  accumulating  evidence  that  density  functional  methods,  in  general,  can  yield 
reliable  one-electron  excitation  spectra. 

Hcwever ,  the  density  functional  method,  at  least  in  any  of  its  present  simplified 
forms,  is  not  particularly  useful  for  total  energy  calculations  of  systems  which  exhibit 
a  large  charge  rearrangement  with  a  change  in  intemuclear  separation.  In  addition, 
the  method  is  not  useful  at  all  for  systems  with  electronic  degeneracies  or  com¬ 
plicated  spin  representation.  In  these  regards,  it  suffers  the  same  limitations 
as  conventional  wartree-Fock  methods.  When  such  considerations  are  present,  and 
unfortunately  they  are  in  the  general  case  and  not  as  exceptions,  ab  initio  techniques 
must  be  relied  upon  for  quantitative  results. 

Our  conclusion  is,  when  properly  used  and  applied  to  systems  involving  simple 
charge  descriptions,  density  functional  methods  are  a  valuable  addition  to  our 
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list  Of  techniques  for  electronic  structure  calculations. 
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Table  I.  Optimum  Atomic  and  Ioni 


c  Density  Functional  Exchange  Parameters 


and  Energies . 


-fecies  Exchange  Parameter, a 


We 

0.73081 

..  + 

We 

0.73684 

Ar 

0.72177 

Ar+ 

0.72249 

Kr 

0.70574 

+ 

Kr 

0.70583 

Xe 

0.69986 

__  + 

Xe 

0.69986 

Xn  energy  (a.u.) 

£ 

HF  energy 

-128.5470 

-128.5474 

-127.8176 

-127.8178 

-526.8176 

-526.8173 

-526.2745 

-526.2745 

-2752.0560 

-2752.0569 

-2751.5679 

-2751.5674 

-7232.1501 

-7232.1500 

-7231.7114 

7231.7118 

a.u. ) 


Table  II.  Long-Range  Force  Interaction  Parameters  for  the  Noble  Gas  Dimer  Ions. 


Species 

aQ  (bohr^) 

£  <r“T  >  (bohr“ ) 

1  1 

I  (eV) 

Ne 

2.66  a’b 

9.699  d 

21.561+ 

Ne+ 

(1.32)  c 

6.81+2 

1+0.962 

Ar 

11.08 

26.11+5 

15.759 

Ar+ 

(6.52) 

20.051+ 

27.629 

Kr 

16.7^ 

39.671+ 

13.999 

Kr+ 

(10.78) 

31.81+7 

21+.359 

Xe 

27.29 

62.511 

12.127 

Xe+ 

(18.90) 

52.022 

21.2 

a 


b 

c 


d 

e 

f 
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The  ion  polarizabil.  'es  are  scaled  from  the  atomic  o’  using  the 

o 

Non-relativistic  Hartree-Fock.  calculation. 

R.  L.  Kelly  and  D.  E.  Harrison,  Jr.,  Atomic  Data  3j_  177  (.1971).. 
C.  E.  Moore,  Nat.  Bur.  Std.  Cu.S.)  Circ.  I+67  (.1958). 


hydrogenic  form. 
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Table  III. 

Long  Range  Interaction  Potentials 

for  the 

Noble  Gas 

Dimer  Ions 

R  (a.u. ) 

U(R)  (a.u.) 

R  (a.u. ) 

U(R)  (a 

N  e-N  e+ 

Kr-Kr+ 

8.0 

-.00034 

10.0 

-.00093 

7.0 

-.00058 

9.0 

-.00144 

6.0 

-.00109 

8.0 

-.00238 

5.0 

-.00326 

7.5 

-.00314 

4 .  6 

-.00393 

7.2 

-.00375 

4.4 

-.00478 

7.0 

-.00424 

4.2 

-.00587 

6.8 

-.00481 

6.4 

-.00628 

Ar-Ar+ 

Xe-Xe+ 

9.0 

-.00092 

11.0 

-.00106 

8.0 

-.00151 

10.0 

-.00158 

7.0 

-.00265 

9.0 

-.00249 

6.5 

-.00363 

8.2 

-.00374 

6.2 

-.00445 

8.0 

-.00417 

6.0 

-.00513 

7.6 

-.00523 

5.8 

-.00595 

7.4 

-.00589 

5.4 

-.00813 

7.0 

-00755 
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Table  IV.  Ab  Initio  Potential  Energy  Curves  for  Ar0+ 
assuming  A-3  coupling.  Energies  in  a.  u. 


R  ( a.u. ) 

A 

b  h 

—  g 

— nu 

D  2l  ' 

-  g 

3.5 

.0869 

.2727 

.3923 

.4707 

3.75 

.0145 

- 

- 

.3440 

4.0 

-.0238 

- 

- 

.2388 

4.25 

-.0409 

- 

- 

.1674 

4.5 

-.0459 

.0264 

.0593 

.1185 

4.75 

-.0454 

- 

- 

.0847 

5.0 

-.0397 

- 

- 

.0596 

5.5 

-.0278 

-.0013 

.0069 

.0299 

6.0 

-.0177 

-.0032 

-.0006 

•  0143 

00 

0. 

0. 

0. 

0. 
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Density  Functional  Potential  Energy  Curves  for  Ne£+  assuming 
iV"S  coupling.  Energies  in  a.u.  relative  to  E(»)=  -256.3646  a.u. 


R  (a.u. ) 

A2E  + 

- u 

B  2TT 

-  8 

— nu 

d2d  + 

—  g 

3.0 

-.04737 

.04756 

•09573 

.16557 

3.2 

-.05068 

.01916 

.05401 

.11261 

3.4 

-.04840 

.00404 

.02933 

.07765 

3.6 

-.04340 

-.00355 

.01487 

.05417 

3.8 

-.03729 

-.  00686 

.00659 

.03820 

4.0 

-.03101 

-.00785 

.00200 

.02720 

4.2 

-.02515 

-.00767 

-.00044 

.01951 

4.4 

-.01998 

-.00693 

-.00161 

.01410 

4.6 

-.01563 

-.00599 

-.00207  ^ 

.01024 

CO 

0. 

0. 

0. 

0. 

Table  VI. 


+ 


Density  Functional  Potential  Energy  Curves  for  Ar 
\-S  coupling.  Energies  in  a.u.  relative  to  E  (°°)2= 


assuming 

-1053.0921. 


R  (a.u. ) 

~+u 

b  2n 

—  g 

—  u 

D  2E  + 

-  g 

4.2 

-.04353 

.06295 

.10399 

.18203 

4.6 

-.05070 

. 01862 

.04412 

.10336 

4.8 

-.04945 

.00756 

.02765 

.07887 

5.0 

-.04659 

.00081 

.01663 

.06072 

5.2 

-.04279 

-.00316 

.00931 

.04707 

5.!+ 

-.03855 

-.00504 

.00478 

.03706 

5.8 

-.02991 

-.00608 

.00001 

.02330 

6.0 

-.02589 

-.00566 

-.00086 

.01890 

6.2 

-.02221 

-.00515 

-.00137 

.01532 

Table  VII.  Density  Functional  Potential  Energy  Curves  for  Krn+ 
assuming  \-S  coupling.  Energies  in  a.u.  relative" 
to  E(»)  =  -5503.6239 


P.  ( a.u.  ) 

Aul 

b  2n 

- g 

0 

c  ~n^ 

D  2Z  ' 

- g 

5.0 

-.04987 

.03373 

.06103 

.11527 

5.2 

-.05063 

.01831 

.04035 

.08833 

5.k 

-.04949 

.00843 

.02613 

.06817 

5.6 

-.04701 

.00210 

.01635 

.05288 

5.8 

-.04373 

-.00243 

.00920 

.04059 

6.0 

-.03999 

-.00420 

.00505 

.03206 

6.4 

-.03211 

-.00620 

-.00021 

.01933 

6.8 

-.02468 

-.00616 

-.00229 

.01159 

7.0 

-.02136 

-.00577 

-.00266 

.00897 

7.2 

-.01835 

-.00532 

-.00282 

. 00690 

00 

0. 

0. 

0. 

0. 
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Table  VIII.  Density  Functional  Potential  Energy  Curves  for  Xe* 
assuming  A.  — S  coupling.  Energies  in  a.u.  relative"- 
to  E(»)  =  -II+I+63.8615 


R(a.u. ) 

-  u— 

b 

—  g 

c  2n^ 

D  2Z  + 

- S 

5.8 

-.05028 

.02817 

.05042 

.09204 

6.0 

-.05061+ 

.01683 

.03519 

.07264 

6.2 

-.01+955 

.00992 

.02505 

.05813 

6.1+ 

-.01+748 

.00359 

.01606 

.04531 

6 . 6 

-.01+1+75 

-.00037 

.00989 

.03544 

6.8 

-.01+161 

-.00311 

.00533 

.02750 

7.0 

-.03825 

-.00455 

. 00239 

.02142 

7.1+ 

-.03139 

-.00548 

-.00080 

.01292 

7.6 

-.02808 

-.00537 

-.00153 

.01000 

7.8 

-.021+94 

-.00510 

-.00166 

.00768 

8.0 

-.02202 

-. 00469 

-.00211 

.00580 

8.2 

-.01932 

-.004i8 

-.00205 

.00451 

00 

0. 

0. 

0. 

0. 
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-able  IX.  Potential  Energy  Curves  for  Ne2+  including  spin-orbit  counling. 
Energies  in  eV  relative  to  c.g.  of  Ne0+  at  R=<». 


R  (a.u. ) 

A  2  + 
i-£-l/2u 

“  ^3/2g 

0  2tt 

— IL3/2u 

S-H-1/2U 

2  + 

2-S-i/. 

2.8 

-0.89497 

- 

- 

- 

- 

- 

3.0 

-1.28953 

1.26194 

1.32586 

2.57256 

2.63768 

4.50598 

3.2 

-1.37990 

0.48908 

0.55283 

1.43735 

1.50265 

3.06507 

3.4 

-1.31808 

0.07755 

0.14109 

0.76574 

0.83129 

2.11391 

3.6 

-1.18232 

-0.12882 

-0.06558 

0.37242 

0.43830 

1.47550 

3.8 

-1.01630 

-0.21895 

-0.15609 

0.14705 

0.21334 

1.04125 

4.0 

-0.84617 

-0.24597 

-O.18362 

0.02207 

0.08890 

0.74236 

4.2 

-0.68726 

-0.24091 

-0.17924 

-0. 044i6 

0.02338 

0.53395 

4.4 

-0.54766 

-0.22074 

-0.15998 

-0.07598 

-0.00748 

0.38739 

4.6 

-0.43046 

-0.19535 

-0.13579 

-0.08868 

-0.01897 

0.28363 

00 

-0.03230 

-0.03230 

-0.03230 

-0.03230 

0.06460 

0.06460 
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Table  X. 


Potential  Energy  Curves  for  Ar 2  including  spin-orbit  coupling. 
Energies  in  eV  relative  to  c.g.  of  Ar2+  at  R=®. 


R  ( a.u. ) 

A  z+ 

— -  1/2 

—  3/2g 

— ni/2g 

— I3/2u 

C  2"  l/2u 

D  V,  , 

4.2 

-I.l86l8 

1.65365 

1.76982 

2.77055 

2.89064 

4.95560 

4.4 

-1.34226 

- 

- 

- 

.. 

4.6 

-1.38227 

O.il+736 

0.56260 

1.14139 

1.26240 

2.81570 

4.8 

-1.34896 

0. 14654 

0.26118 

0.69322 

0.81483 

2.14999 

5.0 

-1.27160 

-0.03714 

0.07676 

0.39334 

0.51563 

I.65681 

5.2 

-1.16918 

-0.14504 

-0.03200 

0.19416 

0.31726 

1.28605 

5.4 

-1.05453 

-0.19620 

-0.08425 

0.07089 

0.19489 

1.01491 

5.8 

-0.82182 

-0.22463 

-0.11562 

-0.05891 

0.06740 

0.64328 

6.0 

-0.71382 

-0.21307 

-0.10599 

-0.08259 

0.04512 

0.52554 

6.2 

-0.61523 

-0.19933 

-0.09464 

-0.09647 

0.03290 

0.43043 

00 

-0.05919 

-0.05919 

-0.05919 

-0.05919 

0.11837 

0.11837 

Table  XI. 


Potential  energy  curves  for  Kr2+  including  spin-orbit  coupling. 
Energies  in  eV  relative  to  c.gT  of  Kr„+  at  R=»  . 


R  ( a.u. ) 

4.8 

A 

— -  l/2u 

-1.28087 

— n3/2g 

— 1 Tl/2g 

C-~h/2u 

^-ni/2u 

2-^  l/2u 

5.0 

-1.28709 

0.69575 

1.09151 

1.43876 

1.91284 

3.18488 

5.2 

-1.41503 

0.27623 

0. 66358 

0.87603 

1.35598 

2.46023 

5.4 

-1.38910 

0.00729 

0.38423 

0.48908 

0.97546 

1.92201 

5.6 

-1.32866 

-0.16482 

0.19960 

0.22295 

0.71627 

1.51846 

5.8 

-1.24730 

-0.28795 

0.06150 

0.02825 

0.52948 

1.19889 

6.0 

-1.15345 

-0.33612 

-0.00457 

-0.08468 

0.42439 

0.98465 

6.4 

-0.95763 

-0.39054 

-0.  j.0324 

-0.22768 

0.30020 

0.68250 

6.8 

-0.77681 

-0.38945 

-0.15500 

-0.28414 

0.26501 

0.52489 

7.0 

-0.69752 

-0.37898 

-0.17192 

-0.27435 

0.26591 

0.48097 

7.2 

-0.62696 

-0.36673 

-0.18626 

-0.29870 

0.27285 

0.45125 

CD 

-0.22197 

-0.22197 

-0.22197 

-0.22197 

0.44393 

0.44393 

147 


Table  XII. 

Potential  Energy  Curves  for  Xe 

Energies  in  eV  relative  to  c.g 

+  including 
.  of  at 

spin -orbit 

R=“ . 

coupling. 

R  (a.u.) 

2  + 

2 

2 

„  2 

2 

2_+ 

A  2  l/2u 

2-n3/2g 

“ — ^l/2g 

— — ^3/2u 

— 1l/2u 

— •  1/i 

5.6 

-1.40193 

- 

- 

- 

- 

- 

5.8 

-1.48344 

0.33100 

0.95693 

0.93646 

1.92269 

2.74968 

6.0 

-1.50867 

0.02237 

0.61482 

0.52197 

1.52369 

2.25511 

6.2 

-1.49362 

-0.16564 

0.38771 

0.24616 

1.26243 

1.89954 

6.4 

-1.45490 

-0.33783 

0.17473 

0.00150 

1.03543 

1.59142 

6.6 

-1.39811 

-0.44565 

0.02261 

-0.16640 

0.88498 

1.36717 

6.8 

-1.33076 

-0.52015 

-0.09755 

-0.29043 

0.77902 

1.19677 

7.0 

-1.25661 

-0.55939 

-0.18338 

-0.37035 

0.71648 

1.07791 

7.4 

-1.10398 

-0.58459 

-0.29772 

-0.45724 

0.66367 

0.93575 

7.6 

-1.03040 

-0.58149 

-0.33514 

-0.47699 

0.66026 

0.89673 

7.8 

-0.95973 

-0.57425 

-0.36461 

-0.48069 

0.67123 

0.87026 

8.0 

-O.89665 

-0.56304 

-0.38755 

-0.49278 

0.67575 

0.85330 

8.2 

-0.83747 

-0.54924 

-0.40162 

-0.49133 

0.69129 

0.84609 

ao 

-0.43549 

-0.43549 

-0.43549 

-0.43549 

0.87098 

0.87098 
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Table  XIII.  Spectroscopic  Constants  for  Ne0+ 


State 

Te  (eV) 

a)  ( cm  ^ ) 

to  X  (cm-1 
e  e 

)  ag( cm 

reU) 

Bg(cm  ''') 

D  (eV) 

p 

D  (eV) 

D 

V 

l/2g 

4.445 

(vertical 

excitation 

energy,  repulsive  curve 

) 

0 

C 

ni/2u 

1.359 

122.5 

8.4 

o . 0166 

2.50 

0.2665 

0.085 

0.078 

C 

n3/2u 

1.287 

103.9 

8.7 

0.0194 

2.54 

0.2579 

0.060 

0.054 

B 

\/2g 

P 

1.195 

246.5 

9.0 

0.0125 

2.15 

0.3626 

0.152 

0.137 

l 

n3/2g 

O  4- 

1.132 

250.8 

9.8 

0.0122 

2.14 

0.3634 

0.215 

0.199 

1 

r 

*  l/2u 

0.0 

597.0 

6.1 

0.0090 

1.69 

0.5840 

1.347 

1.310 
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Table  XIV. 


Spectroscopic  Constants 


for  Arg 


+ 


State 

T  (eV) 
e 

aje(cm_1) 

u)  X  (cm"1) 
e  e 

a  (cm"1) 
e 

re(A) 

B  ( cm"1 ) 
e 

D  (e\ 
e 

2  + 

1  l/2g 

1+.198 

(vertical  excitation  energy. 

repulsive  curve) 

\/2u 

1.1+13 

62.1+ 

9.3 

0.001+7 

3.1+0 

0.0732 

0.087 

2 

n3/2u 

1.285 

53.7 

3.7 

0.001+3 

3.1+0 

0.0732 

0.037 

\/2& 

1.261+ 

11+6.9 

12.0 

0.0030 

3.02 

0.0921+ 

0.058 

\<2S 

1.155 

153.7 

10.1 

0.0021+ 

3.02 

0.0927 

0.167 

V!/2u 

0.0 

297.9 

1.7 

0.0011 

2.1+3 

0.11+28 

1.322 

D  (eV) 
o 


0.081+ 

O.03I+ 

0.01+9 

0.158 

1.301+ 
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Table  XVI.  Spectroscopic  Constants  for  Xe,. 


(cm  1)  01  X  (cm'1)  ct  (cm"1)  r  (a)  B  (cm"1)  D  (eV)  D  (eV 


(vertical  excitation  energy,  repulsive  curve) 
^.5  0.3  0.0002  3.97  0 
(vertical  excitation  energy,  repulsive  curve) 


8.5 

0.3 

0.0009 

1+.00 

9.7 

0.3 

0.0003 

3.91 

r.  2 


0.1+ 


0.0001 


■3.18 


0 


Table  XVII. 


luminary  of  Spectroscopic  Constants  for  the  A  V  State 
of  the  Noble  Gas  Dimer  Ions  1^2u 


Ion 

Dg(eV) 

{-SO} 

Dg(eV) 

<  r  >(a.u. ) 

{-SO} 

R  (a.u. ) 

Se2* 

1.35 

1.37 

0.965 

3.19 

V 

1.32 

l.Uo 

1.663 

^•59 

+ 

OJ 

1.19 

- 

1.39 

1.952 

5.17 

Xe2+ 

1.07 

1.38 

2.338 

5.9b 

{-SO}  :  without  spin-orbit  effects 
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Table  XVIII.  Dissociation  Energies  for  Electronic  States  of  the 
Noble  Gas  Dimer  Ions.  Energies  in  eV. 


State 

+ 

Ar^ 

+ 

Kr 

+ 

Xe.+ 

Present 

Study 

Other 

Work 

C- 

Present 

Study 

Other 

Work 

Present 

Study 

-d—— 

Other 

Work 

- 2 - 

Present  Other 
Study  Work 

n  2  + 

D  y 

l/2g 

repulsive 

repulsive 

repuls ive 

repulsive 

C  \,2» 

0.08 

0.08 

0.10e 

0.04f 

0.17 

0.l4e 

0.10f 

0.21  0.19e 

0.12f 

B  \/2g 

0.l4 

0.07b 

0.05 

repuls ive 

repulsive 

„  2 

C  h/2n 

0.05 

0.03 

0.01f 

0.07 

0.02f 

0.04  0.03f 

B  2ll3/2g 

0.20 

0.07b 

0.16 

0.l4e 

0.10f 

0.16 

o.i6e 

f 

0.13 

0.15  0.l8e 

0.12f 

A  2  l/2u 

1.31 

1.35a 

1.20b 

l.b5C 

1.10d 

1.30 

1.33g 

1.19f 

1.25° 

1.20h 

1.18 

1.17b1 

f 

1.05 

1.0b  1.03J 

0.79f 

T.  R.  Connor  and  M.  A.  Biondi,  Phys.  Rev.  l40,  A778  (1965); 

L.  Frommhold  and  M.  A.  Biondi,  ibid  185,  A244  (1969) . 

Ref.  7  (Cl  calculation) . 

Refs.  5,  b  (SCF  calculation). 

d  Ref.  lb* 

0 

Ref.  l4  (Photoionization  studies). 

f 

Ref.  11  (ci  calculation). 

S  Ref.  36* 

Ref.  10  (SCF  calculation) 

1  Ref.  4-3* 

J  Ref.  39* 

* 

This  reference  contains  a  summary  of  Dq  from  other  experimental  studies. 
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APPENDIX  E 


ELECTRONIC  STRUCTURE  OF  THE  NOBLE  GAS  DIMER  IONS. 

II.  ABSORPTION  SPECTRUM  FOR  THE  A  2E+  +  D  2l+  SYSTEM 

§  g 

ABSTRACT 

A  systematic  study  of  the  electronic  structure  and  chemical  binding  in 
the  dimer  ion  sequence,  Ne2+,  Ar2+,  Kr2+  and  Xe2+,  has  been  carried  out 
using  both  density  function  and  ab  init io  configuration  interaction  computa¬ 
tional  approaches.  From  the  results  of  this  study,  the  absorption  spectrum 

for  the  UV  transitions  of  the  AT  -*■  D  E  system  has  been  calculated. 

u  2 
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APPENDIX  F 


ELECTRONIC  STRUCTURE  OF  THE  NOBLE  GAS  DIMER  IONS. 

III.  ABSORPTION  SPECTRUM  FOR  THE  A  2I  +  -*■  B  ^7g  SYSTEM 

u 

ABSTRACT 

A  systematic  study  of  the  electronic  structure  and  chemical  binding  of 
the  dimer  ion  sequence,  Ne2+,  Ar2+,  Kr^  and  Xe2+,  has  been  carried  out  using 
both  density  functional  and  ab  initio  configuration  interaction  computational 
approaches.  From  the  results  of  this  study,  the  absorption  spectrum  for  the 
visible/lR  transition  of  the  A  Eu+  -*  B  2/7g  system  has  been  calculated. 
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APPENDIX  G 


ELECTRONIC  STRUCTURE  AND  PHOTOABSORPTION  OF  THE  He  +  DIMER  ION 

2 

ABSTRACT 

Calculations  have  been  performed  on  the  electronic  structure,  potential 
energy  curves,  spectroscopic  constants  and  photoabsorption  characteristics  of 
the  Hg^  dimer  ion.  These  studies  indicate  that  the  ground  A%u  state  of  the 
Hg^  ion  is  stable  with  a  calculated  dissociation  energy  of  0.44  eV.  This 
ion  exhibits  strong  absorption  (a~1.0  A^)  for  the  A^u  -*■  B^g  transition  at 
X  =  1.1  p .  The  only  other  predicted  absorption  is  for  wavelengths  shorter 
than  250  nm.  We  conclude  that  photoabsorption  of  Hg2+  does  not  appear  to 
represent  an  important  loss  mechanism  in  excimer  lasers  based  on  Hg  chemistry. 
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APPENDIX  H 


ELECTRONIC  STRUCTURE  AND  PHOTOABSORPTION  PROPERTIES  OF  NOBLE  GAS  TRIMER  IONS 


ABSTRACT 


A  study  of  the  electronic  structure  and  absorption  characteristics  of 
the  Ar^  ion  has  been  carried  out  to  determine  the  relative  importance  of  such 
species  in  the  analysis  of  loss  mechanisms  in  the  noble  gas-halide  excimer 
systems.  This  study  included  detailed  quantum  mechanical  calculations  of  the 
potential  energy  hypersurfaces  for  Ar^"*"  and  prediction  of  the  absorption  bands. 

We  find  this  ion  to  have  trigonal  symmetry  with  a  bound  E'  ground  state 

9  2  2 

+  Al  in  C2V^  *  The  vertical  excitation  spectra  are  k.2'  (.21  eV) , 

^E"  (.63  eV) ,  ^2" (1.06  eV)  and  ^A^ '  (2.12  eV) .  The  short  wavelength  tran¬ 
sitions  (300-500  nm)  are  found  to  be  weak  for  photoabsorption  from  the  ground 

state.  A  strong  photoabsorption  band  (\~590  nm)  is  predicted  for  the 
2  2 

E'  -*■  A^  transition.  This  suggests  that  trimer  ion  absorption  may  be  very 
important  in  the  analysis  of  loss  mechanisms  in  excimer  lasers  operating 
under  high  pressure  conditions. 
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APPENDIX  I 


ELECTRONIC  STRUCTURE  OF  NOBLE  GAS  DIMER  AND  TRIMER  IONS 


ABSTRACT 


A  systematic  study  of  the  electronic  structure  and  chemical  binding  of 
the  noble  gas  dimer  and  trimer  ions  has  been  carried  out.  This  study  includes 
detailed  quantum  mechanical  calculations  of  the  ground  and  low-lying  excited 
state  potential  energy  curves  for  Ne^,  Ar2+,  Kr2+  and  Xe2+  and  the  ground 
and  low-lying  excited  state  surfaces  for  Ar-^.  Both  density  functional  and 
ab  initio  configuration- interaction  computational  techniques  have  been  em¬ 
ployed.  A  regular  progression  is  found  in  the  calculated  spectroscopic 
properties  of  the  dimer  ions  which  results  in  some  remarkably  simple  con¬ 
clusions  regarding  the  nature  of  the  binding  and  the  geometry  of  these  ions. 
For  the  Ar^+  trimer  ion,  studies  in  D^,  C^v  and  D^  point  groups  were 
carried  out.  The  lowest  surface  corresponds  to  the  degenerate  ^E '  symmetry 

O 

with  an  indicated  bond  strength  of  ~0.2  eV  relative  to  Ar2+(A  I+u)  + 

ArC^Sg).  This  trimer  ion  exhibits  a  Jahn-Teller  distortion  arising  from 

the  degenerate  vibrational  modes  of  the  ground  state.  Several  low-lying 

.  2  9 

excited  state  surfaces  of  Ar^  have  also  been  studied  include  A2 ' ,  E"  and 

2 

A^'  symmetries.  A  strong  photoabsorption  band  (X~590  nm)  is  predicted  for 
2  o 

the  E'  ZA^  transition  of  this  trimer  ion.  This  suggests  that  such  species 
may  be  very  important  in  the  analysis  of  loss  mechanisms  in  excimer  lasers 
operating  under  high  pressure  conditions. 
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APPENDIX  J 


VISIBLE  PHOTOABSORPTION  BY  NOBLE  GAS  TRIMER  IONS 


ABSTRACT 


A  study  of  the  electronic  structure  and  photoabsorption  characteristics 
of  noble  gas  trimer  ions  has  been  carried  out  using  quantum  mechanical  methods. 
This  study  includes  detailed  density  functional  calculations  of  the  potential 
energy  surfaces  for  a  prototype  system,  Ar3+,  an  analysis  of  the  spectro¬ 
scopic  properties  of  the  ground  state  of  this  ion  and  prediction  of  the  strong 
absorption  bands.  Studies  in  D3h,  C3v  and  symmetries  were  carried  out. 

The  lowest  state  of  Ar3  corresponds  to  the  degenerate  ^E'  symmetry  with 
c.n  indicated  dissociation  energy  of  0.17  eV  relative  to  Ar„"*"  (A^Z  "*")  +  Ar  (^S  ). 
This  trimer  ion  exhibits  a  small  Jahn-Teller  distortion  which  does  not 
significantly  alter  our  conclusions  concerning  the  ground  state  stability  or 

absorption  characteristics.  A  strong  and  very  broad  photoabsorption  band 

2  2 

(1=496-670  nm)  is  predicted  for  the  E'  +  Aj  transition  of  this  trimer  ion. 

This  suggests  that  such  species  may  be  very  important  in  the  analysis  of 
loss  mechanisms  in  visible  excimer  lasers,  operating  under  high  pressure 
conditions,  with  mixtures  containing  the  noble  gases. 
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